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ABSTRACT 


Assuming that hyperfine structure in spectral lines has its origin in the coupling 
of a nuclear spin with electron resultant J, a comparison between gross structure and 
hyperfine structure is made from the following relation, Av,/Avy=m,/4im,. For 
those atomic systems for which the electron configurations show LS or jj coupling, 
the agreement with the above equation is quite as good as would be expected. For 
electronic configurations which do not show LS or jj coupling it is doubtful whether 
such a comparison of gross structure with hyperfine structure can be suitably made. 
From the meagre data as yet available, one is led to suspect that nuclear spin is 
due to the last electrons, protons, or both (most probably electrons) added to the 
nucleus in atomic construction. These electrons should, in accordance with the 
principles of the wave-mechanics, have a probability density spherically symmetrical 
about the nucleus. The difference between the electrons associated with the nucleus 
and the outer electrons is that the nuclear electrons carry with them in space quantiza- 
tion the total mass of the atom. 


HE splitting up of ordinary energy levels into a number of hyperfine 

structure components is quite generally attributed to a nuclear angular 
momentum. Observations indicate that this nuclear spin, whether due to 
electrons, protons, or to both electrons and protons, carries with it in space 
quantization the total mass of the nucleus. Assuming the nuclear spin to be a 
rotation of the nucleus as a whole, Jackson? has derived the following equa- 
tion for the approximate shift in energy of ordinary energy levels due to 
nuclear spin. 


Av,/ Avy =m,/2m.. (1) 


Here Ay, is a gross or multiplet structure separation, Av; the corresponding 
hyperfine structure separation, m, the mass of an electron and m, the mass of 
the nucleus. Since this formula was derived for a nuclear spin of i=} a 
quantum of energy, it may be made general for all values of 7 by replacing the 
constant 2 in the above equation with 47. This gives 


1 National Research Fellow. 
2 Jackson, Proc. Roy. Soc. A121, 432 (1928). 
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Av ,/ Avy =m;/ tim (2) 


which for 7 =} reduces to equation (1). 

Jackson has shown® that in the simple case of a so called one electron 
system, like caesium, an electron spin doublet separation (Av,), for example 
2 P\2—°P3.(6p), is to be compared with the hyperfine separation (Av,) of one 
_of the terms, ?P. The ratio Av, /Av;, as compared with computed values from 
Eqs. (1) and (2), isin reality in quite good agreement (Table I). 

From a discussion with Dr. Schiiler* it appears that the hyperfine struc- 
ture of \5486, °S,—*P 1.2, in singly ionized lithium may be accounted for by 
assigning an angular momentum, 7=}3, to the nucleus of Li.? From the ob- 
served hyperfine structure, the *°S,(1s2s) term and the *P.(1s2p) term have 
normal doublet separations of about 0.60 cm! and 0.45 cm~' respectively 
(error +0.10 cm~!). The splitting up of these terms is shown in Fig. 1. In 
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Fig. 1. Fig. 2. 


this case of a nucleus and two electrons, it may be shown by the use of Fig. 2 
which gross structure and which hyperfine structure separations are to be com- 
pared with each other. The left hand vector 7, in Fig. 2, may be taken to 
represent the nucleur angular momentum, 3, of Li,’ and 7; and je may be 
taken to represent 1s and 2 electrons respectively. In the space quantiza- 
tion of these three vectors the coupling between 7; and jz produces the gross- 
structure energy levels, whereas the coupling between 7; and 7 and between 
je and i produces the hyperfine energy levels. For the term 5/.°P%: of Fig. 1, 
all three vectors 7, j; and jo, Fig. 2(a), are parallel as shown.‘ When 7 is now 
inverted, to position Fig. 2(b), the shift in energy to the term 3/2°P% is ob- 
served to be about 0.45 cm~!. Owing to the relative nearness of the 1s electron 


3 Schiiler, Zeits. f. Physik. 

‘With the apparent need of the new quantum number f to represent an energy level, 
it is suggested that this quantum number preceed the ordinary term designation as a sub- 
script for example, 3/2°D°,. If the electron configuration is given, it is desirable to include the 
nuclear spin, for example, is;24s3d. A radiated hyperfine structure line would be designated, 
372° Si 5/2°P°2. 
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to the nucleus, this shift in energy is almost entirely due to the coupling be- 
tween j, and 7. When from 5,.°P%s, Fig. 2(a), the vector j; is inverted to position 
Fig. 2(c), the shift in energy to 3;.!P; is observed to be about — 7000 cm-. 
The inversion now of 7 to position 1,2'P,, Fig. 2(d), should produce, although 
not yet observed, a shift of about —0.4 cm~'. The inverted hyperfine struc- 
ture in this 'P, term follows from the hypothesis that, while the separation 
arises from the coupling between j; and 7, space quantization must take place 
between nuclear spin, 7, and electron resultant, J.6 The ratio Av,/Av; is given 
in Table I along with the values obtained for other elements. 

Making quite the same arguments given above for lithium, a correlation 
between gross structure and hyperfine structure may be made for an atomic 
system having a number of valence electrons, for example, manganese, Mn I. 
Using again Fig. 2, vector 7 is taken to represent the nuclear angular momen- 
tum 7=5/2 of manganese, vector j; to represent a 4s electron and, vector j2 to 
represent the resultant angular momentum of all of the rest of the electrons 
outside of the nucleus, for example 3d°, Ds, of Mn II. With all three vectors 
parallel as in Fig. 2(a) the atom is in state 7°Do9/2(t5,2.3d°4s)._ The inversion of 
vector 7 to position (b), 2°Dg/2(%;,23d®4s), produces a shift in energy of about 
0.50 cm~!. An inversion, however of vector j,; from position (a) to position 
(c), 6£Dz/2(is,23d°4s), produces a shift in energy of —6,244 cm~'. The ratio 
Av, /Av;= 12,488 is in quite good agreement with the calculated value 10,100. 

In order to account for the observed hyperfine structure of the principal 
doublet of caesium, Jackson? assigned an angular momentum of i=} to the 
nucleus of caesium. It is interesting to note that the observations are equally 
well explained by assuming greater values of 7. From values of nuclear spin 
found in other elements 1/2, 5/2 and 9/2, and the last two columns of Table 
I, excluding Tl and Cd, the value i = 5/2 seems to be, although the evidence is 
not strong, the more probable for caesium. 
































TABLE I. 
A comparison of gross structure with hyperfine structure. 

Element t Poy Terms Av, Avy Av, /Avy m,/4im, 
Li’ Il i 1s2p sP~—'P,° 7,000 0.45 15,555 6,440 
Mn I 23 3d*4s ‘Dy —*Dy 6,244 .50 12,488 10, 100 
Mn I 2} 3d*4p4s 8. S44° — 6 Sy1° 9,984 .60 16,640 10,100 
Cs I 3? 6p 2P;°—*P,)° 554 01 55,400 122,000 
Cs I 23? 6p 2P\°—?P,\° 554 01 55,400 24,500 
Cs I | 43? 6p 2Py°—*P\° 554 01 55,400 13,600 
Tl I 3 6p 2P\°—2P,\° 7,793 .40 19,480 187 ,600 
Tl Il 3 6s0p sP,°—'P,° 13,935 3.5 3,980 187 ,600 
Cd I $ 5s5p sPP—'P,° 11,864 .25 47,400 103 ,000 








The ratios Av,/Av,;, given in Table I for Tl I and TI II, are decidedly not in 
agreement with the calculated value 187,600. From Zeeman effect data on 
43776 of Tl I® the nuclear spin is quite certainly 1/2. The enormous hyperfine 


5H. E. White, Phys. Rev. 34, 1404,(1929). 
6 Schiiller and Briick, Zeits. f. Physik 56, 291 (1929). 
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structure separations in Tl II as compared with the small separations found 
in Tl I are undoubtedly due to the strong coupling between 7 and the deep 
penetrating 6s electron for Tl II, in contrast with the coupling and penetra- 
tion of a 6pelectron in TI 1. In TI II, it would appear that the difficulty arises 
in the uncertainty of the coupling between the two valence electrons. Since 
the terms do not indicate either LS or jj coupling, it is doubtful whether such 
a comparison of gross structure with hyperfine structure can be suitably made 
for such cases. It has been interesting to note, however, for all observed hy- 
perfine structure, that the space quantization between electron resultant J 
and the nuclear spin 7 always follows the interval rule,i.e., a/7 coupling where 
the intervals are proportional to cos Ji. The exact form of the constant 
multiplier of 7 in equation (2) depends largely upon the distribution of charge 
causing the nuclear spin. Even with the meagre data as yet available, one is 
led to suspect that the nuclear spin is due to the last added electrons, protons, 
or both electrons and protons. From observations electrons alone seem to be 
the most probable.’ These electrons should, in accordance with the principles 
of the wave mechanics, have a probability density spherically symmetrical 
about the nucleus and might be thought of as a charged shell over the surface 
of the nucleus. The difference between the electrons associated with the 
nucleus and the outer electrons is that the nuclear electrons carry with them 
in space quantization the total mass of the atom. 

The author wishes to take this opportunity to thank Dr. Paschen for his 
interest in and criticism of this paper. 
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THE SECOND SPARK SPECTRUM OF ANTIMONY AND A NOTE 
ON THE FIRST SPARK SPECTRUM OF TIN 


By R. J. LanG 
UNIVERSITY OF ALBERTA, EDMONTON, CANADA 


(Received January 30, 1930) 


ABSTRACT 
Some sixty lines of the vacuum spark spectrum of antimony have been classified 
as transitions between the following terms of Sb III: (s*5p)*P, (s*6p)?P, (s*6s)*S, 
(s*7s)?S, (s*8s)?S, (s*5d)*D, (s°6d)?D, (s*4f)*F, (s*5f)*F, (s*5g)*G, (s*6g)?G, (s5p*)*P, 2D, 
*S,*P and (5p*)4S. Corresponding to the deepest term value an ionization potential 
of approximately 24.7 volts is calculated. A note on SN II gives the terms *P and *S 


from the (s5p*) configuration and new value for the terms: (s*7s)*S, (s*8s)*S, (s*7d)?D 
and (s*6f)?F. 


HE vacuum spark spectrum of antimony has been measured over the 

range 6600 to 600A and also the Paschen hollow-cathode spectrum from 
5200 to 1950A. By means of these two measurements it has been possible 
to separate the first and second spark spectra rather definitely from one an- 
other and thus to arrive at a classification for the latter. Work on the first 
spark spectrum is now in progress and will be reported on later. The measure- 
ments of the vacuum spark spectrum below 2500A were made on a two- 
meter grating mounted in a vacuum spectrograph giving a dispersion of 
4.5A per mm. Above 2000A several instruments were used. The vacuum 
spark and hollow cathode spectra were measured in Edmonton on a two- 
meter Rowland grating kindly loaned to the author by Professor Smith. The 
vacuum spark spectrum was also measured in Michigan by the author using 
a Hilger &, quartz spectrograph made available through the courtesy of the 
Department of Physics. The author is also indebted to Professor R. A. 
Sawyer of Michigan for a plate of the vacuum spark in the region from 5300 
to 6600A made on a large Hilger glass spectrograph. A few wave-lengths are 
taken also from Schippers measures given by Kayser. 

In a previous report! the author gave three multiplets for Sb III all of 
which turn out to be real but the first one given was improperly classified as 
(s*6s)?S—(s*6p)?P. It should have been (s5p?)?S—(s*6p)?P. It thus corre- 
sponds very satisfactorily to the doublet in As III at 4200A. A paper on AsIII 
and Sb III appeared later in the Indian Journal of Physics.? 

Some important electron configurations and resultant terms expected in 
the spectrum of Sb III are as follows: 


' Lang, Phys. Rev. 32, 737 (1928). 
* Pattabhiramiah and Rao, Indian Journ. Phys. 3, 437, 1929 
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Configuration Terms Configuration Terms 
(s*5p) 3p (s5p*) ‘P,2?PDS 
(s°6s) 28 (s5p5d) ‘PDF,*PDF 
(s*5d) 2D (s5p6s) ‘p2pP 
(s24f) °F (5p) ‘S, 2PD 
(s*5g) 2G; (s5 pop) 4PSD,?SPD 
TABLE I. Empirical term values in Sb III. 
Odd Terms | Even Terms 
(s°5p)2P} 200272 | (sSp*yDy 123744 
6576 | 1270 
(sS5p)*P 193696 | (s5p*)*Day 122474 
(s*4f)? Foy 118940 | (s*6s)?Sy 107321 
373 
(s*4f)? F3) 118567 (s5p*)?Sy 106852 
(s°6p)?P3 85550 (s5p?)*P, 105624 
1668 5392 
(6p)? Pi; 83882 (sSp*)?P iy 100232 
(s°5f)? F3, 64055 (s*5d)*Dyy 101450 
—55 1567 
(s*5f)? Fo: 64000 (s*5d)?D21 99883 
(5p*)4Sy; 51259 (s°7s)2Sy 57143 
i ditinsinaisiinaatbicatin = (s°6d)*Dy 55589 - 
Even Terms — 237 
a be (s°6d)*D24 §5352 
(s5p*)*P} 145906 (s*8s)?Sy 36882 
3596 
(s5p?)*P, 142310 (s°5¢)°G 35966 
5358 
(s5p*)*P 2 136952 (s*6g)°G 25812 
TABLE II. Comparison of term values. 
Term Element 4 5 6 7 8 
In I 22295 10366 6031 
2S Sn Il 15205 7841 5194 
Sb Ill 11924 6349 4098 
In I 46668 14811 7808 
*Py Sn Il 29426 11553 
Sb Ill 22252 9506 
In I 13775 7620 4832 
*Dy Sn Il 11575 6867 4342 
Sb Ill 11272 6176 
In I 6960 
2 Fo Sn Il 7104 4511 
Sb Ill 13215 7111 
In I 
2G Sn Il 
Sb III 3996 2868 
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TABLE III. Combinations, intensities and discrepancies in the spectrum of Sb III. 











Odd Terms || 





= x) * a x) a = 
oy <3 ad i - Q = x Y 
— a ta] e -_~ i] nN = 
a » = . x -_~ —_ _~ = 
nN a — ~ na aA —~ ~ 
“ wn + be a oOo =) nN wy A 
o iS) a a os ° wy 
| “ %, “ “ “ “ “ ~~ 
~~" — — ~ ~~ 
—- — 





— | 200272) 193696 | 118940 |118567| 85550| 83882 | 64055 | 64000 | 51259 




























































































; i} 12 2 | 10 
(s5p?)*P; 145906 || 2 | 1 | | 9 
| 15 | 10 | | | | 10 
(s5p°)"Pyy 142310 |) —2 | —2 | | | a 
ee |e ee ee es es ee “2D 
(s5p2)*Px 136952 | | 0 | | | | | 0 
| 20] wo} | |a2i] 5] is | 38) 
(s5p*)*Dy 123744 || 1 | 1 | | 5 | 3 i wt ff 
| | 20° } ft 45 | 45 | 3 |) (208 
(s5p?)*Doy 122474 || —6 | | | § | -2°| -3 | 1 
| 30 40 | | 30 | 40 | 
(s%s)*Sj 107321 | 1 | 0 | | of} o| | 
| 2 | 30 | | | 30 530 | 
(s5p*)*Sy 106852 || 0 | 0 | | 0 0 | 
| 10 | 10 | | | 
(s5p*)*Py 105624 || —3 0 | | 
/ 4 | 20 | «| 31 3/ a | | wf] 
(s°5d)2Dy3 101450 | —2 oi «<4 | -15n @ | =! 0) 
| 15 | 20 rf | | | 
(s5p*?)?P14 100232 | —2 0 
| 40 1; | | s{]2of}] s5s|- 
(s*5d)*D2, 99883 || 4 —6 | | 3 | 0 1 | 
10 15 | 15 30 | | 
(S7s)*Sy_ 57143 || -3 —3 | 1 =i | | | 
| 10 8 2 | 15 | 
(s%6d)2Dy «85589 || —3 | —2 | 1 0 | 
| 15 | 50 
(s*6d)*Dxy 55352 | —4 | 0 | 
4  f +2) £f¢ f¢ 
(s*8s)2Sy 36882 0 0 | | 
50 | 50 | 40 40 | 
(s?5g)°G 35966 0 0 | —1 —1 | 
| 5 5 | r et £7 
(s*6g)?G 25812 | 3 0 | | —2 | 3 | 
_ a 














* Classified twice. 
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TABLE IV. Classified lines of Sb III. 





Cc lassification ACL. \. ) | v 








Pete i SP lle a aa ) 
(s*5p)*P 4— (s*6d)* Dy 91 ; 18 10 144680 | 
-Pyy—  * Dy 722.86 15 138340 
2Py = *Dy 724.81 8 138105 
(s*5p)?P 3 — (#75) st 698 .69 10 143126 
*Pyy— 732.33 15 136550 
(s*5p)*P 4 — (s5p*)? Pu 999 62 15 ~ 100038 | 
PF 1056.58 10 94645 
“Py - Ph 1069 .93 20 93464 
*?Py— *P3 1135.43 10 88072 
(S'S PRP 4 —(s"5d)*Dy 1011.94 40 98820 
*Pyy— 2D 1065 .90 40 03817 
Py, 2Diy 1084 .06 20 92246 
(s?5p)2P 4 — (s5p2)2S) 1070.43 20 93420 
*Piyy— 2S) 1151.49 30 80844 ) 
(s5p*)'P — (5p°)'Sy 1056.58 10 94645 
‘Pi - ‘Si 1098 . 34 10 91047 
‘Pi - ‘S11 1166.96 12 85693 
(s*4f)? Foy — (s°6g)°G 1073.76 5 93131 
Fy—- 1078.10 5 92735 
(s*5p)?P 4 — (s*6s)*Sy 1075.82 30 92952 
2Py— 8S 1157.74 40 86375 
(s*4f)? Foy — (s*5g)°G 1205.20 50 82974 
* F334 — 2G 1210.64 50 82001 
(s°5p)2P 4 —(s5p°)*Dy 1306.69 20 76529 
*Pu— 2D 1404.18 20 71216 
2?Py- = "Dy 1429.57 10 69951 | 
(s5p2)2Dy — (5p°)*Sy 1379.58 8 72486 
2D — 4S13 1404.18 20 71216 ) 
(s5p*)2Duy — (S*5f)?Foy 1673.89 15 59741 
2Dy— °F 1710.23 3 58472 ! 
“Dy — *F 3 1711.84 15 58417 i 
(s*5p)2°P 4 —(s5p°)*Py 1725.33 15 57960 
*Piu— ‘Po, 1762.30 12 56744 
2Py— ‘P; 1839 .32 12 54368 
2Py- “Py 1946.13 10 51384 
Py ‘Py 2091.85 2 47791 
(s*6p)*P 4 — (s°8s)?.Sy 2054 .10 3 486068 
2Py- 8S 2127.00 5 47000 
(s5p2)*Diy— (s*6p)*Py 2507.71 5 39865 
2D — *Pr 2590.13 15 38597 
2Du- Py 2617.17 12 38199 
(s*5f)? Fa, — (s?6g)°G 2614.20 1 38241 
* Fy — 2G 2617 .63 1 38191 
(s*5d)*Dyy — (s*5f)? Fay 2669.39 20 37450 
2Dy— *Fy 2785 .87 5 35884 
2D — 2Fy 2790.27 20 35828 














———— 
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TABLE IV (Continued) 


Classification A(1.A.) I v 
(6p)?P 5 —(s°6d)*?Dy 3336.61 20 29962 
*Pyy— 2Dox 3504 .07 50 28530 
*Py— “Diy 3533.45 15 28293 
(s*op)*?P 4 — (s°7s)*S, 3519 .06 15 28408 
*Py— 2S) 3738.90 30 26738 
(s*5f)? F33 — (s*5g)°G 3559.18 40 28088 
2Fy — °G 3566.25 40 28033 
(s*6s)?S 3 —(S6p)?Pi3 4265 .09 40 23439 
°Sy— *Py 4591.89 30 21771 
(sSp?y?S4— (Sop) Py 4352.16 50 22970 
2S *Ps 4692.91 30 21302 
(s*4f)* Foy — (s°5d)* Dy 5247.71 1 19051 
2Fy — *Diy 5717.3 1 17486 
*iy4— “Dy 5845.5 17102 
(s*8d)°Dys — (Op)? Py 5690.8 1 17567 
‘y= *Pis 6246.7 5 16004 
*?Du-— Pp, 6287.6 3 15900 


In Table I the terms which have been located together with the empirical 
term values and separations are given. All of the term values rest upon an 
arbitrary choice of 64000 cm~ for the value of (s*5f,?F4. Experience has 
shown that one can probably arrive at as good or even better values for the 
terms by a comparison such as is given in Table II than by the use of a 
Rydberg formula applied to two P or two S terms such as are available in 
this instance. The deepest F terms are not absolutely certain since the 
combination with the deepest normal D terms is not very satisfactory. The 
same may also be true of the second G terms and the third S terms, the last 
depending upon one combination only since no others fall within the range 
of measurements. Hence until some of these terms can be established more 
certainly it seemed as well to estimate the value of the basic term 

In Table III the intensities of observed combinations are shown and under 
each intensity the discrepancy between the observed wave-number and the 
wave-number calculated from the values assigned to the energy levels is 
given. In Table IV a list of all the lines which have been classified in the 
spectrum of Sb III is given. Throughout this report all wave-lengths above 
2000A are given in I.A. (air) below this value in I.A. (vacuum) while all 
wave-numbers are reduced to vacuum values. Corresponding to the deepest 
term of this spectrum an ionization potential of 24.7 volts is found. 


A Norte on SN II 


The hollow cathode spectrum of tin was also photographed between 
6000A and 2000A in an endeavor to find more terms from the (s5p?) con- 
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figuration. Out of the total of nine terms (*P, *D, 2S, ?P) arising from this 
configuration none but 7D had been located.* 

The (s*5p)?*P—(s5p*)*P multiplet in In I had already been found* and 
those for Sn Il and Sb III were located by the use of the doublet laws. The 
multiplet in Sn II is taken from the spectrum of the hollow cathode. It occurs 
but very weakly in the vacuum spark spectrum. That for Sb IIT is taken 
from the vacuum spark spectrum and the only abnormality is the weak 
intensity of one of the lines in this multiplet. 


TABLE V. Multiplets and terms in the spectrum of Sn II. 


Configuration ' = X(T.A.) I | v Av | Term values 

(s°5p)?P 3—(s°7s)*S) 1158.21 2 | 86340 (@5pyPy 117704 
2Py— 283 1218.19 4 | 82089 4251 °Pyy = 113451" 

(s?6p)2Py — (s°8s)2S3 3030.37 -2-«|-« 25435 (st7s)28y 31363 
2Py1— 28) 4071.79 2 | 24552 883 

(s*5d)*Dyy — (s*6f)° F 3537.56 1 28260 
2D — °F 3620.54 1 27613 «O47 «|| (s°8s)2Sy 20778 

} 

(s5p)?Piy — (s5p°?)2 Diy 1699.50 = 12 58841 4252 || (s*7d)2Dy, 173067 
Py— 2D 1811.23 15 | 55211 622* || 2D», =—-:17327 
2p) — Diy 1831.85 12 | $4589 | 

| 

(s5p2)2Dyy — (s26f)2F 2448.93 1 | 40821 | (s6f)2F 18043 

2D — °F 2486.95 2 | 40197 624 | 
| (sSp*?)*P 3 73196 

(s*Sp)?P 1—(s5p*)*Piy 2151.57 12 46462 ‘Py 71242 
2Py\— ‘Py 2209.70 10 | 45241 1954 ‘Pr, 68210 
2P y— iP 2246.11 10 | 44508 
2Py— ‘Py 2308.26 15 | 42211 1957 (sSp°)2Dy, ——-§8863 
2Pyi— ‘Py 2483.48 12 | 40254 2D», ~—-58240* 

| 

(s*5p)2P y—(s5p2)2Sy | :1223.73 10 | 81717 | 
2Py— *Sy | 1290.89 20 | 77465 4252 || (s5p)2Sy 35986 

' 1] 

(s%6p)2P 4—-(s*7d)Dy, |S «3465.73 = 3«|:s« 28846884 | 

2Pyy— 2D} 3570.09 8 | 28002 | 
3575.31 4 


2Py— 2D: 27962 «40 

















* Green and Loring. 


The (s°5p)?P—(s5p?)?D multiplet seems well established in Sn II and 
Sb III. One peculiarity worth mentioning is that the position of this multi- 
plet in the spectrum of Sn II is almost exactly the same as in Si II°. However, 
there seems little doubt that in Sn II this group has been correctly classified® 
for when a tin spark in nitrogen, excited by a 2200 volt transformer joined 
to the 110 volt A.C. mains, was photographed between 2000A and 1250A but 
nine tin lines were found including the three lines of this multiplet, and every 
one of these lines except one (1290.90A) had been classified as arising from 
transitions involving the deepest ?P levels of the atom of Sn IIT. 


3 Green and Loring, Phys. Rev. 30, 574 (1927). 
4 Lansing, Phys. Rev. 34, 597 (1929); Sawyer and Lang, Phys. Rev. 34, 712 (1929). 
5 Fowler, Trans, Roy. Soc. A225, 26 (1925). 
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In regard to the (s*5p)*P—(s5p?)?S doublet it may be stated that the line 
(1290.90A) together with another line at 1223.79A which occurs in the 
vacuum spark, but unable here to pass the fluorite window, have the separa- 
tion of the deepest *P terms. There seems little doubt that these two lines 
represent this doublet as has been suggested.® 

The remaining (s*5p)?P—(s5p?)*P multiplet in Sn II has not been 
located with any certainty. There is one possibility worth mentioning, 
however. The (s*5p)?P—(s?7s)?S doublet® occurs with too great intensity in 
the vacuum spark and the 2S term is greater than one would expect. Another 
pair of weaker lines, one of which had not been resolved from a stronger tin 
line, probably represent this doublet and the stronger pair represents one-half 
of the missing multiplet. However, it must be said that there is not complete 
agreement yet between Sn II and Sb III in regard to these terms of the (s5p*) 
configuration and it does not seem possible to obtain this at the present 
time. When the arc spectrum of indium has been studied below 2000A it 
should be possible to obtain this agreement. An attempt to do this is now in 
progress. One or two other minor alterations have been made in the spectrum 
of Sn II all of which are shown in Table V. 

The author takes this opportunity of thanking the Research Council of 
Canada for a grant to carry on this research. 
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THE EMISSION SPECTRUM OF BENZENE 
IN THE REGION 2500-3000A 


By J. B. Austin* anp IAN ARMSTRONG BLACKT 
DEPARTMENT OF CHEMISTRY, YALE UNIVERSITY 


(Received January 25, 1930) 


ABSTRACT 
The ultra-violet emission band spectrum of benzene vapor has been observed with 
a Tesla discharge as a means of excitation. The wave-lengths of over 100 bands 
lying between 2400 and 3000A have been measured with an accuracy of +0.5A. 
It isfound that the measurements are in general in agreement with those of McVicker, 
Marsh and Stewart. All the bands except the fainest and most diffuse ones may be 
represented by the following equations: 
1 \ = 37485 + 924m — 161n 
1\ = 37398 + 924m — 162n 
1/A = 37485 — 986m — 162n 
where m and n are integers. 
The evidence supporting the conclusion that the bands are due to benzene and 
not to some decomposition product is summarized. The applicability of the data to 
spectroscopic analysis is also discussed. 


= 37425 — 986m — 162n 
36479 — 991m — 162n 
= 36412 — 991m — 162n 


—_—_ —_— — 
a 
II 


INTRODUCTION 


HE general features of the spectra obtained from the vapor of benzene 

and its simple derivatives under the excitation of the Tesla discharge 
have been described by McVicker, Marsh and Stewart,! who have designated 
them as Tesla luminescence spectra to distinguish them from the absorption 
and fluorescence spectra. The present research was undertaken as an ex- 
tension of their work; it was hoped, on the purely chemical side, to develop 
the suggestion made by them that a means of ascertaining the identity and 
quantity of any impurity present could be elaborated. In this regard it has 
been shown that such a method would be unwieldy in practice, except per- 
haps for routine work on a large number of samples; moreover, the number 
of likely impurities which yield satisfactory spectra under these conditions 
is rather small. On the other hand, it was felt that more comprehensive data 
on the Tesla luminescence spectra, particularly on the spectrum of benzene 
which is noteworthy for the sharpness of its bands and the absence of any 
continuous background, might be obtained, which with the increased ac- 
curacy would lead to some theoretical interpretation of the observed regu- 
larities. In this object some success has been achieved and the present work 
describes the results for benzene; the results for the simple derivatives will 
be published in the near future. 


* From a dissertation presented by J. B. Austin to the Graduate School of Yale University, 
June 1928, in candidacy for the degree of Doctor of Philosophy. 

+ Commonwealth Fellow 1926-1928. 

1 McVicker, Marsh and Stewart, J. Chem. Soc. 123, 642 (1923). 
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EXPERIMENTAL METHOD 


The methods employed differed from those of McVicker, Marsh, and 
Stewart in many ways, the most significant difference being the use of much 
lower pressures of the vapor. They used benzene vapor at pressures as high 
as 65 mm of Hg and were forced to use a steam jacket to prevent conden- 
sation; in the present work the pressures ranged from 0.01 to 0.10 mm which 
yielded satisfactory results at room temperature. In its final form our 
apparatus was fitted with internal electrodes whereas the previous practice 
used external ones. The general scheme for producing and maintaining these 
low pressures was to place the liquid benzene in a side limb of the discharge 
tube, to immerse it in liquid air to freeze the liquid to a solid cake, and then 
to seal the side limb and evacuate the system with a mercury diffusion pump 
until the residual air was so attenuated that a discharge could not be pro- 
duced. The liquid air bath was then partially removed to allow the frozen 
benzene to warm slightly so that its vapor at low pressure could continuously 
diffuse through the discharge tube to a liquid air trap where it was condensed. 
In this manner equilibrium between speed of pumping and rate of sublima- 
tion was soon established and could be modified to give pressures between 
the limits of 0.01 and 0.10 mm. This method of maintaining a constant flow 
of vapor in the tube helped to remove any decomposition products as soon 
as they were formed and in conjunction with the low kinetic energy, and 
consequent small molecular disruptive power of the Tesla discharge gave 
some assurance that the spectrum originated in the benzene molecule. As a 
further precaution against contamination a new tube and liquid air trap 
were used for each exposure. 

The first attempt to obtain the desired spectrum was made with the 
so-called “electrodeless” discharge, produced by placing the discharge tube 
in a helix of stout brass rod through which the current from a 1 KW Thordar- 
son transformer was passed. A few trials sufficed to show that the spectrum 
obtained under these conditions was not that of the undecomposed molecule, 
but was due apparently to molecular fragments. Moreover, this discharge 
was very sensitive to slight fluctuations in pressure; with the pressure above 
a certain limiting value the free electrons and ions did not seem to attain a 
velocity sufficiently high to excite the rest of the vapor; with pressures too 
low the intensity of the glow suffered appreciably. At times this discharge 
showed a tendency to form a shell adjacent to the inside wall of the tube 
which electrically shielded the interior and caused a luminescent cylinder to 
form; the dark core of this cylinder rendered the discharge quite unsuitable 
for photographing. 

It is interesting to note that the exciting coil and discharge tube grew 
quite warm during the passage of the current through the former; at the end 
of a run the latter invariably cracked, the cracks following the helix around 
the tube; whether the rupture of the glass was caused by thermal or electrical 
strain was not determined. 

The next attempt was made with electrodes of lead foil wrapped around 
a Pyrex tube. This set-up yielded a satisfactory spectrum but endless trouble 
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was experienced with the breakdown of the glass along the surface between 
the electrodes. The presence of air streaks in the glass made trouble certain. 

Finally it was decided to install internal electrodes, a procedure which 
had hitherto been regarded as hazardous because of the possibility of 
excessive decomposition of the vapor. Such electrodes, however, proved the 











Fig. 1. Discharge tube. 


most satisfactory of all the arrangements tried and in most cases caused no 
greater decomposition than the other types. This system was used to 
obtain all the results herein reported. Fig. 1 shows the discharge tube which 
consisted of a cylinder of Pyrex, P, with a quartz window, VW, sealed on at 
one end with de Khotinsky cement; the other end was connected to a vacuum 


0.003uf 


{ kw -20,000v. 
Thordarson transformer | - | Tesla coil 


Fig. 2. Electrical connections. 





| to electrodes | 





system and a mercury diffusion pump. The electrodes, 7, were cylinders of 
aluminum attached to tungsten seals; a side bulb, S, held the benzene. 

The electrical system is shown in Fig. 2. When properly tuned this circuit 
gave a potential of about 120,000 volts as estimated from the spark produced 
between metal points. 
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SPECTROSCOPIC AND PHOTOGRAPHIC METHODS 


In this work two spectrographs were employed: 

(a) A Hilger size C quartz spectrograph, giving a photograph 200 mm in 
length for the range 2100A to 7000A. On this instrument exposures ranging 
from 3 to 12 hours were made; a plate exposed for 12 hours was used for the 
final measurements. For comparison, the spectrum of a mercury arc was 
superposed on the middle of the benzene spectrum with the aid of a reducing 
stop on the slit. 

(b) A Hilger E, quartz spectrograph (Littrow mounting) giving a spec- 
trum 140 mm in length for the region 2500A to 3000A. The spectrum of the 
mercury arc was superposed on these plates as on the others. 

All photographs were taken on Eastman 33 plates of a special thin glass 
so that they could conform to the curvature of the plate holder without 
danger of breaking. In the earlier stages of the work Nujol and transfor- 
mer oil were tried as sensitizers but were later discarded since they did not 
appreciably enhance the intensities of the ultra violet lines and were not 
deemed of sufficient value to warrant the extra trouble. 

The spectra obtained with the smaller instrument extended almost 
150A farther into the ultra-violet than those from the larger one; on the 
latter no bands were measured beyond 2600A while on the former they 
were recorded as far as 2471A. This difference is probably due to the greater 
absorption of the thicker quartz prism in the big spectrograph. On the other 
hand the large Hilger resolved some of the bands and permitted greater 
accuracy in measurement; hence, to obtain a complete spectrum of con- 
siderable accuracy the data obtained with the two instruments have been 
combined. 


WAVE-LENGTH MEASUREMENTS 


All plates were measured on a Société Genevoise ruling engine whose 
ruling mechanism had been replaced by a plate holder and a fixed telescope 
with cross hair attachment. The screw head was graduated in divisions of 
0.005 mm with a vernier which allowed readings to be made to 0.001 mm. 
Tests on the accuracy of the screw showed that any deviations or flaws in 
it were negligible compared to other errors. The wave-lengths of the unknown 
lines were calculated by means of the parabolic interpolation formula: 


Av=Aotas; +s? 


where Xo, @ and b are constants determined from the mercury lines? and s; 
is the reading of the measuring engine corresponding to the i** line. This 
formula proved to be as accurate as the observations for ranges of 200A in 
the lower visible and ultra-violet and involved much less labor than the 
Hartmann formula. The whole spectrum was covered in overlapping 
sections of 200A. 

Each plate was measured in the center of the spectrum to avoid any error 
due to curvature of the lines. The further precaution was taken of having 


2 Kayser, Tabelle der Hauptlinien der Linienspektra J. Springer, Berlin 1926. 
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each of us independently measure each plate twice, once when it was run 
along the measuring engine forward and once backward; in this way each line 
was approached from both sides. The agreement obtained in these sets of 
observations was most satisfactory, the difference rarely exceeding +0.5A 
and in most cases being less than +0.25A. The four measurements of each 
line have been averaged and converted to wave numbers in vacuo. 


PREPARATION OF BENZENE 


Merck’s product was shaken six times with concentrated sulfuric acid; at 
the end of this treatment no color was imparted to fresh acid. The benzene 
was then washed successively with concentrated sodium hydroxide, dilute 
sodium hydroxide, and distilled water, until the final wash water was neutral 
to litmus. The resulting product was dried over sodium for 25 hours and 
redistilled from the same. A large fraction came over at a constant tempera- 
ture and this was taken for use. B. P. 80.0°C at 760 mm Hg pressure. 


RESULTS 


The general scheme of the bands is shown in the photograph, Fig. 3. The 
following table gives the average results obtained from the two spectrographs. 
Lines observed only on the smaller spectrograph are marked with an as- 
terisk. The results of McVicker, Marsh, and Stewart are given for compari- 
son. 

While a rough notion of the relative intensities of the bands may be 
gathered from Fig. 3, a more exact measure for the most prominent band 
groups is given in the microphotograph, Fig. 4. This was made on a Koch- 
Goos microphotometer from a plate taken on the large spectrograph. The 
ratio of length on this curve to length on the plate is 2: 1. 

It was found that all the principal bands could be accounted for with an 
accuracy of +3 ecm! by the following equations: 


1 X= 374854 924m —161n 1/X= 37425 — 986m — 162n 
1 A\=37398+924m—161n 1/A = 36479 —991m— 162n 
1 X\=37485— 980m —162n 1/X= 30412 —991m—162n 


where #7 and » are assigned integral values. In one band (2667.4) some of the 
fine structure was observed, the wave-length difference being about +0.5A. 
This agrees with the value reported by Henri* in the absorption spectrum. 
It is interesting to note in this connection that the wave-number differences 
991 and 162 cm! appearing in these equations are found also in the Raman 
spectrum of benzene. It is hoped that before long some definite significance 
can be atti ched to these values. 

It may be well to review the evidence supporting the assumption that this 
spectrum originates in the benzene molecule for in spite of all precautions 
some decomposition did occur. At the end of the 12 hour runs the surfaces of 


§ T’enri, Structure des Molecules, Paris, 1925, p. 109. 
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the aluminum electrodes were found to be covered with a laver of carbon and 
the discharge tube showed a very thin film of a shellac-like substance adher- 


2898 


2809-10-12 


Benzene lines 


2737-9 


2689 
2678 
2673 


2667 
2657 


2614 
2609 
2603 


ing to the walls. However, despite the appearance of these products there is 
little doubt that the spectrum observed is actually that of benzene. Apart 
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from the fact that the decomposition products seem non volatile, the follow- 
ing considerations point to the same conclusion. 

First, it was observed that conditions which produced the greatest visible 
decomposition gave spectra of low intensity. Thus, for example, as the 
pressure is lowered decomposition is decreased while the intensity of the 
glow is increased. This is just what is to be expected since, for a given energy 
input, the larger the amount of energy used in molecular rearrangements, the 
smaller will be the fraction of the total energy as light. 

Second, the regularity exhibited by the spectrum seems to militate 
against any decomposition hypothesis, for this regularity implies a uniformity 
in the system which produces it; if this system is a benzenoid molecule, its 
symmetry is easily understood, but if it is produced by some dissociation or 
rearrangement under the stress of the discharge then it is necessary to postu- 
late an absolute uniformity of decomposition which would be hard to concede 
from a theoretical view point and would certainly not account for the fact 
that several decomposition products are actually found. 

Aside from these considerations, which involve only the emission spectra, 
the results can be checked with the spectra in solution. For instance, Dick- 
son‘ has observed the fluorescence spectrum of benzene in alcohol solution. It 
is well known that a solvent shifts the spectrum of a substance in solution but 
this shift is generally a translation of the system as a whole so that the rela- 
tive positions of the bands remain unchanged. It is found that if 19 units of 
wave-number are added to the values given by Dickson the broad bands of 
his spectrum coincide with the principal emission bands which we have ob- 
served. 

Even more convincing is the agreement between the absorption and emis- 
sion spectra of the vapor, which prove to have identical band arrangements. 
The emission spectrum should include the bands found in absorption but the 
far ultra-violet part of the former is lacking due to the absorption by unex- 
cited benzene molecules in the discharge tube. However, it is found that the 
most refrangible emission band group coincides with the bands of the red end 
of the absorption spectrum as determined by Henri.’ This arrangement is 
shown in the following table. 














TABLE II. 
Absorption Bands Emission Bands Absorption Bands Emission Bands 
of benzene vapor _ of benzene vapor of benzene vapor of benzene vapor 
after Henri (A & B) after Henri (A & B) 
r r ny r 
2528 .6 2528.6 2599.9 2599.7 
2541.7 2541.5 2603 .0 2602.6 
2552.5 2551.7 2610.2 2610.1 
2589.0 2588 .7 2613.9 2613.5 
2595.2 2595.2 











4 Dickson, Zeit. Wiss. Photo. 10, 166 (1912). 
’ Henri, J. de Phys. et Rad. 3, 181 (1922). 
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This check is sufficient to establish the fact that the carrying structure is the 
the same, and, since no decomposition is to be expected, nor indeed is any 
found, in the absorption work, it follows that both spectra are due to unde- 
composed benzene. 

The foregoing considerations, in addition to demonstrating that the ben- 
zene molecule is the carrier of the spectrum, also indicate a very satisfactory 
agreement between our work and the known absorption and fluorescence 
spectra. It may be of interest to conclude with a comparison of our results 
with those of McVicker, Marsh and Stewart. Taking the results from both 
spectrographs we have recorded the wave-length of about 110 bands and in 
addition to more than doubling the known number of emission bands have 
given evidence of two new band systems in the far ultra-violet. We do not 
claim to have measured all the bands in the spectrum, however, for in many 
cases bands could be seen on the plates with the naked eye which were not 
susceptible of measurement under the microscope. We have also shown that 
many of the lines reported by McVicker, Marsh, and Stewart are really 
doublets. A number of our measurements when expressed in wave numbers 
seem to be about 30 units smaller than those of the corresponding lines in the 
work of McVicker, Marsh, and Stewart. This is interesting on account of the 
fact that in the comparison made between their results and Hartley’s values 
for the absorption of benzene approximately 20 units were added to the latter 
in order to make the regularities coincide. 

In conclusion the authors wish to express their thanks to Dr. John John- 
ston at whose suggestion the work was undertaken, to Dr. A. E. Ruark who 
gave valuable assistance with the experimental work, to Professor Uhler for 
much help in the spectroscopic and photographic work and for suggesting the 
methods of measurement and calculation, to Professor McKeehan who kindly 
placed the resources of the Sloane Physics Laboratory at our disposal. The 
authors wish also to thank the Commonwealth Fund for the fellowship which 
one of them held and without which he would have been unable to undertake 
this research. 


























MARCH 1, 1930 PHYSICAL REVIEW VOLUME 35 





NOTE ON THE THEORY OF THE INTERACTION 
OF FIELD AND MATTER 


By J. R. OPPENHEIMER 
BERKELEY, CALIFORNIA 


(Received November 12, 1929) 
ABSTRACT 


The paper develops a method for the systematic integration of the relativistic 
wave equations for the coupling of electrons and protons with each other and with 
the electromagnetic field. It is shown that, when the velocity of light is made infinite, 
these equations reduce to the Schroedinger equation in configuration space for the 
many body problem. It is further shown that it is impossible on the present theory 
to eliminate the interaction of a charge with its own field, and that the theory leads 

to false predictions when it is applied to compute the energy levels and the fre- 
quency of the absorption and emission lines of an atom. 


HE relativistic theory of the interaction of electrons and protons with 
each other and with the electromagnetic field has been developed in two 
| papers.! The theory is developed in close analogy to the corresponding class- 
ical theory: the field is on the one hand determined by the configuration of 
the charges; and the motion of the charges is affected by the field. The in- 
teraction between two charges is not then, on this theory, expressed directly 
as a function of the configuration of the charges, but as the effect on each of 
the charges of the field induced by the other. On the classical theory this pro- 
cedure involves grave difficulties, because each charge reacts also with its 
own field. The proper energy of this interaction is, for point charges infinite; 
and it depends upon the motion of the charge. On the classical theory one 
tried to avoid this difficulty by ascribing to the elementary charges a finite 
size; but it was not possible to carry through the theory in a way that was not 
completely arbitrary; nor was it possible to make the work relativistically 
invariant. Oneof the purposes of the present paper is to see in how far these 
difficulties persist in the quantum theory, and in what measure they render 
impossible the application of the theory. 

We may recapitulate briefly the main points of divergence between the 
present quantum theoretic treatment and the classical theory. In the first 
place the state of the matter is here represented, not by a trajectory, but by a 
wave function. Further, the Hamiltonian for the matter is that derived from 
Dirac’s linear wave equation, and not from the quadratic wave equation 
which would follow from the classical relativistic Hamiltonian. Finally, both 
the material waves and the electromagnetic waves are quantized, the matter 





1W. Heisenberg and W. Pauli, Zeits. f. Physik 56, 1 (1929); ibid. in press. The second 
of these two papers is referred to in this work as LC. I am greatly indebted to Professor 
Heisenberg and Professor Pauli, not only for the opportunity of seeing their work before its 
publication, but also for their very valuable criticism and advice. 
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to make the particles satisfy the exclusion principle, the field to make the 
quanta satisfy the Einstein-Bose statistics. This precedure leads to a formal 
difficulty; for the fourth Maxwell equation 


div E—4rp=0 (1) 


is inconsistent with the quantum conditions, according to which there are 
functions of the electromagnetic potentials which do not commute with div 
E but which must commute with the charge density. The two papers of 
Heisenberg and Pauli are distinguished chiefly by different methods of re- 
solving this difficulty. In the former paper new terms were added to the 
fourth Maxwell equation to make the new equation consistent with the quan- 
tum conditions; in obtaining physical results these terms were to be made to 
vanish. In the second paper a much more satisfactory method has been used, 
which takes advantage of the fact that the left hand side of (1) is a constant 
of the motion for all systems involving matter and radiation; and it is shown 
that this constancy follows from the gauge invariance of the Hamiltonian 
for all such systems. The solution of the dynamical problem thus reduces to 
finding a wave function for the coupled system of field and matter, which 
makes the Hamiltonian for the coupled system a diagonal matrix, and which 
in addition makes the left hand side of (1) vanish. The wave function has 
thus to satisfy not only the Hamiltonian wave equation, but also a series of 
wave equations which express the fact that (1) is satisfied at all points in 
space. It is from this set of wave equations that we shall start in this paper; 
we shall write them first in the form given in LC Eq. (68), in which the wave 
function is taken as a function of the Cartesian coordinates 


qr = (qr , qr , qr) -P=1,2,---,N 
and spin variables 
op; P=1,2,---,N 
of the N particles in the system, the number of quanta 


Mn 
of frequency 
Vy 


vector of propagation K,=(K,', K, K,°); K,=v,/c and polarization \, and 
finally a third set of variables P,3, which are essentially the constant com- 
ponents of the electric field parallel to the vectors of propagation K,. We 
shall first show that it is possible to eliminate the P,s's from the wave equa- 
tions in such a way that (1) is automatically satisfied. The reason why this is 
possible is that the condition (1) determines div E precisely—instead of deter- 
mining only the relative probabilities of different values of div E—when the 
configuration of the charges is known, so that the values of the P’s so deter- 
mined can be put at once into the Hamiltonian. When this is done the 
variables Q,; canonically conjugate to the P’s must disappear from the Hamil- 
tonian, since otherwise the Hamiltonian would not be consistent with (1); 








— ny — 
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we shall show that this is in fact the case, and then proceed to an investiga- 
tion of the resulting Hamiltonian. 

It should be observed that the wave functions must ewe in addition to 
the wave equations in configuration space, the condition that they be anti- 
symmetric in all the electrons of the system, and antisymmetric in all the 
protons of the system. Only wave functions satisfying these conditions are to 
be considered in this paper; and it will therefore be unnecessary to indicate 
the antisymmetry of the wave functions explicitly. The fact that electrons 
and protons satisfy the exclusion principle is largely irrelevant to the difficul- 
ties discussed in this work;for these difficulties persist even in the one-electron 
problem. On the other hand it is to be hoped that the resolution of the errors 
in the present theory will make the heuristic postulate of the exclusion prin- 
ciple unnecessary. 

We shall look for a solution of the wave equation in which both wave 
function and characteristic values are expanded in powers of v/c. It has al- 
ready been shown by Breit? that, when the interaction of the particles may be 
treated as small, and radiation processes may be neglected, that these inter- 
action terms give a contribution to the energy of the system, which, in second 
order in v/c, agrees with that computed from Breit’s equation.’ We shall not 
retain the assumption that the interaction terms are small, and shall show, 
somewhat more generally, that when the proper energy terms are systemat- 
ically neglected, the characteristic energy levels and the wave function are 
determined by Breit’s equation. This equation is, of course, not relativisti- 
cally invariant; and it takes no account of radiation processes. In order to 
remedy these defects one must retain the proper energy terms; and we shall 
show that it is then possible to make a formally satisfactory theory to give the 
shape and position of all spectral lines, and the energy of the normal state. 
The theory, is, however, wrong, since it gives a displacement of the spectral 
lines from the frequency predicted on the basis of the nonrelativistic theory 
which is in general infinite. This displacement arises from the infinite inter- 
action of the electron with itself; this interaction depends upon the state of 
the material system; and the difference in the energy for two different states 
is not in general finite. Thus the present theory gives no more than the non- 
relativisitic theory of Jordan, Klein and Wigner.‘ It seems improbable that 
Breit’s ees gives the energy levels of an atom correct to second order 
terms in v/c; but we shall see that there is ground for supposing that it does 
give the separation of the fine structure levels correctly in this order. On the 
other hand the displacement in the frequency of the spectral lines which arises 
from the proper energy should be of the second order in v/c, and is thus larger 
than the natural line breadth, which is of the third order; and on the present 
theory it is not possible to compute this displacement. We shall return later 
to a consideration of these difficulties. 


2G. Breit, Phys. Rev. 34, 553 (1929). 

3 Breit, ibid., Eq. (6). 

‘ P. Jordan and O. Klein, Zeits. f. Physik 45, 751 (1927). P. Jordan and E. Wigner, ibid. 
47, 631 (1928). 
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1. The condition that the left hand side of (1), regarded as an operator on 
the wave function y of the variables 


adp,9P, Mn, Ps 


shall make the wave function vanish at every point in space, gives a series 
of wave equations 


I Cv= | Pat Leve(ar) ar ,op,M,y, P,3) =0 
P 


which must hold for all values of the K, consistent with the boundary 
conditions. Here the v’ are functions defined in LC (54): 


2 
v°(qp) =—(ck,3L*)—!? sin rk, gp™ sin rk, gp sin rkpqp® (2) 
T 


Here Lis the length of the fundamental cube, or Hohlraum, and is taken finite 
in LC toavoid the introduction of a continuous manifold of normal coordi- 
nates; in all our results we shall make Z become infinite. Furthermore the 
Hamiltonian for the coupled system, regarded as an operator on the same 
wave function, gives the wave equation 


( N 
II \- E+ >| T Malin + arate +2 [se gradp) + mpay?c? ‘|} y 


r A=1,2 


N 


+4 Dd lerAo%(gr) + er(a?- Ava) by 


P=1 


+i a * re (Ma+1)"2Aq—!— Ma! “Sal fy 


r A=1,2 P=l 
+} 2 YurQaby=(+H- E)\y=0. 
r P=} 


Here A 0° is the external scalar potent al, and A® the external vector potential; 
and A, is an operator which transforms 1/,, into M,—1, and leaves all other 
variables unchanged; the w”’s are operators which operate only ong,, and are 
derived from the Dirac matrices a,,' by the definition 


a,"F (op) = Lice p'ppF (pr) ; aol F (op) = Diaep spl (pr). 
ep 


pp 


Further the Q,; are canonically conjugate to the P,3 so that 
h 
[ PrsQr’s | =—F5,, (3) 
21 
Finally the functions u,”, up,” are defined in LC (59): 


m2 S 
up’ = ere(= : ) Deas?2:"*(qp) 


l=] 
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h \'2_3 
up? = cre ) DoaPo(gr) ; A=1,2 (4) 
2rv, l=1 
S 1/2 » 
v1" (gp) = (=) Fy cos rhigp sin rkyge™? sin rky-ge™” ; 
1=1,2,3. 
A=1,2,3. 


The square matrix Fy, is given by the scheme 


A= 1 2 3 
| : 
1 | €2"/(€17 + €2")'/? — ey€3/(€:°>+e2")'/? ey 
2 —e€, /(ex> + €2”)"/? — €n€3/(€:* +2")? ee (5) 
3} 0 —(€,?+6€”)'/? és 


where the €,"’s are the direction cosines of the vector K,. It should be ob- 
served that 


ic. h 
——(a” gradp)epv'’(qr) =——up™. (6) 
2ri 2rt 
This equation, together with (3), shows that 
[HC,]=0 (7) 


for all r, so that all the C,’s are constants of the motion. The equations I and 
II are those given® in LC (68). 

From I we see that the wave function must be singular in the P’s. We may 
avoid the use of singular functions by making a contact transformation from 
the variables P to Q, and writing the wave function as: 


V(gr,or,M,Q;3)- 


For we may then solve (3) by taking 








h 0 
Qi 80x 
If now we set 
v(qP OP, M.,0rs) = em 2ris/ht Zep! (ap) Qragy (qp,oPp Mn Ors) (8) 
the equations I give us 
dg 
~ (9) 
00,3 
for all r sothat 
te) =$(gp,oP, Mn) ° 


5 In LC (68) the A’s are dropped, 
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Further the Hamiltonian II becomes 


en 2tilh Devo gr Qe} -F+> | Xu M yhv,+ 1, Can o”"(qp)0"°(qp’) 


2 P,P’ 


0 h? o° > 
+ ih 0" gp)—— —- — ——-+ | a? gradp) + mpc? 
2 a ” i 00,3 tr a 40,3 x 2 ( 7 ; ’) — - 


—¢ r Ee (a? gradp)o"(gr)Qrs+ Ler Arar) +(e" A°(qr)) | (10) 
+ > Dae > 2 atti 
r r P 





— Mp! An] t o(gropM,.) =0. 


The terms in Q drop out because of (6); the terms in 0/0Q give nothing be- 
cause of (9). The equations I and II thus reduce to the single system 


\~EtHot & Mahe, + i Dow (Mat) dat = Mal “an] fo=0 





rAP 
sl (11) 
Ta 
Ho= r45-j@" gradp) + mpc?ao” + pp|Ao(qr) + (a? A(qgr)) 
P 271 
+r z, Dever-v,0"(qr)o™( gr’) ; ‘ 
at r 
It is this system which we must now investigate. 
The terms 
G(qr,qr’) = > rv,v"(qr)o"( qr) 
may readily be evaluated,® and give for L— x 
1 
G(qr gr) ——— 3 ree = | Qe—Qp’| . (12) 
‘ 2rpp: 


The terms for p and p’ different give the electrostatic interaction of the two 
particles; the terms for p=)’ give the infinite electrostatic proper energy of 
the particles; on the present theory it is not possible, as it was on the non- 
relativistic theory,‘ to eliminate these terms; the physical ground for this 
impossibility has already been indicated, and lies in the fact that the field 
acting on any particle is the sum of the fields induced by all particles; it is a 
consequence of the principle of superposition for the field. These electro- 
static proper energy terms do not, however, interfere with the application of 
the theory, since they are constants, and may be dropped from (11) without 
altering the form of the wave function. We shall find other infinite proper 


6 W. Heisenberg and W. lauli, Zeits. f. Physik 56, 1 (1929). Eq. (115); G. Breit, reference 
2 Eq. (57). Breit has independently evaluated the P,; terms in the Hamiltonian; and I am 
much intebted to him for informing me of his result 
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energy terms in the course of the work; but these will turn out not to be con- 
stants, but to depend upon the configuration of the system; dropping them 
does alter the form of the wave function. 

If we now neglect the coupling between matter and the light quantum 


1 field the wave equation reduces to 
| - B+ Hot E Maly |o=0 (13) 
rh 
and for the case that no quanta are present we have 
—E+ 2X | ta" -gradp) + mpc?ag” + ep|Ao"(gr) + (@”- A(qr) | 


(14) 





epep’ 
- } — \ o(or,gr) =0. 


PP’ 2rpp: 


We shall show that the terms which we have neglected in (12) are small of the 
order (v/c)?; and by neglecting other terms of the same order, (13) can be 
considerably simplified. For consider first the equation 











{ ' os , 
i E+ >} | —(@?-gradp)+mpcta? | >o=0. (15) 
; p L2mt 
For N-free uncoupled particles. If we choose all matrices|| @,,'.,,|| of the form 
0 0 6 a 
Oc d 
are (16) 
b « 0 0 
ad0 0 
and all the ||a,,°%.,|| of the form: 
\ +1 (0) 0 0 
0 +1 0 0 
(16) 
0 0 —-1 0 
0 0 0 —-1 
and satisfying of course 
[a*-Pa”?’|=0 for PAP’; [a Pa’? }+ =i (17) 


then any ¢(¢,) in which of the a,’s have either of the values 3 or 4 will be 
small compared with any of the @’s for which all of the o,’s have the values 
1 or 2 of the order (v/c)". Now the terms 


, e€pep’ 





? epAo'(qp) 
Ppp 2rpp 


in (14) do not involve the a?’s, while the terms 


ep(a?- A(qr)) 
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are small of the order v/c; thus as v/c is made to vanish, all the solutions of (14) 
vanish except those for which all the o,’s have the values 1 or 2; and (14) 
reduces to 


h? €pep: 


r= TON oongr) =0 (18) 





| - E+ >} mpc? + epAo(gr) — 
P 


for all o,’s=1 or 2 and 


o(opgp)=0 for any op=3 or 4 


which is the Schroedinger equation for the N-body problem.’ It would thus be 
possible to take (18) as the starting point for our systematic solution of (11) 
in powers of v/c. We shall not do this, however, as it would complicate the 
analysis, and lead to no new results. We shall thus give up the assumption 
that the a’s are written in the form (16), and take for our zero’th approxima- 
tion to (11), the solutions of (14).* 

2. We may write the equations (11) seriatim: 


EM a= M=0;(—E+H0)0=—i Dur (1a); 6=9(0n) (11.1) 
ry rAP 
a (OO Et het Hoon) = +3 Cunro—| iD Teta ten) | 
— . ‘ ail (11.2) 
( — 1/2; > uo(2,) 
P 
( (-— E+ h(vp- +r) +Ho)o( 1a ’ 1,°) = +i } i [up™(1--r*) 
P 
M=2 | tur" o(1n)]—i Do” Dou’ bla, dear, bern) (11.3) 
| rr PP 


| = 227 Do [up (2, Lene) Fur’ (1a, 2ra°) | 
\ P 


| (— E+2hv, + Ho)o(2n) = +2"2i Doup™( 1a) 
hs 


) 
—i dD Yue’ o(2n, tea) —3"/24 Do upO(3r) 
\ r’d’ P (11 4) 


etc. 
Now for fixed 7, and fixed K,, 
up =0(c-¥2) 
whereas 
vp =O(c). 
7 It is possible to write the two component wave equation when the magnetic interactions 
are retained up to the order (v/c)?, as has been shown by Breit, reference 2. 


* We shall not make explicit use of the fact that the a’s are in the form (16); we shall, how- 
ever, retain the assumption that the a’s are small of the order v/c, to obtain u»™ = O(c”), 
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Thus we should expect ¢(1,,) to be small of the order c~*/?, #(1,,,1,’x’) to be 
of the order c~*, and so on, and we should try to find a solution of (11) of the 
form 


E=E©+EM4E®... 


(19) 
o(M pn) =o (Mn) +O(M ap) «> - 
with 
E™ =0(c-*) 
@")( Myx) =O(e-*- 3/24) , 


It should be observed that in general there will be certain frequencies for 
which $(1,), (1, 1,’,’) ete. will not converge uniformly for c—>© ; and that 
for these frequencies the expansions (19) will be illegitimate. The frequencies 
for which this convergence is non-uniform are those for which 


+E-— }Malv, 
rh 


is a characteristic value of the homogeneous equations 


Such frequencies will not occur if (20) has no solutions for \<&, i.e. if the 
material system is in a normal state; but in general the expansions (19) must 
be modified; we shall return to this modification later, and shall see that it 
gives a satisfactory theory of the absorption and emission of radiation; but for 
the present we shall assume that the atom is in a normal state, so that (19) 
is justified. 

On the present theory there is no normal state for the matter, because 
states of infinite negative energy are possible; one may in fact show that, on 
the present theory, Dirac jumps to such states from states of positive energy, 
jumps in which the energy and momentum lost by the matter are taken up by 
the field, are not only possible, but infinitely probable. But that the theory 
should predict this is a token of an error in the theory; and since the Dirac 
jumps do not seem to be directly responsible for the difficulties with which 
we are, in this work, most concerned, we shall for the present neglect them. 

We shall first give a complete solution for the case that we drop all proper 
energy terms, for the case, that is, that in all double sums of the form 


>F (gr »OP ,P’ op’) 

ad 
we may set the terms with p equal to p’ equal to zero. This solution is not 
unique beyond terms of the second order in v/c; for in the higher orders it is 
no longer possible uniquely to separate proper energy and interaction energy. 
But we may readily obtain a possible solution: 


wif > up™ My 
P 


¢(Mn)= J] }] ————-| 9 (21) 


rh h Vy 
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where 


| - na~t oe = “| =0. (22) 


rh PP’ h Vy 


For if we put these values, for example, in (11.2) we get 


+i > up™ 
at Noanes > ee +o 


hv, rh’ PP hv, 
(23) 
jan Can |o =1 Dour. 
hv, P 
Now in 
he : a fitian 
>| H ,ep”™ = ; i —— -gradp) + mpc?ag? + ep(a?: A(qp) ,ap™ (24) 
P PP’ 271 
and in 
aa 
PP’ 
we may put 
> ED ay 
PP' Pp’ 
so that 
> [Hue] =0 (25) 
P 


and (11.2) is satisfied. In a similar way it may be shown that all the equations 
(11) are satisfied. 
We may evalute the terms 
up” up” c Pp? v:”*(qp)0v"™(gp’) 
= = €pep»— doa ay’ — 
} 2 2 
rh iV, iT il’ rr Vr 











in (22) by observing that for /¥1’ 
vi" (qp)0v™( gp: ) 0 0 




















us 
> = ——F(gpqp’) ; ApPF =——— 
” vy ~ aq'p ' Op. Crpp: 
and 
v ™(qp)vr (gp) —T eo 
} ie - = —— ep gr’) 
rh v;" 2c? “T pp’ dq'p 
so that 
wp™up" — erer” | ee) eee or) \ an 
te (27 
rh hv, 4 \ Tpp’ r pp: { 


This gives for ¢o 





1 , pep: 
‘ —E+ mt. , —- [CaP a )rnp +(e epe)(ae”-top.)|bo=0 (28) 


pp T"pp’ 


This is the equation used by Breit.’ It is patently not relativistically invariant ; 
this means that the proper energy terms are not invariant, and forces us to 
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retain these terms, at least in part. Furthermore, we have not, in the deduc- 
tion of (28), used the fact that the atom is in its normal state; in spite of this 
there is no sign, in the solution, of processes involving the emission or ab- 
sorption of radiation; for these processes arise from the interaction of the 
particles with their own field. We have, therefore, to consider the solution of 
(11) when the proper energy is not neglected; the retention of these terms will 
preserve the invariance of the theory,and give usan account of radiation proc- 
esses, but it leads to results in contradiction with experiment; it makes the 
validity of (28), even to the second order in v/c, doubtful. 

3. Wecan readily find a solution of the form (19) when E°and¢ correspond 
to a normal state for the matter; but we cannot find this solution in closed 
form; nor is there an equation in configuration space, corresponding to (28), 
for ¢. If we put (19) in (11.2) etc. we get 


( rh 
rh | Don = > ° f an eee dqx‘*? 
A bomdm ”? | P 
o(1,) =i Zz aie 


m hv, + En — Ey . | ps Sno | (29) 


— 





\ (Ho— En Som =0 
and further 
rh _#’h? (0) 


bom ban Pn 
(1, len) = — , 
’ \ . . 2X [hv + Em — Eo | [hyp +0) +E, © — Eo | 


























r’d’ rh (0) 
> bom ban Pn ( 30) 
mn [hve + Em — Eo |[h(v-+y,) + En — Eo | 
et = 
Om YUmn 
o 2; j=-—- 21/2 ; etc. 
-_ 2X [hy + Em — Eo | [2hv, +E, — Eo | 
Further E(o') and ¢(o!) vanish, and 
_ Pe 
On 
E, (2) = — 
: “ ~ hv,+ Em) — Eo 
_ a rey 
Om Umn Yan 
(ws ; 
” x X X [hy + Em — Eo |[E, — Eo | 
Moreover @#'(,1,), @(0*) and E(o*) vanish, and 
r’\’er'dh’ rh (0 
bom Gas bak Pk 
$° (In) =i — 
(ln »y u y [hve + Em — Eo | [hy + Ex — Ex |[E, — Eo | 
(32) 


Tr’ er’d’ rh ofr 


Ey? = — > i > > bom Omn Oni Deo 
rh “'y 


mk nO [hv + Em —E, ] [ hv, +E, = E, | [E, —E,| 





etc. 
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The terms 
— Epo (1,,) 
and 


—§ Dar’ G (aden) — 2447 Dur (2a) 


Pr P 
in the equation for ¢“)(1,,) cancel, so that 


o® (1) =0 
and 
oo =0, Ey) =0. 


The expansion for E‘® and ¢‘ can be continued, and only terms of even 
order in v/c appear. 
It will be observed that the interaction terms in (31) 


- ] 
= > > [hv + =. -— Ey | | f 200 urroorda | Jon umroonaal (33) 
ry m | il 


| we write fea for f tes faa - ++ dqy® p | 
o’:+-0N 


differ from those computed for the second order from (28): 
P | ma | 
~X Lon] fdrmrovindal fsomrrerda| — 8 
rd m i 


It is possible to express 


up” up:™ 


rd hv, + En — Ey 





in terms of the confluent hypergeometric functions, but the expressions are 
too complicated to be suitable for calculation. 
The proper energy terms 


9 





-xz >! a- $14, 


m hv, +E, — Ey rAPm 





do not exist, although both 


rh 
p 2 7 Om ,P 
rh 
and 
rh 
be Tom r 
m 
converge; for large p,, 


> Tim p= 0(1) 
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and for large E 


rh 
> Tom,p =0(1). 
rr 


The energy level of the normal state is thus infinitely displaced by the 
interaction of the particles with the field; the question which we have now to 
consider is whether or not the energy differences between two states are dis- 
placed by a finite or an infinite amount. 

In order to answer this question we must treat the case of energies which 
do not correspond to a normal state for the matter; and we have to modify 
(19), and take account of the emission of radiation by the system. And for 
this purpose it is convenient to make the dimensions L of the Hohlraum in- 
finite, because that makes the physical interpretation of the solution more 
immediate. For then vy, becomes continuous and we may normalize the 
v:"(q») to the intervals AvAw, where Aw is the element of solid angle of the 
unit vector €,. Furthermore we may treat here, to simplify the writing, the 
case that E corresponds to the first excited state of the atom, so that there is 
only one energy £» lower than E for which (20) is soluble. 

We define the energy of the normal state by 


Eg=Ep™+ Ep? + Ey - - - (35) 


where E, and Eo“ are given by (31) and (32), and we define the corres- 
ponding wave function 


Up = 0? +40? +9” ae (36) 


where ¢o” is given by (31). We can then extend this definition formally to 
obtain the energy and wave function of excited states: 


En=EnO+En@+En™ «as 
TA 19 
| binn | 
- i», f dw, 37) 
f adv Ww u X hv, +E, © — En ( 


bm +om™ +om™ -+- etc. 


En® 





Um 
But in the expressions for £,,°) etc, and ¢,,°’ etc, the integrals over v are now 
improper, and we have to displace the path of integration around the singu- 
larities. This is equivalent to replacing 
1 


bw +E, —E, 





for E,<E, 


by 





+4 C, () — Ey, 6 38 
oes ak ix6(hy+ En) (38) 
and then taking the principal value of the integrals over v. Then in general all 


the E,,’s except Ey are complex. We now transform the ¢(M,, ¢,0,) by the 
formulae 
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o(Ma.m)= >> f sae f age > ++ dgy®un(qr,op) (39) 
eo} 


o(Mx,gP,op) 


mt p= rf. Sf age - ++ dqgy@tippp™n, . (40) 
oP 


Then if the @(.1/, m) satisfy the equations, which follow from (11): 


-E 
(— E+ Fn (0,m) =~ i fel, dw, es Cun of ors (v4) 
n 
— fu du, >, Spe fa! dw,’ >> F pnon pF a( co, 00," , 07’ ,d,2’) (41) 
M P é' 
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and introduce 





’ P 
Eo-E 
i(», +p, +=) a oi? 
h 
, 4 ry ‘ 1 176 mo I m “)! 
(1,1) =t 2D bam. pO(O,2)- _— +750, += 
' Letom t \-E+u+E.. & ait 


Pr 


+ Fy(w,r)6 (> + os) —4% p f dy,’ fe / pr P (42) 
d’ 


P 





Eo—E 
F(a," , wr Yr, XN N’)6 (> +y,/+ ) er 


7 


F r’d’ rr 
(Lin, Lene ym) =4D, | Bam PO(Lar, 2) + pn POL, 2) } 


Pm 


i) 1 ; 176 mo s( + +")! 
oO erry I? 
\-E+h(,.+,)+En fh h JS 


-—E 
+ Pale! 0! W,X)0( 1 ety,’ +") oes 





1 
where now the functions F are completely arbitrary, then the ¢(Mn, gp, o>») 
satisfy (11); for each choice of these F’s we can obtain a solution of (11). Now 
the terms 

rh 
Mn m POCO, n) 
—E+b,+En 
give a radiation field which on not extend to infinity; and the terms 
+ Sipe, p Ey-—E 
ae ,2)5 mod ; etc. 


U U 





represent outgoing® electromagnetic waves, so that by the choice of the F's 


* P. A. M. Dirac, Zeits. f. Physik 44, 585 (1927). 
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we may determine the radiation incident upon the system. The simplest case 
is that in which only quanta of the single frequency 


1 
ym B— Be) =), (43) 


are incident upon the system, so that all the F’s vanish except 
F,(w,d) . 


For this case 


i > of diction pF (eo ,X) ° 
» ?P IG m” 


0, =+ = . 
o(0,m) E-E. E-E. 





(44) 


Now by hypothesis E£ is to be chosen that only one E,,—E,that for m =1, say, 
is to be small, so that only 


iG,” 


E-E, 





¢(0,1)= (45) 


is large. The probability of absorption to this state is thus proportional to 














| G,’| 2 | G,’| 2 
- = (46) 
|E-E,|? | w+E.—E,|? 
so that the shape of the absorption line is given by 
const. 
(47) 


2 


1 
v+—(Eo— Ej) | 
h 





since G’,, varies slowly with v. 

If we evaluate Eo, E, to the second order in v/c,and drop the higher terms— 
and this is equivalent to neglecting transitions in which more than one quan- 
tum plays a part—(47) reduces to 


const. 





. (48) 


1 1 1 i 1 ie 
natal Pee ae A | (2) quetens Ce) am 2 
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The absorption line is thus of the same shape as that predicted on the basis 
of the correspondence principle, and that found, for this case, by Dirac®, and 
the half-breadth of the line is 


E,®—E,® 
| 2h 


1 








(49) 
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where 7; is the natural life time of the state 1. 
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The center of the absorption line is displaced to the red from (1//) (4, — 
Eo‘*’) by an amount 





rh | 


| b,, | | | Bon | 2 \ 
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Here the principle values are to be taken for all improper integrals over v. 
The terms for p# p’ are just those to be expected from (33) for the displace- 
ment of the energy levels by the magnetic interaction of the particles. The 
terms 

rr 


are YEO Ww+E,O— Ey f 





must be ascribed to the effect of the interaction of the particle with its own 
field. They may be compared with the formula obtained’ for the same effect 
by Dirac, who finds a displacement 


phe fu pp Lak (52) 
hy+ Eo — 7X0 - 


There does not appear to be any justification for this result, because in its 
derivation terms were neglected that are of the same order as those retained. 
But it is of interest to observe that the integral in (52) exists, and gives a 
finite displacement of the line of the second order in v/c. This displacement 
is thus larger than the natural line breadth, which is of the third order. It can 
be computed® when the w’s are known. Thus for the first Lyman doublet of 
hydrogen we find the same displacement, in this order, for both components; 
it turns out to be 





ch 1 ch 


ee ea eee ae 
32r7e" 71 327°e 








which is about forty times the line breadth. The fact that this term, and the 
similar terms in (51), are the same, in second order, for the two components 
of the doublet, suggests that the formulae (33) in which the proper energy is 
neglected, will give the atomic fine structure splitting correct to the second 
order. 

If we try to compute the displacement from (51), we find that the in- 
tegrals over v diverge logarithmically for high frequencies. One can readily 


9 The calculations of the displacements predicted by the results of Dirac were carried 
through in collaboration with Harvey Hall; and I am indebted to him for permission to quote 
them here. One must use the retarded potentials to obtain a convergent integral. 








INTERACTION OF FIELD AND MATTER 477 


see that this is not the result of the neglect of higher order terms, nor of any 
of the approximations made in the work. The theory thus leads to the false 
prediction that spectral lines will be infinitely displaced from the values com- 
puted by the Bohr frequency condition. The behavior of the expression (51) 
calls for some comment. As the formula stands, the integral over v diverges 
absolutely; this may be verified by evaluating the terms for a free particle. 
But the question arises whether it is possible so to rearrange the order of 
the integration over vy the summation over n, and the two integrations involved 
in the evaluation of the um,’s by (40), that the limit »—*, v->= exists. 
This cannot be effected by an interchange of the sum over m and the integral 
over v; but there is a procedure which, when the £’s in the resonance de- 
nominators of (51) are dropped, does give an absolutely convergent result. 
This procedure was suggested by Heisenberg, who showed that, if we first 
perform the integration over vy and w and the summation over A, then sum 
over all the states » of the same energy, then sum these up to some large but 
finite energy E, take the difference of the two terms for the state (0) and the 
state (1), then perform the two integrations over the configuration space, 
and finally allow E to become infinite, the limit E- = exists, and (51) tends 
to zero. But if we try to apply this procedure to (51) when the E terms are 
not dropped, we get for the leading term the divergent result 


de? /he- (Ey —Ey™)/h - [dx ‘x (53) 


Nor is there any method for obtaining an absolutely convergent expression 
for (51). It should be observed that (53) gives us another justification for 
using (33) to get the fine structure separations correct to the second order. 

One can see quite simply that (47) ought not to give a finite line dis- 
placement. For consider two states of a free particle; in one let the particle 
be at rest; in the other let it have the velocity v. Then if the energy of the 
particle at rest be E, the energy of the moving particle in proper coordinates 
moving with the particle will also be E. But we know how this energy 
transforms under a Lorentz transformation; in the original coordinates it 
will be 


EB B=[1—(e/c)*]-"? 
and in the same coordinates the difference in energy of the two states, which 
gives the line displacement, will be 
E(g—1) 
But this can only be finite if £ is finite which, by (31), it is not. 

We have treated these difficulties in some detail, because they show that 
the present theory will not be applicable to any problem where relativistic 
effects are important, where, that is, we cannot be guided throughout by the 
limiting case c—*. The theory can thus not be applied to a discussion of the 
structure of the nuclei. It appears improbable that the difficulties discussed in 
this work will be soluble without an adequate theory of the masses of electron 
and proton; nor is it certain that such a theory will be possible on the basis of 
the special theory of relativity. 
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ABSTRACT 


The leads of a tungsten filament in vacuum cool the ends of the filament and so 
affect the voltage, candle power, electron emission and other properties of the fila- 
ment. For long filaments, where there is a central portion at a uniform temperature 
T,,, the temperature distribution near the lead is derived. A method for determining 
Ty, the temperature of the lead-filament junction, is given. Tables and formulas 
are presented which allow ready calculation of the effect of the leads on the properties 
of any long tungsten filament for which the current and diameter are known. From the 
more general results it has been found that the decrease in voltage due to the cooling 
of one lead may be represented by AV =0.154 (T,,,/1000) —0.081 (7, /1000) —2.1- 10-8 
T. Tm —90.056. There is an extension of the theory to cover the cases of filaments in 
gases, filaments of other materials, etc. 

Part II of the paper gives figures from which may be found the properties of 
filaments so short that the first theory does not apply. Some experimental checks of 
the theory are given. 

In general the results and the methods of application have been placed first, 
and the mathematical derivations have been placed at the end of each part. 

For a short filament with leads cooled in liquid air a negative slope of the volt- 
ampere characteristic when the central temperature is much smaller than T,, is 
observed. 


Part I. THe LONG FILAMENT 


HE extensive use of tungsten filaments in research and industry makes 

it important to consider how the cooling effects of the leads influence 
the characteristics of such filaments. For wide ranges of temperature the 
characteristics of hypothetical filaments which are not cooled by leads may 
be found from tables of the properties of tungsten'-? The magnitudes of 
the lead losses have been evaluated experimentally': *:* and by theoretical 
methods.!” It is felt that there is still a place for a systematic treatment 
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2 C. Zwikker, Royal Acad. Amsterdam 34, No. 5 (1925). 
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8 A. G. Worthing, Journ. Frank. Inst. 194, 597 (1922). 

® T. H. Amrine, Trans. Ill. Eng. Soc. 8, 385 (1913). 
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7, 5 (1927); V. Bush and K. E. Gould, Phys. Rev. 29, 337 (1927). 
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which may be conveniently applied to a filament operated under any ordinary 
conditions." 

Temperature distribution. For most tungsten filaments the cooling effect 
of the leads does not extend appreciably to the central portion of the filament. 
The absolute temperature, 7,,, of this portion may be calculated from the 
diameter of the filament and the current through it.!* The temperature of 
other parts of the filament is best expressed as a fraction, 0, of 7. Thus 
6=T/T», where T represents the absolute temperature of any point of the 
filament. 

We may consider the effect of each lead independently, since the two 
effects do not overlap. Of fundamental importance is the variation of 6 
with x, the distance from the lead. This is shown later to be governed by 
the equation, 

@ = ad6/dx (1) 


where ¢ is a function of 0, and is given by Eq. (30). a is a parameter of the 
dimension of length, and depends on 7, and the’filament diameter D. 


TABLE I. Values of ao for various values of Tm. (For ‘a’ use Eq. (2).) 














Sel x) | ao(cm) | | do » | ao 
600 =| —Ss«S 84 1700 0.646 2800 | 0.275 
700 ~=Sti«‘|Cts« LOS; 1800 "582 2000 261 
300 | 3.01 1900 "527 3000 "247 
900 | 2.33 2000 "481 3100 235 

1000 =| «1.863 2100 441 3200 223 
1100 1.524 2200 "406 3300 | 213 
1200 1.274 2300 "377 3400 209 
1300 «=| s«1.084 S| S200 "351 3500 | 195 
1400 936 ~=—|~Ss«2500 329 3600 187 
1500s 921 2600 309 «| 3655 183 
1600 | 724 2700 21 | 

| | 














do, the value of a for D=0.01 cm (4 mil approx.) is given in Table I. a for 
other values of D may be found from 


a = ao(D/0.01)'?. (2) 


Integration of (1) gives 
8 
(x/a)8 = f d0/@. (3) 
0 


The values of (x/a)é for various values of 6 are tabulated in the second col- 
umn of Table II. Note that a is effectively a unit of length. The symbol 
(x/a)§? will in general represent the distance, expressed in a-units, from a 
point at temperature 6,7,,, to a point at temperature @.7,,. To obtain dis- 


J], Langmuir, Trans. Faraday Soc. 17, 634 (1922). See also references (7) p. 210, (6) 
p. 312, (3) p. 356. In these papers formulas for AV and AVy are given. They were derived by 
methods similar to the ones used in this paper. 

® Ref. (3), p. 312, Table I, column 4. 
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tance in ad-units, divide the distance x in cm by the value of a, in em, found 
from Table I. 


TABLE II. Values of (x a)? and 8(60). 
For n>4: 3(@9) =(x a)%o, For n =1.2, see columns 3 and 6. 





* ) 
0 (v/a )° Fy ), 0 (x/a)? faye 
0.0 0.000 0.000 0.7 0.7394 0.464 
a 0419 040 8 9766 532 
s 1110 102 85 1.1354 
25 1522 137 9 1.3592 598 
3 1974 172 95 1.7260 630 
4 2999 245 99 2.5535 054 
5 4200 320 | 999 3.7224 659 
6 5028 ~~ | 302 | 1.000 | | 660 





The temperature distribution near a cooling lead is represented by the 
curve farthest to the right (labeled 0.995) in Fig. 2 (Part II). The incre- 
ment in the abscissa (x/a)% from the ordinate 62 to ordinate @; gives the dis- 
tance along the filament from a point at temperature 6,7", to a point at tem- 
perature 0.7 ,,. 

In practice the junction of lead and filament will be at a temperature 
1) =6 T,,. The exact determination of this temperature is rather difficult. 
If the leads are short and fairly heavy we may assume 7) = Tr, where Tp is 
the room temperature. The error due to this assumption in the value of any 
filament-property for the whole filament!’ computed by the methods of this 
paper will be less than 1 percent when the length of the lead is less than a 
certain maximum length /). We find that /» is given approximately by 


ly = 0. 32(x/a)(Dz/0.1)? (Az/0.586)/A. (4a) 


A is the filament current in amperes and (x/a) represents the half length of 
the filament. D, is the lead diameter in cm. \, is the thermal conductivity 
of the lead in watts cm~! deg~'. For nickel leads \,/0.586 = 1, for tungsten 
leads \,,/0.586 = 2.73, for molybdenum \;/0.586 = 2.49. 

Thus with nickel leads for which D; =0.1 cm and /<1.6 cm, used with a 
20 cm filament of D =0.02 cm for which the highest operating temperature 
is 7',,=2400° we may assume T)=T7'p, since we find A =4.02," x/a=20.2, 
and thence from Eq. (4a) J) =1.6 cm. 

For leads such as those used in incandescent lamps it is sufficiently 
accurate for many purposes to assume that 7')=(1/4)7T». 

If desired we may evaluate AT =7)—T >, in terms of the lead length 7. 
We find that 


AT =1- A (0.586/Az)(0.1/Dz)? (ATo). (4b) 
(A7), is the value of AT for a nickel lead for which D; =0.1 cm, /=1 cm 


and A=1 amp. It is given in Table III as a function of the value of T,, for 
the filament. Note that these values are for a filament of constant current, 


8 Such as the voltage or the candle power of the whole filament. 





= 


i- 
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and hence that the diameter of a filament for which A7’=(A7T)p¢ is smaller 
for the higher temperatures in the table. 

The data from Table III are not dependable above 7)=1000°, where 
radiation loss is appreciable, nor when the resistance loss in the lead is large. 


TABLE III. (AT)o=T0—Tpr for nickel leads when A=1, DL=0.1, /=1. 





=> — ——— 














Ta(K) | (AT) | T, | (aT || Tm | (Deo 
1000° 21° | 1800° 45° | 2800° 80° 
1200 20 2000 | 51 | 3000 87 
1400 32 | 2200 58 | 3200 95 
1600 38 2400 65 3400 102 


2600 72 





The actual temperature distribution of the filament is given by (cf. 


Eq. (3)) 
8 6, 
(x/a)6, -{ dé/ -f d0/o. (5) 
0 0 


To find the distance from the lead, x, of a point of known temperature @, 
we may evaluate the two integrals above by means of Table II, and then ob- 
tain x from (x/a)$, by using Table I. Conversely we may find the temperature 
6 of a point at any given distance, x. 

By the use of data on the characteristics of tungsten filaments as functions 
of temperature,': * and from the temperature distribution along a filament as 
found above, the properties of the filament at each point can be calculated. 
For example, we can determine the electron emission, the radiated energy, 
the luminous intensity, et cetera, at each point. 

The effect of lead losses on filament characteristics. Many filament prop- 
erties which are functions of the temperature would be strictly proportional 
to the length of the filament if the temperature were everywhere uniform. 
Let 4 be a quantity which measures some one of these properties per unit 
length at any given absolute temperature 7. For example, 4 may represent 
the voltage drop per cm or the electron emission per cm of length. Let 
h,, be the value of h at the temperature 7,,.. Nearly all the properties of 
tungsten which we shall need to consider vary in proportion to some definite 
power of the temperature over rather wide ranges. Thus we may put 


h=h,,6" (6) 
n=(dh/h)(T/dT) (7) 


where ” is approximately constant." 
If a filament of length 2x were all at its maximum temperature 7°, 
the value //,, of any property for the whole filament would be 


H. = 2xha. (8) 


44 For values see Ref. 3, p. 354 Table II or Ref. 1, p. 153 Table I-B. 
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The cooling effect of the leads makes H, the actual value of the property 
for the whole filament, less than J/,,. AH, the amount of this decrease due 
to one lead, may be thought of as the total value of the property over a 
short length of uncooled filament. Designate by AV the voltage drop across 
this length. AVy is then the volt-equivalent of A//J, and the fractional 
decrease of //,, can thus be expressed as a fraction of the total voltage V,, 


AVu/Vm = AH/H ». (9) 


The ratio ///J//,, is a measure of the extent to which the cooling effect 
changes the property. We have 


H/Hm = (1m — 2AH)/Hm = (Vm — 2AVn)/ Vm. (10) 
The factor 2 accounts for two leads. 
TABLe IV. Values of By=f\(1—6")d0/¢. 
| 























| 
n B, n B, ] n B, 
E Ps Cheese : 
1 0.583 9.0 1.626 | = 20, | 2.032 
1.2 | 0.660 10 1.682 | 22 2.079 
2.0 | 0.882 11 1.728 | 24 | 2.124 
3.0 | 1.076 12 1.772 | 25 | 2.145 
4.0 | 1.217 | 13 1.813 | 26 2.165 
5.0 | 1.329 | 14 1.850 | 28 2.203 
5.1 | 1.339 | 15 1.885 | 30 2.238 
6.0 | 1.421 | 16 1.918 35 2.315 
7.0 | 1.500 | 17 1.949 | 40 2.384 
8.0 | 1.566 | 18 1.978 | 50 2.497 
| 19 2.006 60 2.589 
| | | ; 

















If V is the actual voltage drop (in volts) and AV the value which A// 
has when the property measured by JJ is voltage 


H/Hm=(V+2AV —2AVy)/(V+2AN). (11) 
It will be shown later that the value of AV» is given by 
AVy=1.812-10-5T,,'-3[B,—B(8) |. (12) 


B, is given in Table IV. It is a function of m, the temperature exponent for 
the property JZ in question. For n>5, B, is given by the equation 
B,=0.5170+1. 1660 logig m»—0.0591/n+0.2224/n? 
+0.0140/n?—0.468/n*+---. (13) 
The coefficient of Eq. (12), 1.812-10-°7’,,!3, is given in Table V. 
B(@o) in Eq. (12) is a function of @. It is independent of » if 7 >4 and 
6,<0.5, and is given by 
n>4,00<0.5 B(Oo) =(x/a)oe. (14) 


This value of x/a is to be taken directly from Table II for @=6o. For n=1.2 
(the exponent for resistance and the only important small value of m) 
B(@o) is given in the third column of Table II. 
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TABLE V. 1.812-10-5T7,,!4 

















a | 1.812 -10-T,,!3 Tm 1.812 -10-5T,,!3 
1000° 0.1439 | 2300° 0.4250 
1100 1629 2400 4493 
1200 | 1825 2500 4737 
1300 2024 | 2600 4985 
1400 2229 2700 5236 
1500 | 2438 | 2800 5490 
1600 2651 | 2900 5744 
1700 2869 | 3000 6003 
1800 | 3091 ! 3100 .6266 
1900 3316 ] 3200 6528 
2000 3544 1 3300 .6797 
2100 3777 3400 7066 
2200 | 4011 3500 7337 

| | 





Since the current is constant, the voltage has the same temperature ex- 
ponent as the resistance. Hence AV may be found from Eq. (12) by setting 
n=1.2. In many cases to apply Eq. (11) it may be easier to find the theoreti- 
cal voltage V,,=V+2AV directly from the resistivity at 7, and the fila- 
ment dimensions. The wattage imput depends on the resistance of the 
filament and hence AVy in this case is to be found for n»=1.2. The wattage 
radiated on the other hand depends on m=5.1 or thereabouts." 

Method of application. In finding the value of B, for Eq. (12) from the 
value of m, we notice that for most properties of tungsten m is not constant 
as was assumed, but varies slightly with the temperature. We must take 
a mean value of , that is, its value at some effective temperature 7',. 
This temperature is roughly that at which h=h,,/2. This temperature 
and the corresponding value of may be found directly." 

As an alternate method to find 7, we note that for some properties 4 may 
be quite accurately expressed as 


h=CT*e-*/7. (15) 


Thus for candle power Wiens’ law (using a Crova wave-length) gives k =0, 
b=25200°. The Richardson-Dushman equation for the electron emission 
from pure tungsten has k=2, b=52600°. The rate of evaporation of a tung- 
sten filament is expressed by Eq. (15) with k=0, b=94100°. Setting h=h,,/2 
in Eq. (15) and using Eq. (6) to evaluate the term in k, we find approximately 


b/ Tr =b/Tn+[1—k/n] log, 2. (16) 
Differentiating Eq. (15) and comparing with Eq. (7) we see that 
n=k+b/T. (17) 


Hence from Eq. (16) and the values of the constants given above 
we obtain the following equations for effective values of m in terms of T,, 


candle power #=25200°/7,,+0.7 (18) 


% Ref. 3, p. 312, Table I, column 5. 
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electron emission 2 =52600°/7,+2.6 (19) 
evaporation 7 =94100°/7,,+0.7. (20) 


The decrease in the rate of evaporation near the leads is ordinarily not 
a matter of experimental interest, but under certain conditions its effects 
may be directly observed. Nitrogen or carbon monoxide in the presence of 
a tungsten filament at very high temperatures gradually disappears because 
every atom of tungsten which evaporates combines with a molecule of the 
gas to form a stable and non-volatile compound. Thus the rate of “clean- 
up” of the gas depends on the total amount of metal that evaporates. 

The direct application of Eq. (12) as outlined above is the most accurate 
method for the evaluation of AVy. In many cases where only approximate 
results are desired AV may be found from the following empirical equations, 
which were found to fit the data calculated from Eq. (12). The deviation 
from the results of Eq. (12) is less than the amount tabulated for the given 
range. The actual error of the results may in some cases be larger than this, 
due to approximations in the derivation of Eq. (12). 


Tasie VI. Alqg=P(T,,/1000) —Q(7To/1000) — R volts. 


IT is Pp Q R_ | Range—- Range— Range— |Max.2rror 
r, " T» | Vo | (volts) 
Voltage* | 0.154* 0. 0s1* | | 0. 056" | 1000 2500° any values 0. 1 0. 5 | 0. 008 
Candle Power 338 | .182 |—.004 | 600-3500 | 300-1400° | .1- .5 |: 
Electron emission 4400 | 6158 | .072 | 1000-3500 | 300- 900 | .07-.5 ‘009 
Evaporation | 480 .100 | .060 1500-3500 | 300- 900 | .07-.5 .008 


Watts radiated 293 | 160 | .084 lh inal 3000 | 300— 900 | .1- al 01 


* For voltage (Watts nie’ aterm —2.1-107 'T.Tm i is to “— added to the row —_ side 
of Eq. (21). 


In many cases with short, heavy leads T)=300°. In these circumstances the following 
approximate equations hold. 











TaBLe VII. Aly =P)(T,,/1000) —S volts (22 
Yn —— _ “em —- — ———— 
H is Po 5 om Te | | Max. error (volts) 
Voltage | 0. 148 | 0.080 1000- 2500 0.004 
Candle Power . 338 | .051 600-3500 | .O1 
Electron Emission .439 | .119 1000-3500 .007 
Evaporation | .477 .103 1500-3500 .001 
Watts radiated 287 121 1000-3100 | 01 


{ 








Computation of T». If we knew the diameter of a filament and the cur- 
rent through it, 7,, may be obtained directly from Tables which give tem- 
perature tabulated against current divided by d*/?.% 

If the diameter is not known, but if the length 2x is known, the volt- 
age V and amperage A corresponding to the temperature we wish may be 
found. Assuming that the filament is all at the maximum temperature 
Tm, we compute VA'/$/(2x) and find a first approximation for 7,.1° For 


% Ref. 3, p. 312, Table I, column 6 gives VA'/*/(2x) as a function of T. 
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this value of 7°, there is a certain voltage correction AV. This gives us a 
much better value, 1} +2AV, for the voltage if the filament were all at 7,,. 
(V+2AV) A'/*/(2x) then gives us a second approximation for T,,.. As many 
approximations as desired may be made. 

Shorter filaments. \With shorter filaments the cooling effects of the two 
leads overlap, and the temperature at any point may be found approximately 
by adding the cooling effect of each lead at that point. With still shorter 
filaments these temperatures and the values of J/////,, found as above 
are in error. The amount of the error depends on n, the temperature expo- 
nent of the property in question. Part II, Table X gives in column 2 the 
maximum value of the half length (x/a)( for which the error in H/H,, is 
less than 1 percent. Column 3 gives similar information for 5 percent error. 
For details see Part II. 


DERIVATION OF THE EQUATIONS 


The fundamental differential equation giving the temperature distribu- 
tion near a cooling lead is® 


A®rt [N(d?7/dx*) + (dd/dT)(dT/d x)? |xD?/4 = w. (23) 


The symbolism is explained in Table VIII. The terms A*r and w correspond 
respectively to the rate of production of energy and the rate of radiation 





TaB_e VIII. Symbols. 
A Filament current in amps w =h for power radiated 
D Filament diameter (cm) r =h for resistance 
Dr Lead diameter (cm) v =h for voltage drop 
Length of lead (cm) H Value of any property for the whole 
Distance along filament (cm) filament 
Unit of length (cm) Table I Ba Value of H if the whole filament 
dy a for D=0.01 cm were at 7), 
sub m Value at the uncooled central por- V =H for voltage drop 
tion of the filament HH. Value of H if the whole filament 
sub c Value at the center of the filament were at 7, 
(Part 11) n Temperature exponent for any 
a Absolute temperature property =d log hd log T 
Tr Room temperature p =n for resistance 
T, = )T,, Lead-filament junction temperature =» =n for radiation 
AT =1T,—Tr k =n for thermal conductivity 
0 =i /Tw ry =adé/dx 
(x/a ie: Distance in a-units from point at @ Aly see Eq. (12) 
to point at 4, B, Table IV 
r Thermal conductivity of filament Bis Value of B, for w=5 
XL Thermal conductivity of lead B(%) Table Il 
h Value of any property per cm of — 1% =T7,/T. 


filament length 


of energy per unit length of filament, while the expression involving A corre- 
sponds to the net rate of conduction of energy into an element of the 
filament. 
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An inspection of tables giving the characteristics of tungsten filaments 
as functions of the temperature shows that it is possible to express Eq. (23) 
in a much simpler form.’ 

The resistance of a tungsten filament can be expressed quite accurately 
over a wide range of temperatures by the equation r=c7*, where c isa con- 
stant and p=1.20 (for the range between 600°-3,000°K). Similarly the 
radiated power may be expressed approximately by the relation w=c’'T*, 
where w is fairly constant, having the values, 5.65 at 1,000°K, 5.12 at 1700°, 
4.93 at 2,000°, 4.71 at 2,400° and 4.48 at 3,000°. In the majority of experi- 
ments in which it is desired to calculate the cooling effect of the leads the 
temperature of the hottest part of the filament will probably be below 2,400°. 
By averaging the values of this exponent from 2,000° to 400°, weighing 
each in proportion to the corresponding value of w, the effective exponent 
is found to be 5.1. We shall, therefore, take this to be the value of w. Even 
at very high filament temperatures, where the effective exponent would be 
about 4.7, we shall see that the error made by using w=5.1 is practically 
negligible. 

The heat conductivity of tungsten at temperatures from 1,300 to 2,500° 
has been given by Forsythe and Worthing.’ It ranges from 0.93 watts 
cm! deg at 1,300° to 1.21 at 2,500°. We find that the empirical equation 


4=0.840(7/1000)°"4 (24) 


expresses the values of A at the 13 observed points given in their table 
(at 100° intervals) within an error of 0.0022 or about 0.2 percent. 
This equation is used throughout this paper. 

In the central uncooled portion w,,, the power radiated per unit length, 
is equal to A?r,,, where 7, is the resistance per unit length at this place. Since 
A is constant throughout the length of the filament, the temperature ex- 
ponent of A*r is the same as that of 7, that is p. Hence we can replace the 
first term in Eq. (23) by w,,0*. Similarly w may be replaced by w,,6°. 
From Eq. (24) we obtain the relation \=A,,0°*, where A,, is the thermal 
conductivity at temperature 7,,. Using the values of \ and dA\/dT from this 
relation, we obtain from Eq. (23) 


d6/dx?+0.4(d0/dx)?/0 = (92-9 4 —ge-9"4) /q? (25) 

where a is a parameter defined by 
@ = rD*\nTn/4wm. (26) 
We can replace w, by its value v»?/rm where vm is the voltage drop per cm 
at temperature 7,,. The factor 7,,D? which then occurs in the equation is 


independent of D and varies as T,,'? (it is in fact the function R’ given by 
Jones and Langmuir).* Thus 


’mD*? =7.89-10-°T,,!*?. (27) 
From Eqs. (24), (26), (27) 
a=1.812-10-57,,'°?/om. (28) 
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In computing a from this equation v,, was found from Tables of charac- 
teristics.% Since v, « D-'/?, we obtain Eq. (2). 
In Eq. (25) set w=5.1, p=1.2 and substitute @ as defined by Eq. (1) 


for x 
od¢/d0+0.4¢2/6 =04-7—9"°8, (29) 


By taking as new variables ¢? and log @ this equation becomes linear and may 
be solved in the usual way to give 


$? = (3 — 562-64 296-5) /6. 598, (30) 


The constant of integration, 3, is fixed by the condition that d6/dx =0, and 
therefore ¢=0 at @=1, the center of the filament. 
For values of @ close to unity (30) may be expanded in terms of Z =1—8 


¢? =3.9Z*—4.81Z%+3.715Z4—0.8886Z5+0.4068Z°+0.27Z7+ ---. (31) 


Temperature distribution. When @<0.6 the integral in Eq. (2) may be 
found by expanding 1/¢ in the series 


1/o=1.4720°'4[1+(5/6)02°8+(25/24)05? 
— (1/3)0%5+1.447078—(5/6)0%!+ --- |]. (32) 
Integration gives 
(x/a)§ =1.05146'-4-+0. 306764+0.23230°* 
—0.062167-°°+0.23109?—0.1179!5+4 ---. (33) 
When 020.6 Eq. (31) gives 
1/¢=0.5064[1+0.6167Z+0.0941Z?—0. 180928 
—0.227Z4—0.172Z5+ --- ]/Z. (34) 


Integrating we obtain the indefinite integral, and find the integration con- 
stant by comparison with Eq. (33) at 6=0.6, where both series are sufficiently 
convergent.to give results accurate to 1 part in 1,000 


(x/a)y? =0.2247 —1.1660 logio Z—0.3123Z—0.0238 Z? 


+0.0305Z*+0.0287Z4+0.0174Z5+ ---. (35) 


Column 2, Table II was calculated from Eqs. (33) and (35). 

Lead-filament junction temperature. The rate of flow of heat, Q, in watts, 
past any point of the filament is equal to the integral, from that part to 
the center, of the difference between the heat generated by resistance and 
that lost by radiation. 


Q= wae f (0!:?—95:')d0/o. (36) 
8 


Setting w,,=Av,, substituting for a from Eq. (28), and integrating 
Q=(1.812-10-57,,'°3)6. 5~°-5(3 — 5926+ 296 '5)°-5.4 (37) 


where @ corresponds to the temperature of the point in question. 
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For short fairly heavy leads at temperatures below 1,000°K the resist- 
ance loss and radiation loss in the lead are negligible. All the heat that flows 
into the lead from the filament must flow out at the cold end, which may be 
assumed to be at room temperature. Taking the heat conductivity of nickel 
leads as constant at 0.586 watts cm ' deg ', the leads have a constant tem- 
perature gradient given by 


AL(dT /dl)rD1?,/4=0. (38) 
Consequently for the total temperature difference A7’=7,)—7, for the 
leads we obtain 

AT =4Q/ (xD 7X1). (39) 
Eqs. (37) and (39) give Eq. (4b) for A7 in terms of (A7))o. To tind the latter 


we notice that the lead temperature @) has little effect on Q. Setting 4 
=(0.24, A=1,/=1, D,=0.1, A, =0.586, we obtain from Eqs. (37) and (39) 


(AT) o=145[1.812-10-57,,'3], (40) 
Eq. (40) was used in conjunction with Table V to compute Table ITI. 
Filament characteristics. We have ° 
vl 
24H =Hy—U=20 | (hn—h)d0/¢. (41) 
Oy 


Applying Eqs. (6) and (28) and substituting for h,,, v7, the equal ratio 17»; Vim, 
where V,, is the voltage drop that would exist between the ends of the 
filament if it were all at the temperature 7,,, 


1 
AV yn =AH(V'n/H») =1.812-10 cae | (1—0")d0/o. (42) 
80 


AVy is a convenient symbol for the expression AJ/(V,,/J/,,). It represents 
the voltage acrossasection of uncooled filament of such length that Z/,, 
for this section would equal the decrease caused by the cooling effect of the 
lead. The advantage of this nomenclature is that it requires no knowledge 
of filament diameter or length. 

By breaking up the integral of Eq. (42) into 


l 
B,= f (1—0")d0/o (43) 
0 
and 
99 
a(o.)= f (1040/6 (44) 
0 


we obtain Eq. (12). 
In evaluating B,, we meet the difficulty that we must use two different 
series for 1/@, one for small and the other for large values of 6. Let 1/¢, 
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be the value of 1/@ given by the first six terms only of the series in Eq. (32) 
and 1/¢» by the first six terms only of the series in Eq.(34). We have 


0 =<¢50.5 = >1 
0.5 <6< 0.65 1 = ¢ = ode 


Since the series are equivalent between 0.5 and 0.65 we may put 
t 1 
B,= f (1—0")d0/o,+ f (1—0")d6/¢z 
0 t 

where 0.5 </<0.65. A simple transformation gives 

t 1 

B, = f (1—0")(1/¢1—1/62)d0+ f (1—0")d0/¢z. (45) 
0 0 


The value of 1/¢,—1/@2 is*practically zero for values of 6 between 0.5 and 

0.65, and therefore the value of the first integral is not dependent on the 

actual value of ¢. Evenat@=0.1 the value of1/¢;—1/¢:2 is only 0.117, but 

at @=0 it becomes 0.573. Thus the larger part of the first integral is for values 

of @<0.1. In this range 6” may be neglected in comparison with unity for 

all large values of m. Even if »=1 the error in neglecting 6" will be small. 
Calling the first integral of Eq. (45) Fi, we have 


t t 
Fy = f dé, ‘b1 —_ J de, ‘be 
0 0 


t 
=(x ‘a)i— f d0/¢z. 
0 


Designate the right hand side of Eq. (35) by F2. From the derivation of this 
equation 


(46) 


F:+C= f d8/¢2 


C being a constant of integration. Hence 


z=l1—¢ 


t 
f d0/o2= [F,+Ch.= [Fo+C Jen1 =F,(1 —t) —F,(1). 
0 


But from Eq. (35) Fe (1-2) = (x, ‘a) >: Hence from Eq. (46) 
F,=F,.(1) = —0.0348 


the numerical value being found by setting s=1 in the expression for F2 
in Eq. (35). Due to the definition of ¢ this value is not affected by the missing 
terms of the series. 

Putting the value of 1/¢. from Eq. (34) in the second integral of Eq. 
(45), we find 
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6 1 

B,= —0.0348+0. 5064 Yay f (1—8")(1—6)?-*d6 
p=1 0 

where A;, Az, As are the coefficients, 1, 0.6167, etc., in the series in Eq. 

(34). This reduces to 


1 
By = 0.2247+0.5064 f (1—0")d6/(1—6) —0.3123/(n+1) 
0 


—0.0477[(u+1)(m+2) |-!+0.1833[(a+1) -- + (u+3)]-'--- (47) 


the coefficients of the next two terms being 0.6894 and 2.092. 
If 2 is an integer we have 


1 
f —onay/(t—0) =14 444444 +1/n (48) 
0 


and for other values of m the integral can be expressed in terms of gamma 
functions. For »=1.2, its value is 1.1216. For large values of m the integral 
is given by the series 


1 
f (1—0")d6/(1—0) =0.5772+log, n+1/2n 
0 


— [12n(m+1)]-'— [12n(m+1)(n+2)]-1. (49) 


Inserting this value in Eq. (47) and expressing as a series in reciprocal powers 
of n we obtain Eq. (13). 

In computing Table IV for values of m less than 5, Eq. (47) was used 
but Eq.(13) was found more convenient for the larger values. 

To evaluate B(6)) we note that in general 6)<0.5 and hence that for 
fairly large » the term @” is negligible, and Eq. (14) holds. For AV, when 
n=1.2, B(6) for the third column of Table II was obtained by using the 
series expansions of 1/¢ given by Eqs. (32) and (34). 

A pplication under other conditions. For filaments in the presence of gas, 
or for filaments of materials other than tungsten there will be changes in the 
values used in Eq. (29) for w, p and k, where k is the temperature exponent 
of the thermal conductivity. By methods similar to the derivation of Eq. 
(30) we find that the general expression for @ is 


go? = 26-**[(1— BP 1) /(p k+1)— (10°F H1)/(w+k+1)]. (30a) 


The integral B, in Eq. (12) may be found by numerical integration or by 
direct integration insome cases. Thus, when w+k+1=2(p+k+1), ¢? is 
a perfect square, and 


By =(w—p)-!2{ p[(n+k+1)/(o+k+1)]—v[(R4+1)/(o+k+1)]} — (50) 


where ¥(x) =dl/nI'x/dx is the logarithmic derivative of the gamma function. 
Table IVa gives some values of B, obtained in some of these ways. 
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; TABLE [Va. 
B, = fo (1—6")d0/@ for various exponents w, p, k, n 











@ p k n B, 

5.1 0.0 1.0 5.8 1.695 
4.3 1.85 —.4 20. 4.079 
6.0 1.0 1.0 1.0 0.371 
4.0 1.0 1.0 1.0 0.428 
4.0 1.0 1.0 5.0 1.118 
4.0 1.0 1.0 10.0 1.486 
4.0 1.0 1.0 20.0 1.871 
4.0 1.0 1.0 40.0 2.264 
3.8 1.2 0.4 27.2 2.522 








Data such as those in Table [Va show that the variation of B, with w 
is small and may be represented thus 
Forn=p B, = Bys[1+0.077(5.0—w) ] (51a) 
For n=20 Bi = Bys5[1+0.100(5.0—w) | (51b) 
where B,; is the value of B, forw=5.0. B,; may be found by these equations 
from tables of B, such as Table IVa. Fig. 1 is a plot of B,; as abscissa against 


p as ordinate for constant values of k. The full lines are for voltage cor- 
rection, =p. The dotted lines are for n =20. 


> 0.2 Oo “0.2 A 0.2 Q 0.2 
20 ; } ; ; ; 
/ / / / i 
2. r ‘ 
} { / ae , 
; / 
26 / / C6 ; Py 
/ / / / / 

e + / -f $ 4 ‘ ° oJ . ‘  é 

- / / J vA / - ’ / / / 
| / ; ; : ; / 


0.8 





10 44 48 22 26 30 aq 
8s 


Fig. 1. Plot of By (a =5) against p for constant k. 
Full lines » =p. Dotted lines » = 20. 

The deviations of B,; found from Fig. 1 from the value of B,, for the 
same n, found in Table IV may be expressed as a fraction, N, of the latter 
value. The variation of N with , for constant values of w, p and k may 
be approximately expressed by 








N(n) = N(20) —a(20—n) [N(20) — N(p) ]/(20—p) (51c) 
where 
nip| 5 | 10 | 30 | 40 
a{1| 0.42 | 0.27 | 0.15 | 0.11 
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We may find B, for any exponents n, w, k and p by finding N(20) and 
N(p) for the appropriate values of w, k and p from Fig. 1 and Eqs. (51a) 
and (51b). Eq. (51c) and Table IV then give the desired By. 

Other metals. From Eq. (26) and the derivation of Eq. (28) we see that 


@=(AmT mR m/Vm)*? (28a) 
4 


where R,, is the resistivity of the metal in question at the temperature 
T,. We then have 


AV in = (Xin TmRm)!/? [Bi —B(Oo) J. (12a) 


The Wiedemann-Franz law states that, at a given temperature, AR, 
and hence the coefficient in Eq. (12a), is approximately the same for most 
metals and alloys. This coefficient is thus given by 1.812-10-°7,,,) in Table 
V. Consequently the magnitude of AVy may be found for any metal which 
obeys this law with the help of Fig. 1 and the assumption that the change 
of B(@0) is similar to that of B). 

End losses from filaments in the presence of gas. Gas around a filament 
causes a loss of heat by conduction and increases the voltage required 
to reach a given 7. Eq (12a) shows that for a fixed T,,, AVy is the same for 
filaments in gas and vacuum except for the variation in the values of B, 
and 6(@)). The latter may be eva'uated by considering the conduction loss 
to be part of the radiation loss w in Eq. (23). Thus if at T,, the conduction 
loss is 1/4 of the radiation loss, and if the conduction loss varies as 7':7,!7 
the effect on B,; may be represented as a change in the effective value of w 
from 4.6 to 3.8, which by Table IVa means an increase of about 15 percent in 
the values of B,; given in Table IV. In general the values of AVy for vacuum 
hold with fair accuracy for small gas pressures, but the temperature distri- 
bution is altered. 

For new filament materials or other new conditions there will be no 
accurate knowledge of the filament characteristics and the temperature ex- 
ponents for the application of the above method, which method nevertheless 
will, we hope, still be capable of indicating whether or not lead losses are 
important in any given case. 


Part II. Leap LOssEs IN SHORT FILAMENTS 


Temperature distribution. Shorter filaments do not admit the assumption, 
made for most of the results of Part I, that the central portion of the fila- 
ment is not cooled by the leads. Even the calculation of temperature distri- 
bution by adding the effects of the two leads is not very good in these cases. 
Thus when x/a, the filament half length, is 1.53, for which the cooling effect 
of one lead at the other gives 6 = 0.994, the long filament case gives a central 
temperature corresponding to @=0.85, while the true value is @=0.8. Theo- 
retically the maximum temperature 7, can only be attained in an infinitely 


17 These conditions are approximately those for a filament of 0.007 cm diameter at 
T» =2900° in 10 mm of Ng. See I. Langmuir and G. M. J. Mackay, J. Amer. Chem, Soc. 36, 
1717 (1914). For larger filaments the conduction losses are relatively smaller. 
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long filament (cf. Eq. (35)). A “short filament” is one for which this fact 
invalidates the conclusions of Part I. 

Let 7. be the actual temperature at the center of the filament, and let 
the larger value 7,,, still indicate the temperature of a hypothetical portion 
uncooled by leads as calculated from the current and diameter of the fila- 
ment. Let 6.=7)./ Tm. 

Eq. (1) holds as before with a more complicated value for ¢. If (x6 is 
the distance from the center to a lead at temperature 7) =6@o7',,, we have, 
as before 


6, 
(x/a)e= f d0/ > (52) 
99 
a has the same values as before, and is given in Table I. 
Table IX gives the values of (x/a)( for various values of 6. Fig. 2 
gives the plot of (x/a)% as abscissa against @ as ordinate, for constant values 
of 6... The value of 6. for any curve is of course the intercept on the @ axis. 





Fig. 2. Plot of (x ‘aye against 8 for constant values of 0. 


These curves are temperature distribution curves. Thus if 6.=0.8 the curve 
with that intercept gives us the temperature (@) at any distance, (x/a)§, 
from the center. To find for any filament the value of 6., we need know only 
one point on the temperature distribution curve. Thus for a filament of 
known length and of known lead temperature, a point having the coordinates 
(x, ‘aye, 6) is determined on the plot. This point lies on some temperature 
distribution (constant 6.) curve—probably not one of these drawn. The 
value of @. may readily be found, however, by interpolating between the 
two nearest curves. By continuing this interpolation down to the x/a scale, 
the value (x/a){ is found. This is the value that would obtain if the given 
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filament were prolonged until its leads were at 0°K, and were left otherwise 
unchanged. If desired, (x/a)f may be found from @, by interpolating in 
Table IX. 

Of course the value of any property / at any point of the filament may be 
found by first finding the temperature at that point, then applying published 
data on filament characteristics.!4 


TABLE IX. Relation between the temperature T.=0-T» at the center of a short filament with 
leads at °K and the length (2x) of the filament. 


























0. | (x/a)% 6. (x/a)% 6. (x/a)% 
0.005 | 0.7325 0.4 1.2065 || 0.94 2.031 
01 .8262 5 1.2510 | .95 2.118 
.03 9221 || 6 1.3077 | .96 2.225 
.05 9704 | E 1.3905 | .97 2.365 
A 1.0401 8 1.5280 | .98 2.565 
2 1.1154 .85 1.6399 | .99 2.912 
3 1.1645 9 1.802 995 | 3.261 
92 1.899 999 | = 4.074 








The effect of lead losses on characteristics. If all the filament were at the 
actual maximum temperature, 7.., the value H. of any property for the 
whole filament would be 


H, =2h.x = 2ah(x/a)§. (53) 


h, is the value of the property in question for 1 cm of filament at temperature 
T.. The ratio of the actual value for the whole filament, J, to the hypo- 
thetical value is H/H.. The value of this ratio when 6)=0 we designate 
as (11/H.)o. If n>4, H is independent of 65 for constant 6. for all the prac- 
tical range. Hence, applying Eq. (53) 


H =2ah,(x/a)((H/H.-)o. (54) 


H, and hence H//H., is a function of n, the temperature exponent of the 
property in question. In Fig. 3, ordinates (H/H-.)o are plotted against 
abscissae (x/a)¢ for constant values of m (the full lines). The scale for 
(x/a)ge at the top, and the scale for @, on the bottom may be used inter- 
changeably. They correspond as in Table IX. 

For a given filament, knowing x, a, and 0, (x/a)¢ is determined from 
Fig. 2. Fig. 3 then yields (J7/H-.)o for the value of ” corresponding to the pro- 
perty in question. H may then be found from Eq. (54). 

For n $4, H is not independent of 0. For =1.2, the resistance-exponent, 
the dotted lines at the top of Fig. 3 give the value of H/H.[not (H/H.)o| 
for constant values of to=60/8.._ Note however that this is still given in 
terms of the abscissa (x/a)¢. 

For values of x/a less than 1.0, and hence not in Fig. 2, H/H. is constant 
at the value for x/a=1.0. 

It is to be remarked that for 7.<1000°K none of the given data apply 
accurately, nor is the temperature distribution accurate. This is because 
of the uncertainty in the value of J here. 
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Limits of long filament case. Table X gives the smallest value of 
(x/a)¢ for which the errors in 17/H, made in applying the long filament 
instead of the short filament case are less than the percentage at the head of 
the column. The errors of course depend somewhat onthe value of m for 
the property under consideration. A (+) after the value of indicates that 
the long filament case gives too large a value of H7/H., or of H; that is, too 
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Fig. 3. Full lines (47/H,), for constant m against (x, a)? and 6... Dotted lines 
H/H, for n=1.2 constant 7. against (x/a)ée and 6. 


small a value of B,. In the neighborhood of »=5 the sign of the error 
changes. For n=1 (temperature distribution) the long filament case gives 
too high temperatures, even when the cooling effects of the two leads are 
added. 

Lead temperatures may be found with sufficient accuracy for most cases 
from Table III of Part I. For high accuracy or for very short filaments 
Eqs. (36) and (39) may be applied directly, using Eq. (55). 

TaBLE X. Errors in H/H,. made in assuming a filament to come under Part I. This 


able gives approximately the smallest value of (x/a)o® for which the error made by that assumption 
is less than the percentage given at the head of the column. 











n 1 percent 5 percent n 1 percent 5 percent 
1.0(+) 1.9 1.6 10(—) 2.8 2.3 
1.2(+) 1.9 1.7 15(—) 2.9 2.6 

4(+) 1.9 1.8 20(—) 3.0 2.7 
5(+) 1.6 1.6 30(—) 3.2 3.0 
7(-) 2.3 2.0 40(—) 3.3 3.0 
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Computation of temperature distribution. Eq. (30) holds with « new con- 
stant of integration determined by the condition @=adé dx =0 when 6=8@, 
(the maximum temperature): 


1.3¢2=60 "(0,2 5-0. 40,65—62 +0. 4985). (55) 


The integration in Eq. (52) is best carried out by series approximations. 
Case 1. r=0 6. is small 


t 
(x/a)h =6, f dr/o. (56) 
0 


Expressing Eq. (55) in terms of 7, then expanding 1,¢ by the binomial 
theorem and integrating the resulting series term by term 


(x/a)o=0. 81444-9907 4[1 0.175472 
+0.0795.42752—0.0354.10,2 97854 --- | (57) 
where A =1/(1—0.4 6.5-°). 


In Eq. (57) 7°? has been neglected. If @<0.60. the error is less than 1 
percent. If in this series we set 0.=1, we get the series of Eq. (33) of Part I. 
Case II. c=1—7T=(0.—86)/0. is small and @, is also small. 

Eq. (55) may be expanded in terms of ¢ to give 
o? = 20,'86(1 —8,3-9) (1 —a)~"'8[1—0. 80(1 — 2. 4380,3°9) 

+0.1607(1— 24. 780,3°9) +0.01603(1+4+2290,3°9)+ --- |. 
Terms of the form (1—06,3-°)/(1—63-°) were expanded by the binomial 


theorem, and @.7-° and higher powers neglected. Expanding 1/¢ and inte- 
grating 


(x/a)ee=0, f do/d =." (1+ 40,3"9)(2¢)°"5[1 —0. 330.386 
0 


—0.02402(1—100,3°*) —0.00340°(1+5.380.5°)]. (58) 


This is accurate when o* and 6.7: may be neglected in comparison with 
unity. At @.=0.56 this error is about 1 percent. 
Case III. o is small and 8, large. 

Let z=1—6 and z.=1-—8@,. z and zg, are both small. Expanding (55) by 
the binomial theorem 


1/p=0.5064(1—0.42—0. 1222+ - - -)(s?—2,2)-°5[1+1.0167Ci+Ce24+ - - -]. (59) 
where 
Ci=(s8—2.8)/(2?—2.2) =2[1+(1/2)(z./z)?] 
C2=1.5504C,?—0.930(s4—2z.4)/(2?—2.2) 
= [0.621+1.008(z,/z)?]z? 
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substituting, multiplying Eq. (59) out, collecting like terms and integrating 


(x/a)je = [Pave =().5064(1+0.852z.*) arc cosh (z/z,) +0. 3123(22—2,?)""5 


+0.2574s. arc cos (./z) +0.02392(2?—2,?)5 (60) 


This series is good if 6.2 —0.7; 620.5. 
Case IV. @. is nearly unity. 

We may choose z small enough for series (60) to converge, but at the same 
time much larger than z.. The small quantity z, then has little effect on the 
integral from 0 to 1—z. Toevaluate this integral we may neglect z, entirely 
and use Eq. (35) of the long filament case, which, by the choice of the integra- 
tion constant, gives the value of (x/a)’. z as chosen above is sufficiently 
small for this series also to converge. Adding Eqs. (35) and (60), expressing 
arc cosh (z/z.) as a logarithm and arc cos (z,/z) as a series, and remembering 
that z is much larger than z,., we obtain 


(x/a)o° =0.5757 —1. 16596 logis z-+0.40432,+1.92.. (61) 


The coefficient of the last term was chosen empirically so as to compensate 
for the missing terms. Note that the result is independent of the specific 
value of zs, as long as it has such a value as to make the derivation valid. 

Thus we can obtain (x/a)’ by Case I and (x/a)§ by Case II or III. 
These two methods overlap in the admissible values of @ except for inter- 
mediate values of 9., when neither Case II nor Case III is very good. Numerical 
integration was used for accurate results in this region. The value of (x/a)8 
may be obtained by Case IV or by the formula 


(x/a)§e=(x/a)o+ (x/a)pe , 


Table VII and Fig. 2 were obtained by the above methods. 
Filament characteristics. Assuming as in Part I that the value of the 
property in question varies as the mth power of the temperature, we have 


H=2h. f (T/T.)"dx =2ahJ (62) 
0 
where 


Be 89 
J =(1/0.") f 6"d8/¢—(1/0.") f 6"d0/¢. (63) 
0 0 


If n>4, 6" becomes rapidly smaller with decreasing 6. Hence we may 
neglect the second integral to obtain 


Be 
s=(1/0e) f 6"d0/o. (64) 
0 
From Eqs. (53) and (62) 
(H/H.)o=J/(x/a)t. (65) 


There are two limiting cases to be considered. 
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Case I. 6. is small. 
We may neglect 6.65 in Eq. (55). Introducing 7 =0/0. we obtain from 
Eqs. (63), (53), (62), (52), and (55) 


1 
rrt0-4d7( 1 — 72°8)—0°5 


H/H,=— 


7°-4dzr(1—7?°8)—0°5 





To 


The value of 6. has cancelled out. If to=00/0. is held fixed, then H/H. is 
independent of 6... Hence in Fig. 3 the values for x/a<1.0 are the same as 
those given for x/a = 1.0. 

Case II. @. is large. 

As 0. approaches unity we get the transition from the short filament to 
the long filament case. When z.=1—6. can be neglected entirely, Eqs. 
(64) and (43) give 

(x/a)o¢—0."J = By. (66) 
From Eq. (62) 

H =2ah(0."J). 
Thus when (66) is satisfied the value of 7 increases in direct proportion to 
the increase of length. This is because the central portion, to which the 
increased length is added, is at practically constant temperature. 

The degree to which Eq. (66) is satisfied is a measure of the approximation 
involved in assuming that the filament is long. To construct Table VIII the 
true values of J as found below were compared with the value calculated 
from Eq. (66). For n=0 the short filament temperature distributions were 


compared with those obtained from the long filament case by adding the 


cooling effect of the two leads. 
Evaluation of J. If @. and hence D=0.46,3-° (1—0.46,3-°)— is small, we set 
y= (0/0.)?* in Eq. (64) and obtain, using Eq. (55) 


J =0,9'M(1/5.2)*/2(1—0. 40,3°%)—"/2 (67) 


where 
1 
u=f yPd y(1— y)-¥?[1—Dy(1 — y*")(1— y) 71? (68) 
0 


and »=2.6p+1.2. M may be expanded as a power series in D 
M =Go+ (1/2)G,D+ (3/8)G2D?+ (5/16)G3D*+ iB doe (69) 


where 


l 
G,.= f yPt *(1 — y)-¥2[(1 — y8/2)(1— y)-*] Ady. (70) 
0 
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Each of these integrals may be expanded in terms of 1=1—y, and then 
integrated by means of gamma functions. 





Go=2(—1/2)4(p)/x(p+1/2) (71) 
3x(— 1/2)x(pt+ 1) 
G,= 1—(1/8 5/2)-! 
nlp +3/D [1—(1/8)(p+5/2) 


—(1/32)(p+5/2)-"(p+7/2)-1+ - +> | 
a reret [1—(1/4)(p+7/2)" 
—(1/64)(p+7/2)-"(p+9/2)-1+ +--+ | 
_ 2in(—1/2)4(p+3) 
— 8x(p+7/2) 
+ (3/64)(p+9/2)-"(p+11/2)-'+ --- J. 








[1—(3/8)(p+9/2)“ 


As 6. and D become larger the series of Eq. (69) does not converge rapidly. 
A better series may be obtained by setting 


M =G,(1—giD—g.D*— - - - )-*/2 (72) 


and determining the value of each g; in terms of the G;’s by expanding Eq. 
(72) and equating coefficients with Eq. (69). 

If 6. is very small, D may be neglected entirely. For the case n=0, 
Eqs. (69), (71), (67) give 


J =(«/a)§=1.309(0,)°. (73) 


This is useful in determining the smaller values of Table VII. 

When 0.9<6@,<1 even series (72) does not give accurate results. The 
values of J here were determined by means of Simpson's rule. The difficulties 
due to the infinite integrand at @ =8@, were avoided by an integration by parts. 

If m is small, 8, may not be neglected as in Eq. (64). The most important 
case is m=1.2 (the exponent for resistance). The values of J may be found 
from Eq. (63), the first integral being evaluated by the methods above, 
and the second integral by a series similar to that of Eq. (57). This holds 
over the useful range 0)<0.6 @.. 

A similar method might be used for other small values of n. The calculations 
for the second integral may be simplified by neglecting 6°*° in Eq. (55). 
At 6)=0.606. the error made thus is less than 3 percent for any short filament 
(9.0.94). Since this integral is a correction term, the approximation is 
justified. By substituting “=1—0.4603-9—(6/6.)?*, the integral is reduced 
to a form which may be evaluated in terms of elementary functions for 
n=1.2, 2.5, 3.8, 5.1 and 6.4. Interpolation may be used for intermediate 
values of . 

For n24 Eq. (64) is valid, unless high accuracy is desired. When n=10, 
6. =0.7, 09 =0.5, the error is less than 1 percent. 
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Experimental checks. Several lamps were made up to test the short fila- 
ment theory by experiment. Filament L was a straight filament 1.6 cm 
long, D=0.0256cm, cut froma length of wire that had been aged 24 hours at 
2400°. It was welded to nickel leads 5 cm long, D; =0.254 cm. The tempera- 
ture distribution along this filament was measured by means of an optical 
pyrometer mounted on a carriage that could be moved along a horizontal 
scale which gave the position of the pyrometer accurately to 0.001 inch. 
The pyrometer had previously been calibrated against a standard lamp. 

The curves in Fig. 4 show the temperature distribution for three different 
currents. The solid lines indicate the observed values, and a comparison of 
7. in each case with the corresponding value of 7’, listed in the corner of 
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Fig. 4. Temperature distribution for filament L, D=0.0256 cm, 2x=1.6 cm, 
wire aged throughout its length. Experimental and theoretical curves. 


the graph, shows that the center of the filament is greatly cooled by the 
leads; in other words, filament L is a “short filament.” 

The temperature distribution was calculated from the short filament 
theory and the results shown in the curves marked by triangles. In Curve 
1 where 0,=0.912 and 7,—T,=234°, the agreement with experiment is 
fairly good. In the more nearly extreme cases of Curve 2 where 0. =0.836, 
T»— 7T.=408°, and Curve 3 where 6, =0.687, 7,,—7.=725°, the tempera- 
tures given by the theory are too high. We attribute this departure from the 
observed values to an error in the value of \ at the cool ends of the filament. 
The heat conductivity of tungsten is not accurately known at temperatures 
below incandescence, and if we have used values of \ in this range that are too 
low, so that the calculated temperature gradient near the leads is too steep, 
the resulting 7, will be too high. An error from this source becomes important 
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only when the short filament theory is put to the severe test of predicting 
T.. in cases where 7, is much less than 7°». 

Since 7, must be known before either the voltage or the candle-power 
of a filament can be calculated, the following empirical addition to the theory 
has been devised, to be used in those cases in which 7, can be determined only 
by calculation. To compensate for the error that results from using values 
of \ that are too low, we proceed as if the length of the filament were shortened 
by an amount Ax such that the heat loss by conduction will be increased 
above Q,, the value corresponding to A, by an amount ~=4Q,Ax/(7D?*). 
Qy is given by Eq. (37) which for 69 =0.24 becomes approximately 


Q, = 0.6654 A 0.(1.812-10-7,,! 4) (74) 


It has already been pointed out, in connection with the derivation of Eq. 
(40), that changes in 4) have little effect on the factor 0.6654. 6, supplies 
approximately the factor by which Q, must be reduced when the integration 
in Eq. (36) is carried only to 6, instead of to 1. The term in brackets is 
given in Tab'e V. Values of ¥ derived from the data of Fig. 4 are given in 
Table XI. They were calculated using those values of Ax which made the 
theoretical and observed curves coincide at 1500°; the theoretical curves are 
indicated by crosses. yW is tabulated as a function of Jo, since the error is 
greater at low lead temperatures. 





TABLE XI. 
To= 300° 400° 500° 600° 
y= 471 367 263 159 


From Table V and Table XI one can calculate 


Ax = rD>p/(4Q)) (7 


wn 


) 


and subtract Ax from the actual half length x of the filament before calcu- 
lating 7... This correction applies to calculations of 7. only, and the maximum 
values to be used are 0.15x for leads in air, and 0.22x for leads in liquid air. 

Fig. 5 is a plot of the volt-ampere characteristics of filament G, 1.928 cm 
in length, D=0.0103 cm. The curves labelled air, using the bottom scale 
for voltage are for the bulb at room temperature, 300°K. The curves labelled 
liquid air, using the top scale for voltage, are for the bulb immersed in liquid 
air. 

The course of the liquid air observations for low voltages is interesting. 
For 450°<7°..<1200° and 7,,,=2,000°, the current decreases with increasing 
voltage. This phenomenon has been obtained with all short filaments which 
have been tried in liquid air. 

Thus for one value of the current there are in some cases three possible 
values of the central temperature of the filament. With a low temperature 
the heat generated is small and so the small temperature gradient is suffi- 
cient to carry away the heat and maintain equilibrium. Likewise, with a 
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higher central temperature the heat generated is greater and the larger 
temperature gradient is necessary to preserve stability. Hence it is possible 
that all three central temperatures may be stable. The phenomenon may be 
explained in more detail by assuming that below 1,200° X is larger than the 
value given by Eq. (24). Analysis shows that the general form of Eq. (73), 
which gives (x/a)* for small @., is 


(x /a) 2 8.0489) /2, (73a) 


If k, the thermal conductivity exponent, is less than 0.2, then by Eq. 
(73a) the temperature distribution curves of Fig. 2 will cross near @=0 and 
x/a=1. Thus (x/a)’ and hence also 7,, and A decrease with increasing 6. 
for small values of the latter. This is essentially the phenomenon observed 
above. 
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Fig. 5. Volt-ampere characteristics of filament G, D=0.0103 cm, 2x =1.928 cm, wire aged 
throughout its length. Experimental and theoretical curves. 


To obtain the calculated voltages, marked by crosses on the curves in 
Fig. 5, a correction of 0.6Ax subtracted from x was found sufficient to com- 
pensate for the decreased voltage drop along the ends of the filament due to 
the lower temperature which we assume to exist in this region. Thus for 
n=0(, (temperature distribution) the correction is Ax; for n =1.2 it is 0.6Ax; 
and for higher values of 7 no correction is needed. 

Example of the calculations. To illustrate the method of calculating lead 
losses, consider filament G running at a current of 1.295 amps. From 
D=0.01030 cm we find D*/? =0.001045, A /D*/? =1239, thence T,,, =2222°. 
Table I gives by interpolation ay=0.400, whence, as D!/?=0.1015 we find 
from Eq. (2) that a=0.406. To determine the lead junction temperature, 
we find in Table III that (A7))=59° and hence with /=5 cm and D,=0.254 
cm Eq. (4b) gives AT =59°. Thus 7)>=359°, and by interpolation in Table 
XI y=410. Table IV gives (1.812-10-57’,,!*) =0.4064, so that from Eq. (75) 
we have Ax =0.100 cm if we put 6,.=1. The half-length x=0.964, therefore 
x’=x—Ax=0.864. Dividing by a=0.406 we have (x’/a)f=2.128. Since 
To =359° we find 4) by dividing by 7, and have 6)9=0.1616. We then find 
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the point of coordinates 2.128, 0.1616 on Fig. 2, and drawing through it a 
curve parallel to the given curves we obtain the intercepts 6,=0.959'* and 
(x’/a)* =2.220. Multiplying 6. by Tm we find T,=2131°. 

Eq. (54) shows that the voltage of a short filament is given by 


V =A(H/H.)8xp./(aD*) (76) 


p. is the resistivity at 7. and is 61.12-10-* ohm cm!'*. To find H/H, we must 
first obtain To by dividing 359° by T. to obtain 0.168. From Fig. 3 for abscissa 
(x’/a)¢* =2.220 we find H/H.=0.766. We substitute for x in Eq. (76) the 
corrected value x—0.6Ax =0.904 and obtain V=1.315 volts. The experi- 
mental value was V = 1.330. 

To find the candle power we use Eq. (54) to give 


Candle power = H, = 2aDC,.'(H/H.)o(x/a)® (77) 


C.’, the specific candle power of tungsten at 7., is found to be 47.3 inter- 
national candles cm~?"'*). From Eq. (18) using 7. instead of 7, we find the 
effective value of m to be 12.53. From Fig. 3 for abscissa (x’/a)§* =2.220 we 
then find (H/JI.)9=0.362. (x/a)* is obtained from (x’/a)% by multiplying 
by x/(x—Ax) to give 2.477. We thus obtain 7, =0.355 international candles. 
We did not measure the candle power in this instance, but in other instances 
similar calculations of candle power gave results that were in good agreement 
with experimental values. 

Discussion of experimental checks. In ordinary practice the ends of a 
filament are “unaged,” that is, they can never be heated to the temperature 
which a filament must once have before the properties become those of 
“aged” tungsten. It is known® that heating tungsten wire for one minute to 
a temperature of about 1500° causes the cold resistance to be lowered 15 to 
20 percent, and presumably the heat conductivity would be increased at 
the same time by about the same amount. The filaments used for experi- 
mental checks were made from wire that had been previously aged through- 
out its length, in order to obtain uniform results, for the properties of unaged 
wire vary from sample to sample and depend on the history of the filament. 
It has been shown that in the case of aged tungsten the heat conductivity 
at the cool ends of the filament is higher than the values used in the theory. 
But where the filament is unaged in this region, the heat conductivity is 
lower than in the case of the aged filament which tends to restore A to the 
values used in the theory so that in general the short filament theory may 
be used without the correction Ax to calculate lead losses from filaments 
that are made in the ordinary way. 


18 Where @. comes out to be less than 0.95 it is well to substitute @. equal to the calculated 
value instead of unity in Eq. (74), and thus arrive at a second approximation for Ax from Eq. 


(75). 
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ABSTRACT 

The variation with the electron velocity of the absorption coefficient, a, or the 
effective collision cross-section has been measured in the vapors of cadmium and 
zinc and found to follow a curve of the type previously found for mercury. The 
cadmium curve has a maximum at about 40 volts, a2=130, and a minimum at about 
25 volts, a=126, followed by a steady rise with decreasing velocity to the limit of 
accurate measurements, 1 volt, a=700. The zinc curve has a very flat maximum at 
about 50 volts, a=76, a minimum at about 36 volts, a=74, followed by a steady rise to 
a=500 at 1 volt. The magnitudes of the maximums are in the order: Cd, a=130; 
Zn, a=75; and Hg, a=60. Other related properties of these atoms, such as the 
molar refractivity and the critical potentials, show this same irregularity in order. 


HE measurement of the absorption coefficient for electrons in the noble 

gases, Ar, Kr and Xe,' has shown that the shape of the curves obtained 
is the same for all three of these. Similar measurements in the vapors of the 
alkali metals? have shown that these elements also have a characteristic type 
of curve. Mercury has been studied by Brode,*:?7 Maxwell,4 Beuthe,> and 
Jones.® All of these observers except Beuthe have found a curve with values 
of the absorption coefficient in excellent agreement. Cadmium and zinc 
were also studied at the same time with the first measurements of mercury.* 
These observations were shown by subsequent experiments’ to be unreliable, 
due to secondary and reflected electrons which were present in the form of 
apparatus used for the first measurements. 

The apparatus, Fig. 1, used for the measurements of cadmium and zinc 
was of the same design as that described in the measurements of mercury.’ 
About half of the measurements in cadmium vapor were made with the ap- 
paratus used for the mercury measurements. The mean tfadius of the path 
was 15 mm and the final slit was 0.5 mm wide. The rest of the observations 
in cadmium were made with an apparatus in which the mean radius of the 
path was 10 mm and the width of the final slit 1.0 mm. All of the measure- 
ments in zinc were made with this same apparatus. 


'C. Ramsauer, Ann. d. Physik 72, 345 (1923). 

2 R. B. Brode, Phys. Rev. 34, 673 (1929). 

3 R. B. Brode, Roy. Soc. Proc. A109, 397 (1925). 

41. R. Maxwell, Proc. Nat. Acad. Sci. 12, 509 (1926). 
5H. Beuthe, Ann. d. Physik 84, 949 (1927). 

5’ T. J. Jones, Phys. Rev. 32, 459 (1928). 

7 R. B. Brode, Roy. Soc. Proc. A125, 134 (1929). 
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The control of the temperatures, the neutralization of the earth’s field, 
the deflecting magnetic field, and the arrangement of the electrical circuits 
was the same as previously described in the mercury’ and alkali? metal 





Fig. 1. Diagram of apparatus. 


papers. The pressure of the gas was obtained from the vapor pressure at 
the temperature of the lower furnace which contained the supply of metal. 
This pressure was corrected for thermal effusion and the observations were 
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Fig. 2. The absorption coefficient a for electrons in cadmium vapor as a 
function of the velocity of the electrons. 


reduced to 0°C and 1 mm of Hg pressure. The constants of the vapor pressure 
equation were taken from those given in the International Critical Tables.* 


8 International Critical Tables, Vol. III, p. 205. 
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These constants for zinc and cadmium are much more accurately determined 
than those for the alkali metals. The principle source of uncertainty in the 
values of the absorption coefficients is not due to the vapor pressure con- 
stants in this case but to the determination of the temperature. From the 
uncertainty in the temperature, not over 2°C, the resulting absorption 
coefficients might be in error by about 5 percent. Small differences in the 
temperature could be measured to 0.1°C. 

Samples of the metals were used which were over 99 percent pure. 
These were distilled several times in vacuum and finally into the apparatus, 
which had previously been baked at 500°C for several hours, and the metal 
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Fig. 3. The absorption coefficient a for electrons in zinc vapor as a 
function of the velocity of the electrons. 


parts glowed by an induction furnace. This treatment was effective in re- 
moving the gases absorbed in the cadmium and zinc as measurements of 
the absorption coefficient were constant over a period of two months in a 
closed apparatus. 

By plotting the log of the ratio of the current at the end of the path to 
the current from the slit in the cylinder, S., against the pressure, a straight 
line was obtained. The slope of the line gave the value of the absorption 
coefficient for unit pressure. Figs. 2 and 3 show the resulting values of the 
absorption coefficient for cadmium and zinc. Each point is the result of 
measurements at from 3 to 5 different pressures. Due to increasing leakage 
currents above 300°C, the data obtained with zinc were not as consistent as 
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those obtained with cadmium. The data taken with the two different exper- 
imental arrangements in cadmium are shown on the graph and are seen to 
be in excellent agreement. In the upper right portion of each graph the low 
velocity measurements are plotted. In cadmium the measurements were 
taken to 0.5 volts but, due to the leakage currents at the higher temperature 
used with zinc, the measurements could not be extended to 1.0 volt in zinc. 
There is an indication that the absorption coefficient in cadmium ceases to 
rise as rapidly below one volt as above one volt. 

The opening of the final slit determines the maximum angle through which 
an electron may be deflected without being measured as absorbed. In 
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Fig. 4. The absorption coefficient @ for electrons in cadmium, zinc and 
mercury vapors as a function of the velocities of the electrons. 


cadmium two systems of slits and chambers were used with different limiting 
angles but the values of the absorption coefficient were not changed. Meas- 
urements by several observers on mercury’ gave results that are quite con- 
sistent although the chambers used were geometrically quite different. 
This would indicate that the number of electrons scattered through small 
angles is only a small fraction of the total number of scattered electrons. 
These results are not in agreement with the measurements of Arnot® which 
indicate that a large fraction of the electrons of 82 volts velocity in mercury 
vapor are scattered in small angles. An apparatus with an adjustable 
aperture is being constructed to test this point. 


® Arnot, Roy. Soc. Proc. A125, 660 (1929). 
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In Fig. 4 the curves for cadmium and zinc are compared with that for 
mercury. The values of a for mercury’ have been changed from those 
due to vapor pressure data in Landolt and Bérnstein’s tables used in the 


original report to the data given in the International Critical Tables!® from ° 


which the vapor pressure data for cadmium and zinc were also taken. This 
involves an increase of about 10 percent in the values of a. The character of 
all the three curves is the same; a very large effective cross-section for slow 
electrons which decreases rapidly with increasing velocity to a minimum, 
then rises slightly to a faint maximum followed by a steady decrease with 
further increase in velocity. 

The curves indicate that cadmium is the largest of the three, zinc is 
next, and mercury is the smallest. This is the order of size to be expected 
rather than that of increasing atomic weight. The classical scattering of 
electrons due to the polarization of the atom by a passing electron was shown 
by Zwicky" to give the right order of magnitude for the absorption coeffi- 
cients in the noble gases at velocities above the maximum in the curves. 
The absorption coefficient was proportional to the molar refractivity. The 
molar refractivities for these elements are Cd = 20.0 Zn = 14.6, and Hg = 13.7. 
This same agreement between the absorption coefficients and molar refractiv- 
ity is found in the alkali vapors. The effective size of an atom is also indicated 
by its ionization potential. Helium, with the highest ionization potential, 
has one of the smallest effective cross-sections of any observed atom while 
caesium, with the lowest ionization potential, has the largest effective cross- 
section. The ionization potentials of the elements studied in this paper are 
Cd 8.9 volts, Zn 9.4 volts and Hg 10.4 volts. In the curves of elements not in 
the same row of the periodic table and therefore of different shape, it is 
not always easy to designate one as larger than the other. In general, how- 
ever, it seems to be true that the magnitude of the absorption coefficient or 
the effective collision cross-section of an atom is proportional to the molar 
refractivity and inversely proportional to the ionization potential. 


10 International Critical Tables, Vol. III, p. 206. 
uF, Zwicky, Phys. Zeits. 24, 171 (1923). 
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ABSTRACT 


The theory of metallic structure, of Sommerfeld, Heisenberg, and Bloch, is 
carried far enough to explain cohesive forces, and calculations are made for atoms 
with one valence electron, particularly metallic sodium. The numerical results, 
though rough, are in qualitative agreement with experiment. It is found that the 
forces in general are of the same nature as those met in ordinary homopolar binding, 
discussed by Heitler and London; except that the purely electrostatic force from 
penetration of one atom by another is relatively more important, the valence effect 
from the exchange of electrons relatively less important, than in diatomic molecules. 

As a preliminary to the calculation, the relations of the methods of Heisenberg 
and of Bloch are discussed, and it is shown that they are essentially equivalent in their 
results when properly handled. Remarks are made both about conductivity and fer- 
romagnetism. In connection with conduction, it is shown that a definite meaning can 
be given to free electrons, that they are necessary to conduction, and that a method 
can be set up for computing their number, which is rather small compared with 
the number of atoms. Ferromagnetism is discussed in connection with a recent 
paper of Bloch. It is shown that a metal like an alkali cannot be ferromagnetic, for 
atoms at such a distance that the interatomic forces keep the metal in equilibrium, 
are too close to be magnetic. For ferromagnetism, rather, it seems necessary to 
have one group of electrons responsible for cohesion, and another group, of smaller 
orbit and therefore relatively farther apart, producing the magnetism; a situation 
actually found only in the iron group and the similar groups. 


I. INTRODUCTION 


CRYSTAL of a metal is an enormous molecule, with electronic energy 

levels depending on the positions of all the nuclei, just as the electronic 
energy of a diatomic molecule depends on the internuclear distance. In this 
paper, in which we are interested in cohesive forces, we must find this energy 
of the lowest state in terms of the size of the crystal. We limit ourselves to 
geometrically similar arrangements of the nuclei, with changing scale. From 
the minimum of the curve, we find the heat of dissociation, grating space, 
and compressibility of the metal. But also we can investigate the wave func- 
tion of this lowest state, and obtain information about the electric and mag- 
netic properties of the metal. In this way we are naturally led to a discussion 
of the calculations of Heisenberg! and of Bloch on these subjects; in order to 
be sure that we really understand the arrangement of energy levels, we discuss 
the relationships of their methods, and arrive at a consistent picture combin- 
ing them. 


1 W. Heisenberg, Zeits. f. Physik 49, 619 (1928); 
F. Bloch, ibid. 52, 555 (1929). 
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As for the results, one naturally asks first, what are the forces holding a 
metal together? Are they ordinary attractions on account of penetration of 
atoms, or valence forces, or electrostatic forces of ionic attraction, or van der 
Waals forces, or some special sort not found in other cases? This question 
cannot be answered categorically; no doubt all the forces are simultaneously 
present, and the problem is to find the relative magnitudes. The tentative 
result at which we arrive is that the simple penetration of one atom by 
another is the most important part of the effect. But valence effects are also 
present, although weakened by having the valences shared by many neigh- 
bors, and are responsible for a considerable fraction of the attraction. 
Although these actual magnitudes may not be verified by more accurate 
calculation, still we have discussed the problem in enough detail so that the 
general relations can be understood in any case. 

The other question one will ask is, what is the situation of the electrons 
in the metal? Can one give a meaning to the question, how many free elec- 
trons are there? The answer, from whichever side we look at the question, 
seems to be the same. Most of the valence electrons are at any time attached 
to their atoms. These electrons cannot take part in conduction; they could 
do it only by having a whole file of such electrons simultaneously jump to 
the next atom in line, a most unlikely occurrence. But a few electrons at 
any time—calculation suggests a few percent —will be detached from their 
atoms, leaving an equal number of positive ions behind them; and they are 
what, by all rights, one should call free electrons. These electrons, and the 
positive ions left behind, can take part in conduction. First, the free elec- 
trons can move easily from one atom to the next. Second, a bound or associ- 
ated electron on one of the atoms next a positive ion can jump to that ion, 
leaving its own atom ionized. We are thus led precisely to the dual theory of 
conduction, by free and by associated electrons, which Professor Hall* 
has suggested and elaborated. When we look at the metal by the method of 
Heisenberg, these results become clear. In that method, a wave function 
consists of the assignment of electrons to atoms. We find that we must go 
beyond Heisenberg, in assigning sometimes two electrons to one atom, 
sometimes none, instead of always one; for we need such states to solve the 
problem of the stationary states of the metal. That is, we introduce free 
electrons. And when we consider transitions from one state to another, it 
is easy to see that these transitions can result in conduction only when such 
free electrons are present. On Bloch’s scheme, where we describe directly 
the velocity, rather than the position, of the electrons, it is less easy to see 
the relation; but here too one can show that, if there are no free electrons, the 
velocities of all electrons must compensate, so that there is no net current. 
Since this paper is not primarily about conduction, we do not go into these 
points with any detail. 

The only metals specifically treated are those with one valence electron 
per atom, and that in an s state; that is, the alkalies. And it is assumed that 
they can be replaced by single valence electrons moving in non-coulomb 


2 E. H. Hall, Proc. Nat. Acad., 1920-1921. 
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fields. This can be easily justified. It is to be noted that the other metals 
are more complicated, not merely by having more electrons, but by having 
them in ~ or d orbits, thus introducing new degeneracies. The actual cal- 
culations of cohesion have been carried through for sodium, with satisfactory 
results. They are only done roughly, however; the primary purpose of this 
paper is to make clear the general relations, rather than to attempt accurate 
calculations. , The work is being carried further by Dr. Bartlett, and I wish 
to thank him for help on some of the calculations used in this paper. The 
work described here has been done while the writer was on leave, working in 
Leipzig. He wishes to thank Professor Heisenberg for his courtesy in extend- 
ing the privileges of his laboratory, and for a number of illuminating con- 
versations on the subject of the paper; and also to thank Harvard University 
for granting leave, and the Guggenheim Foundation for the assistance of a 
fellowship. 


2. COMPARISON OF HEISENBERG’S AND BLOCH’s METHODS 


The problem of a metal must be attacked by perturbation theory, and 
the unperturbed functions which we use can be set up in two quite different 
ways, one used by Heisenberg, the other by Bloch, either giving us a finite 
set of unperturbed functions. We regard the perturbation problem in the 
following way: we seek those linear combinations of these functions which, 
in the sense of the variation method, form the best approximations to solu- 
tions of Schrédinger’s equation. This problem is solved by computing the 
matrix of the energy operator with respect to these functions, and solving 
the equations 


>> (A(i/k) —8(i/k) W)S(k) =0 


for the coefficients S(k) to be used in making the linear combinations, 
and the energy values W of the resulting terms. (The term 6 (z7/k) must be 
given a slightly different form if the unperturbed functions are not orthogo- 
nal). This differs from the more conventional method: there one starts with 
an infinite, complete set of unperturbed functions, instead of our finite set, 
but solves only as a power series in the non-diagonal terms of the energy 
matrix, breaking off after the second power in all ordinary applications. 
It resembles more closely the quite different method ordinarily used with 
degenerate systems, where one takes only very few unperturbed states, but 
correctly solves the problem of combining them. For a nearly degenerate 
problem like the present one, with a great many states near together, the 
conventional method of developing in series will not work well, {for the series 
do not converge well, and we are forced to use something like the present 
method. The justification comes simply from the assumption that the lowest 
states can be well approximated by such a linear combination of Heisenberg’s 
or Bloch’s functions (which correspond to having the atoms in their normal 
states). Surely this is not exact; for better results we should have to consider 
also the excited states of the atoms. But also certainly it is a fair approxi- 
mation for the lowest states of the metal. 
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Heisenberg’s functions, amplified in a simple way, form good approxima- 
tions when the crystal is extended, for they are derived from the separated 
atoms. Bloch’s functions on the other hand come by analogy with the free 
electron theory of Sommerfeld, and are good approximations when the 
crystal is compressed. The actual solutions of the perturbation problem are of 
course linear combinations of either Heisenberg’s or Bloch’s functions, not 
individual ones, and one gets the same final result whichever set one starts 
with (for the two sets of functions can be written as linear combinations of 
each other). But the fact that in the limiting cases the functions of one of the 
two sets become rather good approximations can be used, along with inter- 
polation, to derive the general nature of the real stationary states. This 
comparison is made in'the present section, and is illustrated by the interesting 
case of He, where the calculations can be made exactly. At the outset, we 
must recognize two facts: first, that we must amplify Heisenberg’s method 
by including polar states, to make it general enough to agree with Bloch’s 
and to permit conductivity; second, that although Bloch has the proper set 
of functions, he has nowhere attempted to solve the perturbation problem, 
but has merely taken his unperturbed functions as being correct, which 
amounts to getting the energy to the accuracy of the conventional “first 
order perturbations.” 

The first step in either Heisenberg’s or Bloch’s method, as we apply them, 
is to write an approximate solution as a product of functions of the individual 
electrons. Heisenberg takes, for these separate functions, the wave functions 
of electrons attached to individual nuclei; the number of such functions is 
the product of the number of nuclei, multiplied by the number of different 
sets of quantum numbers we consider for an individual nucleus. If we 
restrict ourselves to s states, there are then only two states per nucleus, 
corresponding to the two orientations of the spin. For nucleus a, we denote 
these two® by u.(a/x), us(a/x), and we have such a function for each nucleus 


ab...n. Bloch takes, on the other hand, combinations of these functions: 
Ua(klm/ xyz) = > e2tikorlGrttorlGot moslGa) 4, (g gog3/ x2), 
919293 


where gigeg3; are the coordinates of a particular nucleus, G,;G.G; the dimen- 
sions of the rectangular crystal, and wa(gigeg;/xys) the wave function (as 
used by Heisenberg) for an electron moving around the nucleus situated at 
gigeg3. The function with k, ], m represents an electron, in general moving 
in the direction k, 1], m, but pausing at the various atoms on the way. There 
are as many sets k, /, m allowed as there are atoms in the crystal; for larger 
k, 1, m the function proves to be merely a repetition of one already counted. 


3 We use here for convenience in writing Pauli’s notation uw, ug for the spin, rather than 
the more explicit but more cumbersome notation u(n/x;) 5(m;/m.) used in a previous paper. 
See J. C. Slater, Phys. Rev. 34, 1293 (1929). The method used in the present paper is described, 
as applied to atoms, in the paper referred to; it should be understood that, although we speak 
here of using Heisenberg’s and Bloch’s methods, our actual procedure is quite different from 
that of these authors. 
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Now we actually set up the product of functions mentioned in the 
previous paragraph: we pick one out and let it be a function of the coordin- 
ates x, of the first electron, a second for the coordinates x2 of the second, and 
so on to the mth, and multiply them all together. By the exclusion principle, 
no function can be chosen more than once. Then we form an antisymmetric 
combination, by permuting the indices of the electron coordinates, and add- 
ing the permuted functions with appropriate signs, obtaining essentially a 
determinant. These antisymmetric functions are the ones with which we 
start our perturbation calculation. Many such functions can be set up: 
there are 2 functions of a single electron, of which only m are to be chosen 
for each antisymmetric function, so that there are (2m)!/(m!)? different 
functions. Our perturbation problem is that of finding which linear combin- 
ations of these functions most nearly satisfy the wave equation. We may note 
the restriction of Heisenberg’s method as he uses it; he does not include polar 
states. That is, he does not allow for example the two functions u,(a/x), 
us(a/x) to appear together in any product. This greatly limits the number of 
functions; but although the terms obtained by it certainly represent the 
lowest energy levels, since it requires energy to form a positive and a nega- 
tive ion from two neutral atoms, we do not make this limitation. 

Having set up the unperturbed functions, we next make linear combin- 
ations of them, by the method described in a previous paragraph. This 
process can be simplified by using a property of the spin. Every unperturbed 
function has a certain definite component V/s of spin along the axis, equal to 
(t%a—ng)/2, where n, is the number of electrons with positive component of 
spin, ms the number with negative. If now we neglect the magnetic inter- 
action between the spins and the orbital motion, the problems with each 
value of Js can be handled separately: the components J//(i/j) from a 
function with one value to a function with another are zero. The states with 
a given Js include, as one readily sees, all those states whose total spin S 
is equal to or greater than M/s (for just these S’s can be so oriented, on the 
vector model, as to give a component Ms along a fixed axis). Thus by solving 
each such problem, and comparing, we can identify the spin of each term.‘ 

The two methods can be illustrated by the case of Hz. Here there are 
4!/(2!)?=6 different wave functions. On Heisenberg’s method, the four 
functions for an individual electron can be symbolized by (aa), (8a), (ab), 
(8b) ; two of these are to be picked out for each antisymmetric wave function. 
Thus the six are (aa)(ab); (aa)(Ba), (ab)(8b); (aa)(8b), (Ba)(ab); (Ba) (8d). 
They are arranged, first, by Ms: the first has the value 1, the next four the 
value 0, and the last —1. Thus the terms consist of one triplet and three 
singlets. Among the four terms with 1/s=0, the first two are polar (and not 
considered by Heitler and London, or Heisenberg), the last two are non- 
polar. Immediately one finds that the sum of these non-polar functions is 
the component of the triplet. We are then left with three functions: the 


‘ This is essentially the method used in the paper already. quoted It has already been 
applied by Bloch to problems in the theory of metals. See F. Bloch, Zeits.f. Physik., 57, 545 
(1929), 
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two polar ones, and the difference of the non-polar ones, from which to find 
our three singlets. The difference of the polar ones is antisymmetric in the 
nuclei, giving one state; their sum, and the difference of the non-polar func- 
tions, give two functions symmetrical in the nuclei, between which we finally 
solve the simple perturbation, resulting now in a quadratic secular equation, 
and obtain the two remaining singlet states. The energy levels as a function 
of the distance of separation are plotted in Fig. 1. The energy level of the 
lowest, 'XS*, is almost exactly as given by Heitler and London, but its wave 
function contains quite an appreciable contribution from the polar state. 
The triplet is just the repulsive state of Heitler and London. The other two 
levels are essentially polar. They go at infinite separation to the energy of 
H*++H_-, greater than the other limit by the ionization potential less the 
electron affinity of H (this rough approximation gives — } Rh for the electron 
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Fig. 1. Energy levels of He. 


affinity, so that the terms go to 5/4Rh). The lower of these has a minimum; 
it is presumably the polar part which, by combination with other functions, 
leads to the experimentally known B state of the molecule. We notice that 
at large separations the functions behave just like Heisenberg’s (extended) 
unperturbed functions: a triplet and a singlet are non-polar, and go to the 
lower energy; while two singlets are polar, and go to the higher level. 

Next we consider Bloch’s method for the same problem. His functions 
for one electron, for this case, are 


Ua(O/ x) = Ua(a/ x) + a(b/ x) 

Ua(1/x) = Ua(a/x) —Ua(b/x), 
with similar functions for 8. (These do not follow quite directly from the 
general formulas given above; Bloch’s functions must be slightly modified 
for finite systems, for they apply rather to infinite but periodic ones.) 


The discussion of multiplicity given above goes through without change, if 
we only substitute 0, 1 for a, b. We can easily show by direct calculation 
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that the resulting unperturbed antisymmetric functions are linear combin- 
ations of those found by Heisenberg’s method. For example, for M/s =1, 
there is only one function by either method, so these must be identical, 
except for a numerical factor. By Heisenberg’s method the function is 


Ma(a/ %1)Ua(b/ x2) — ta(b/ x1) Ma(a/ x2). 
By Bloch’s it is 
Ua (O/ X31) Ma(1/ x2) — ta(1/ 21) tta(0/ x2) 
= [ta(a/x1)+tta(b/ x1) | [tta(a/ x2) — Uq(b/ x2) | 
— [ta(a/ x1) — ta(b/ x1) | [ta(a/ x2) + ta(b/ x2) | 
= —2[ta(a/21)ta(b/ x2) —tta(b/ %1)q(a/ x2) J. 


We can set up the whole perturbation problem in these functions; and the 
solution can be carried out as easily as before, leading of course to just the 
same answers. The interesting question now is, how closely do Bloch's 
individual functions approximate the correct ones, for small values of R? 
The functions are respectively as follows: a singlet with both electrons in the 
state 0; a singlet and triplet with one in the state 0, the other in the state 1; 
and a singlet with both in the state 1. The state 0 corresponds to the lowest 
vibrational state on Sommerfeld’s theory, the state 1 to the next higher one, 
so that the first state has on the simple interpretation only the zero-point 
vibrational energy, the next two have each one quantum, and the last 
two. Examination of the actual wave functions shows that they agree quite 
closely with the functions of Bloch: the lowest one is made, it is true, by 
combination of the (00) and (11) states, both being S*, but the coefficient of 
the first is about eight times as large as that of the second, when R is such that 
the energy is at its minimum. The next two are made up of the (01) states. 
The highest is about eight parts of (11) to one of (00). The energies also show, 
for high compression, the behavior expected: the two states which should 
have one quantum of vibrational energy draw together, and the one with two 
quanta is just about twice as far above the lowest state as those with one. 
Even the spacing of these levels is just about what would be calculated on 
Sommerfeld’s theory for an electron vibrating in a region the size of the mole- 
cule. Thus we see that Bloch’s unperturbed functions form fairly good ap- 
proximations to the real functions for the compressed state, as Heisenberg's 
do for the extended state. 

We can now return to the general case, and make use of the fact that 
Heisenberg’s functions approximate the real wave functions well for large 
separations, Bloch’s for small. First, for the extended system, the energy 
is the ionization energy, on account of having many ions as well as neutral 
atoms. For a metal, it requires about 6 volts to form a positive and negative 
ion from two neutral atoms. Thus if all the atoms were ionized, we should 
have /2 such pairs, or an energy per atom, or per electron, of about 3 volts. 
This measures the extension of the group of terms, for large R. It is a simple 
problem in permutations to find the number of terms of each multiplicity 
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with each energy value in the limit. The one term of highest multiplicity 
will approach the lowest limit, for large R; it must be entirely nonpolar, for 
all spins point in the same direction, so that no two electrons can be in the 
same atom. For the next lower multiplicity, only one electron has a reversed 
spin; it is the only one which can be in an atom with another electron, so 
that there can be just one pair of atoms ionized. Following out, we easily 
see that terms of lower and lower multiplicity, in the limit of large R, lie 
higher and higher, and at the same time are more and more spread out. They 
spread in such a way that there are terms of each multiplicity way down to 
the bottom limit, although not to the top. As we shall see later, for R large 
but not infinite, in the normal case, the really lowest terms have small spins; 
but near them are many terms with large spin. 

_ For the compressed system, the arrangement is as given by Bloch’s 
theory. The total extension of the group of terms increases with 1/R?*; for 
ordinary values of R, it is of the order of the mean zero-point energy, times 
n, which is decidedly larger than 3 volts. Thus not only do the curves 
tend upward for decreasing R, giving repulsive energy levels, but they are 
definitely doing this at the actual size of the metal. The general physical 
interpretation of this repulsion is obvious: the valence electrons act here 
approximately as a perfect gas, and the energy levels are those of such a gas 
as it is compressed adiabatically against gas pressure, the energy varying 
therefore as V~** or as 1/R*. Here the terms of high multiplicity lie in the 
center of the pattern; those of lower spin also average in the center, but are 
more and more spread out. Since the terms of high spin are so low for large 
R, but not for small R, they must be even more repulsive than the others. The 
possibility seems very remote that any terms except those with very low 
multiplicity could be so low as to have minima, and come into the question 
for the normal state. We see that for cohesion we are interested only in the 
very lowest fraction of the whole set of terms. These terms almost all will 
go to the lowest energy level at infinite separation; they become in this 
limit non-polar. And the accuracy with which one can compute the lowest 
states of H. from Heitler and London’s non-polar functions suggests that 
here too this may be possible. Accordingly for our actual calculation of these 
lowest states, we shall use Heisenberg’s method with only non-polar func- 
tions. We shall find here, as we expect from our qualitative discussion, that 
the terms of low multiplicity really do lie below, some of them being attrac- 
tive; while those of high multiplicity are repulsive, the highest spins lying 
highest. Finally we shall consider the effect of polar terms, and conclude 
that it is really small on the low energy levels, although not on the wave 
functions; for it is the polar character of the wave functions which makes 
conductivity possible. 


3. ELECTRIC AND MAGNETIC PROPERTIES 


Conductivity. In the introduction we have mentioned the interpretation 
of electric conduction on Heisenberg’s and on Bloch’s scheme. One notices 
that a single one of Bloch’s functions implies conduction—the diagonal term 
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of the momentum matrix is different from zero—whereas with Heisenberg’s 
functions we must have a continual change from one stationary state to 
another. But it is particularly important to notice that, without polar 
states, or free electrons, no conduction is possible; we cannot set up combina- 
tions of non-polar states with a resultant momentum. For example, with two 
electrons, we can set up an arbitrary non-polar function cu(a/x)u(b/x2) 
+c.u(b/x,)u(a/x2). If now we compute the momentum, whose operator is 


h/2mi(0/0x,+0/0x2), the only possibly significant terms are the cross terms, 
like 


~f , r) rs) b/s ded 
ts u(a/ x) u( /9(— +) / x1) U(a/ X2)dv,dv2 


2 


Tt 


=o} f u(8/x:)u(a/xa)doe f u(o/s)-—u(b/ ade 


ad f u(a/s)u(8/xi)dor f u(8/2)—ula/s)dos} 


On account of the penetration of one atom by the other, the integrals 
fu(b/x2)u(a/x2)dve are not zero. The integral fu(a/x1:)(0/0x:)u(b/x1)d2 is also 
different from zero. But it is exactly cancelled by [u(b/x2)(0/0x2)u(a/x2)dvs, 
as one can show by Green’s theorem, so that the whole is zero. On the other 
hand, if we set up a polar combination like cyu(a/x;)u(a/xe) +c2u(b/x1)u(b/x2), 
we again get two terms, but now they add, and give a current. As another 
example, we can take the term of maximum multiplicity in any system. In 
this term, we have seen by Heisenberg’s scheme that each atom has just 
one electron, so that we expect no conduction. But in Bloch’s scheme, each 
value k, /, m has just one electron. Since each such value is balanced by one 
with —k, —/, —m, having opposite momentum, the total momentum is zero, 
and there is again no current. 

We can now see the importance of considering exactly the wave functions, 
as well as the energy levels, of the lowest state. In the ordinary low states 
there will, of course, be no current. But near the lowest state, if there is to 
be conductivity, there must be combinations of polar states, having a cur- 
rent, which are assumed in the presence of a field, and whose added energy 
comes simply from the kinetic energy of the electrons and the self-induction. 
Such states are possible only on account of the presence of positive and 
negative ions, with the resulting free and associated electron conductivity. 

Magnetism. The lowest state of He is the non-magnetic '2, and we have 
found such a situation in general. In the region where the lowest states have 
their minimum, the metal must surely be in a compressed state, Bloch’s 
arrangement of energy levels must be a good approximation, and the states 
of large spin must lie very high. We are thus led to the quite general conclu- 
sion that the outer electrons, which are largely if not entirely responsible 
for both cohesion and conduction, cannot produce ferromagnetic effects. If 
a metal is to be ferromagnetic, there must then be other electrons than these 
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outer ones which are responsible for it, and these others must have smaller 
orbits, so that at the equilibrium distance of the outermost ones, the inner 
ones will be relatively further apart, and can be treated as an extended 
rather than as a compressed system. It is a very attractive hypothesis to 
suppose that in the iron group the existence of the 3d and 4s electrons provides 
in this way the two electron groups apparently necessary for ferromagnetism; 
for it is only in the transition groups that we have two such sets of electrons, 
and this criterion would go far toward limiting ferromagnetism to the metals 
actually showing it. 

We next ask just how such inner electrons could be ferromagnetic. Cer- 
tainly the general trend of the terms of high spin to the low energy values 
at large R is an essential part of the question: there will be terms of large 
spin near the lowest level. Bloch® has discussed the problem, concluding 
that for large R’s the terms of high spin actually lie lower than those of 
smaller spin (he does not specifically discuss the dependence on R, but his 
energy formulas all contain it parametrically). This conclusion, however, is 
not correct; Bloch has merely computed diagonal values of the energy, with 
respect to his functions, and for large R these by no means form approxima- 
tions to the actual energy values. From the correct treatment of the problem 
as we have given it, it is plain that at all R’s there are terms of low multiplicity 
as low as those of high spin, or lower. It may be, however, that the mere 
presence of so many low terms of high multiplicity may be enough, on ac- 
count of their high @ priori probability and large number, to insure that the 
terms of large spin should be well represented at ordinary temperatures, 
even though there are low terms of zero spin, and so produce ferromagnetism. 
If, however, this should prove on calculation not to give the right effect, we 
should be led to consider Heisenberg’s assumption that the normal order of 
terms is inverted in ferromagnetic atoms, the terms of high multiplicity lying 
lowest. He has shown by a general argument that electrons of large total 
quantum number (which the 4s electrons of iron have) have an exchange in- 
tegral of the opposite sign to that found in hydrogen, so that the order of the 
non-polar terms would be reversed. This we should fit into our scheme in the 
following way: although this exchange integral is anomalous at large R, it 
presumably changes sign and becomes normal at smaller R; for first, Heisen- 
berg’s general argument only applies at large R; and second, our condition 
that the energy levels should approach those of Bloch at small R, with the 
terms of large spin lying high, seems quite general. Thus we should assume 
that the terms at small R lie as in Fig .1: but that at a considerable value of 
R, there is a crossing over (in this case the *SA¥ crossing and lying under the 
'¥S*), described by a change in the sign of the exchange integral K used 
in the next section from negative to positive. By assuming the existence of 
an inner group of electrons with these properties, we seem to secure a con- 
sistent picture of ferromagnetism. On the other hand, of course it is always 
possible that ferromagnetism is connected with the fact that the valence 
electrons of iron have an orbital angular momentum different from zero. 


5 F, Bloch, Zeits. f. Physik, 57, 545 (1929). 





a ct 














a ee Nee 











COHESION IN MONOVALENT METALS 519 


4. COHESION 


We are now prepared to begin the actual calculation of the lowest sta- 
tionary states. We make several simplifications, which we later remove. 
First, we consider only non-polar states, in Heisenberg’s scheme, and disre- 
gard exchange integrals except between adjacent atoms; this is the approxi- 
mation also made by Heisenberg. Finally, for the present, we consider a linear 
lattice, 2 atoms uniformly spaced along a line, rather than a space lattice. 
Our problem, of course, is to compute the matrix of the energy with respect 
to the wave functions we have chosen, and then solve the problem of making 
proper linear combinations. The computation of the matrix is simple. By a 
fundamental formula of the previous paper mentioned above, the diagonal 
components are a sum, first, of the energies of the separate atoms, which we 
need not consider; next, a sum over all adjacent pairs, as the pair of atoms 
a and 3b, of integrals J(a/b), which is essentially the diagonal energy F, of 
Heitler and London; finally, a sum over all adjacent pairs which have the 
same spin, of terms—K(a/b), where K is the exchange integral £2 of Heitler 
and London. Further, it is easy to show that all non-diagonal terms are 
zero, except those for which the distributions in the two states differ only by 
the exchange of an adjacent @ and §; in such cases, the term is —K. In the 
normal case, to which we shall refer specifically, J and K as functions of R 
are both negative, K numerically greater than J. But in Heisenberg’s case, 
K must be taken to be positive for large R, although presumably negative 
for small R. 

To illustrate by Hz, we have one state with both spins parallel; then the 
energy is J/—K. Next we have the problems with one parallel, the other 
anti-parallel; there are two such states (the two polar ones being omitted). 
Each has the diagonal energy J, and the non-diagonal energy between them 
is —K. Thus the equations for the linear combinations are 


(J —W)S(1) —KS(2) =0 
— KS(1)+(J—W)S(2)=0, 


giving energy values W=J+K, the first evidently being the singlet, the 
second the component of the triplet. 

In the general case, the computation of the matrix is no more difficult; 
the real problem is the solution of the linear equations for the S’s. We cannot 
do this exactly; but we adopt two methods of approximation, one holding 
for larger spins, the other for smaller spins. We first discuss the former. 

Method for large spins. First we take the problem where all spins are 
parallel, n,=n, ns=0. Here there is but one state. Since with our linear 
lattice there are (n—1) adjacent pairs, and all spins are parallel, the energy 
is simply (n—1)J—(n—1)K. Since J and K are normally both negative, 
but K numerically greater than J, this is a positive energy for all values of R, 
and results in a repulsive term. For Heisenberg’s case, on the other hand, K 
is positive, and this term is attractive. Next we take the problem n,=n”—1, 
ng=1. There are now » unperturbed wave functions: the one electron 8 can 
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be attached to any of the ” atoms. We number the functions by the number 
of the atom where the electron is, only decreased by }: we have t1/2 - - - Up_1/2. 
Each of these functions will have the diagonal energy (n—1)J—(n—3)K, 
since two of the adjacent pairs now have opposite spins, except for the two 
functions 1/2 and “,_1/2 where our 8 electron is at an end of the lattice, and 
the energy is (n—1)/—(n—2)K. Also, all non-diagonal terms will be zero 
except those between terms of adjacent number, as for example between 
those symbolized by 


-aaapaaa:: 
and 
-aaaaBpaa:::, 


and which differ by just one interchange of anaand8. Asaresult, the pertur- 
bation equations will be 


[(m—1)J —(n—2)K— wis(—)- xs(=) =0 


-Ks(—)+ [(m—1)J —(n—3)K— whs(>)- xs(=) =0 


3 : 1 
-K3(n->)+[(n- 1)J —(n—2)K— W]s(n-—)=0 


These equations are easily solved; they occur, for example, in the problem of 
a string weighted at equal intervals,® the S’s being the displacements of the 
weights. To solve, we merely assume S(k)={i (ak). The first and last 
equations give boundary conditions. They become like the others if we 
introduce an S(—}3) and S(n+3), the first equation becoming 


—KS(—3)+ [(n—1)J—(n—3)K-—W]SQ)—KS(3)=0, 


and if we further set S(— 3) = S(3) and S(n+3) =S(n—}3). These are then the 
boundary conditions; and to satisfy them we must take 


S(k/p)=cos prk/n, where p=0,1,--- , m—1. 


Now we substitute this form in our difference equations; and we get for W 
pr pr 2 L prk 
=K(cos=(4—1) +005 (b+1))+ [(m—1)J —(n—3)K—W(p)] cos ——=0, 
n n n 
from which in each case 
T 
W(p) =(n—1)J —(n—3)K—2K cos Pn ‘ 
n 


6 See, for example, Rayleigh’s “Theory of Sound.” 
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We now have the transformation coefficients S(k/p) and the energy values 
W(p) for the rotation of axes to the pth stationary state; the W’s are the 
exact energy levels. They are evidently distributed between the values 
W(0) = (n—1)J—(n—1)K,and W(n—1)=(n—1)J—(n—3—2 cos r(n—1)/n) 
K=(n—1)J—(n—5)K, almost, for large n. Obviously W(0) is the energy 
of the level of highest multiplicity, which we have found before. Thus the 
levels 1 - - - (7—1) are those of next to highest multiplicity. 

Next we take the problem with two electrons of spin 8. There are 
n(n—1)/2 such terms: each of the two indistinguishable #’s can be on any 
one of the 7 atoms, so long as they are not on the same atom. Now it is con- 
venient to denote states by the two atoms, say 7 and j (each going from } to 
n—%) on which electrons 8 are. Our problem becomes analogous to that of a 
square membrane loaded at equally spaced points. The diagonal terms 
of the energy are all (n—1)J—(n—5)K, unless one of the @’s is at an end of 
the lattice, or unless the two §’s are adjacent. There are four non-diagonal 
terms for transitions from each wave function: for z~7+1, or for j->j+1. 
A typical equation can be written 


— KS(i,j—1) 
—KS(i-—1,j)+[(n—1)J —(n—5)K—W ]S(ij) —KS(i+1,j) 
— KS(i,j+1)=0. 


This we satisfy by a product of cosine functions, S(ij/pq) =cos (pai/n) cos 
(qmj/n). We easily find that these exactly satisfy the boundary conditions 
when i or j=} or n—}. There remains the condition when 7 is nearly equal 
to j. If <=j+1, the diagonal energy is (n—1)J—(n—3)K, since the two f’s 
are together; on the other hand, since the 6’s cannot be on the same atom, the 
coefficients S(jj) and S(j+1, j7+1) vanish, so that only two transitions, 
rather than four, are possible. If now we define an S(jj) and S(j+1, 7+1), 
we can make the equations of the same form as the general one, if only 
Sj) +SG+1, 7+1)=2SG+1, 7). This furnishes our second boundary 
condition, which is evidently along the diagonal of our square “membrane.” 
Unfortunately we cannot satisfy this condition exactly with our cosine func- 
tions; closer investigation shows that one must have much more complicated 
functions, with hyperbolic cosines, to satisfy it exactly, and one cannot carry 
the method through for the general case. Approximately, however, we can 
easily take care of our condition. If the and q are not too great, so that the 
“wave-length” of the waves in our membrane is large, we can replace our 
condition by a differential one: it states that the amplitude at a point next 
the diagonal is the mean of the two adjacent values on the diagonal, and this 
very nearly means that the normal derivative of the function, at right angles 
to the diagonal, is zero. This we can satisfy by making our function symme- 
trical about the diagonal, or using cos (pri/n) cos (qmrj/n)+cos (qri/n) cos 
(pmj/n). We may expect this to hold best for small » and g, not so well for 
large values. It is clearly not right; for example, it yields n?/2 functions, in- 
stead of the correct number n(n —1)/2. 
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Our function is an exact solution of the difference equations, if not of the 
boundary condition; and we find for the energy 


—_ : pr qm 
W( p,q) =(n—1)J—(n—S5)A =2K (cor +eos), 
n n 
where ~p, g go from 0 to n—1, but each pair is counted only once. The term 
of highest multiplicity comes from p=q=0; the (7 —1) terms of next highest 
value are those with either p or g=0, but the other not; the remaining terms 
are of multiplicity smaller by two. 
This result can now be generalized without trouble: if we have many 
B’s, the energy levels are given by 





n3 pir 
W=(n—1)J —(n—1—2n;)K-—2K p>» cos » pi=0---n—-1. 
m1 n 


The terms where one or more /,’s equal zero are those whose total spin is 
greater than (7.—m3)/2; those with all p’s different from zero are those whose 
total spin equals (m_—3)/2. The latter value is evidently enormously 
greater than the other: every spin has enormously more terms than any 
higher spin. Thus the terms of a given component of spin along the axis, and 
those of the same total spin, are approximately the same. We can at once 
find the distribution in energy of the terms of a given spin. They evidently 
cluster about the value (n—1)/J—(n—1—2n3)K; they are distributed about 
this value like the displacements of a point simultaneously acted on by a sum 
of m3 periodic vibrations of equal amplitudes but arbitrary phases. This gives, 
of course, approximately a Gauss distribution. The width of the distribution 
curve can be derived very easily: we compute the mean square deviation 
of the energy from its mean, (|j'—]| )?=4K"%[S cos (p;z/n)}?, the average 
being taken when each # varies independently from 0 to m. We can take 
this variation to be continuous rather than discrete. Then the product terms 
in the square of the sum of cosines all average to zero, the square terms 
average to 3, and the result is 2K°n3. These results may be compared with 
those obtained by Heisenberg on the group theory, and which as Bloch has 
shown can also be found from the present method. In the notation of the 
present paper, putting the number of neighbors of each atom equal to 2, and 
leaving out the terms in J, Heisenberg finds 


W = —(n—2n,+2n,2/n)K 
(Ww W)?=2K2n,(1—n3/n)(1+2n,/n—2(n;/n)?). 


Our formulas agree with these exact ones to terms in m3 but no further, as we 
expect from the fact that our approximations hold only for small as. For 
small p’s, as we have seen, our results should be good even for large m3; for 
the case of ferromagnetism, when on Heisenberg’s hypothesis the terms are 
reversed, these are the lowest terms, so that this result should be very useful 
here. In the normal case, however, the lowest terms are those of large p, 
and these are the ones we need for cohesion, About these lowest terms, we 
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can be fairly sure that they lie higher than the lowest ones we have found, 
or (n—1)J—(n—1—4ng)K, since the mean lies higher than the mean we 
found. For zero spin, for example, we can be fairly sure that the term lies 
above (n—1)J+nK. But this value need not be a very good approximation; 
we actually find, by the method of the next section, that the lowest term for 
zero spin is about (n—1)J+0.290”nK. Fortunately even this has a positive 
coefficient for K, and is so an attractive rather than a repulsive term. 

Method for zero spin. For zero spin, n,=n3=n/2, and there are n!/(n/2!)? 
terms. We adopt quite a different method of classifying them. Before, most 
of the terms of a given mg had nearly the same diagonal energy; but now the 
range of energy is large, from (m—1)J for the state with alternating a’s and 
8’s so that there are no parallel spins, to (n—1)/—(n—2)K for the state 
where all the a’s come at one end of the lattice, all the 8’s at the other. With 
this large range, we find it convenient to classify terms by their diagonal 
energies; and we find as we should expect, that for the lowest states of the 
perturbed system we must consider most the low unperturbed states. We do 
not need to take into account all states: we find that approximately (though 
by no means exactly) the terms can be divided into a number of non- 
combining sets, and we set up one such set in the following way. We com- 
mence with the lowest state, of energy (n—1)J, where a’s and @’s alternate. 
Next we consider the —1 states which combine with it, coming from inter- 
changing one pair, and each having the energy (~—1)/—2K, except those 
from the two end pairs, with energy (n—1)J/—K. We leave these two out, 
retaining for our set the »—3 states which have the energy (n—1)J—2K. 
Each of these has »—3 states with which it combines, coming from inter- 
change of one of the n —3 adjacent pairs with opposite spins. Of these — 3, 
the two in which the new interchanged pair is next the one already inter- 
changed have the energy (w—1)J—2K; the one in which the pair already 
interchanged is changed back has the energy (~—1)J; the two where the 
end pair is interchanged have the energy (n—1)J—3K; and the remaining 
n—8 have the energy (n—1)J—4K. We retain for our set only these n—8& 
terms of energy (7—1)J—4K. So we proceed, asking which terms combine 
with those already set up, and retaining just those whose energy is —2K 
greater than for those with fewer interchanges. We find that a term of our set, 
with the energy (7»—1)J—2pK, has non-diagonal terms to p terms of the 
set of energy (n—1)J—2(p—1)K and to (n—3—5p) terms of the set of 
energy (n—1)J—2(p+1)K. Evidently so long as p is small, the terms we 
leave out of the set and yet which combine with terms of the set are com- 
paratively few. It is only for the large p’s that we make serious error by 
leaving out these terms, and for large p the diagonal energy is high enough 
so that for the lowest states of the perturbed system these unperturbed 
states are unimportant. Thus we may reasonably believe that the low 
energy levels found by solving this restricted problem will be approximately 
some of the low levels of the actual problem. We can at least be sure of the 
following: by the variation principle, they can be no lower than the actual 
stationary states. 








524 J. C. SLATER 


The other sets of non-combining terms which we can set up are easily 
described, and are of considerable physical interest. Instead of starting 
from the state with spins alternating, we start from a state where the spins 
alternate up to a given point; there the sequence is interrupted, and alter- 
nation commences again, so to speak, in the opposite phase, as 


--aBaBBaBaBa:::. 


With a few such interruptions in the course of the crystal, the energy is very 
little above the really lowest state; yet a great many individual inter- 
changes would be required to pass to the lowest state. With such a state to 
start with, we proceed just as we did before, and construct a whole system 
of states; and the non-diagonal terms between this and the first system come 
only from high values of /, involving many interchanges, and can be neg- 
lected. Physically, at the interruption of phase, one essentially has.a slight 
interruption of crystal structure. Our catalogue of all possible states of the 
metal includes not only that where it is one perfect crystal, but also where 
it is composed of many smaller crystals not perfectly joined together. Ob- 
viously each problem can be treated separately; physically it would take a 
very long time to change from one to the other. And obviously each problem 
will give us essentially the same set of energy levels. 

We now take our set of wave functions, and try to solve the perturbation 
problem between them. For each value of », we have many wave functions; 
and we look for those particular solutions for which all these functions have 
the same coefficient S(p). Afterwards we shall show that we really find the 
lowest solutions this way. Then, remembering the number of transitions with 
non-diagonal term K from a given state, computed above, we have for a 
typical equation ; 


— KpS(p—1)+[(n—1)J —2pK —W |S(p) —K(n—3—5p)S(p+1) =0. 


This set of difference equations for the S’s is somewhat similar to what we 
had before; it also corresponds to a weighted string. But now the properties, 
and hence the wave-length, change from point to point, and we seek the vari- 
ous overtones. The equation is a close analogue to Schrédinger’s equation, 
in many ways; the fact that it is a difference equation rather than a differen- 
tial one is quite immaterial. To solve, we assume S(p) =e/*?, where a is to 
vary slowly with p. Then S(p) =e*S(p—1), etc., so that we have 


—Kp+e*|(n—1)J —2pK —W]—eK [n—3—5p]=0, 


_ —[(@-1)J-2pK-—W] + ([(n— 1) J—2pK — W}*—4K2p[n—3—-Sp])"? 
- —2K|n—3-—5p] | 





The equation expresses e* as a function of p, for any particular W. Now we 
must remember that there are essentially boundary conditions; the S’s must 
remain finite for =0 and p=an extreme value. To tell how to apply this — 
condition, we must investigate the solution we have found. 
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The ratio e* of successive coefficients is real or complex, according as 
|(n —1) J —2pK — W]? is greater than or less than 4K2p[n—3—5p]. Regarded 
as a function of p, the limiting cases, where the two are equal, come from 
(n—1)J—W=2pK+2K|p(n—3—5p|"*. The right hand side, plotted as a 
function of ~, forms an ellipse; the straight line represented by the left side 
intersects the ellipse in two points, or in none, depending on the value of W. 
The region of W where it intersects can be found by computing the maximum 
and minimum ordinates of the ellipse; that is, the values of (n—1)J/—W 
for which (d/dp)(2pK +2K|[p(n—3—5p)]*)=0. This gives p=(n—3)/10 
[1+ (1/6)"?] =(n—3) X (0.0592, 0.1408). At these two limits, substituting, 
W=(n—1)J+(n—3)K X(0.290, —0.690). For values of W between these 
limits, there is a range of p for which e* is complex, and the solution is 
oscillatory; outside this region, which is closed, the solution is in any case 
exponential. To satisfy our boundary conditions, now, we have a problem 
much like that with Schrédinger’s equation in one dimension; and boundary 
conditions can be satisfied only if there is an oscillatory region. As a result, 
the actual energy levels of the problem must lie between the limits given. 
Closer examination shows that a “quantum condition” can be applied, and 
that between these limits there are just the number of energy levels there 

should be. We now have the lowest level: it lies arbitrarily close to our lower 
limit, or is 


W =(n—1)J+0.290(n—3)K, 


as we stated in the last section. In this lowest state, we can show without 
trouble that the unperturbed wave functions with p near 0.0592(m—3) are 
represented most strongly. Thus the value of p is really quite small; relatively 
few pairs are interchanged, and we are safely in the region where we can treat 
the different systems separately. 

We have solved our problem for the lowest state in which all terms of 
the same p have the same coefficient. We can now investigate the effect of 
removing this assumption, varying the coefficient of one function of a given 
p in one direction, varying the rest to keep the same total representation for 
functions of this », and calculating the change in the energy. When we do 
this, we find the energy to be a minimum with respect to such variation; in 
fact, the changes of energy compensate each other to a higher order, showing 
that the problem is nearly degenerate with respect to these coefficients. 
Thus we may be rather confident that we have a good approximation to the 
lowest non-polar states. It is of course obvious that this method becomes 
worse as we go to higher states. 

It is instructive to ask what ordinary perturbation theory would give us 
for the lowest state. The lowest unperturbed state has the diagonal energy 
(n—1)J; this represents the ordinary first order perturbation calculation. 
Now we pass to the second order calculation. The lowest state is not de- 
generate, so that we can use the power series development method. The next 
term in the expansion is 2,(/7;;H;;)/Ei—E,;), summed over all excited states 
j. Now there are non-diagonal terms 4); only to the n—1 states with p=1. 
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Thus the //,;’s are all equal to — A, and the energy differences are all given 
by £;=£,—2K. Thus we have as second approximation 


K? 1 
We=(n—-1)J+(n—- * (n— +e 1)K. 


This differs from our result in having the factor } rather than 0.290; we 
have only the first term of a series development, but it is reassuring that 
agreement is as good as it is. For finding the order of magnitude, we could 
use this term alone; we shall find this simple method useful with the space 
lattice. 

Effect of polar states. One can make an estimate, by a method like that 
used here, of the effect of the polar states in depressing the non-polar ones, 
which alone we have considered so far. We can build up a series of states by 
starting with a given non-polar state; then removing one of the electrons to 
an adjacent atom, producing a positive and a negative ion; then removing a 
second; and so on. The series of states so found behave formally like those 
used above. If we solve the problem by the. previous method, or by the 
second order perturbation method, we get a further depression of the lowest 
state, which again can be written as 


(n—1)X >> 


square of non-diagonal term 





energy difference 


The non-diagonal term which comes in here is presumably of the same order 
of magnitude as before, although it is a somewhat different integral. But 
the diagonal energy difference is now essentially an ionization energy, which 
is of the order of several volts, rather than the fraction of a volt that K is. 
Thus the effect on the energy is a number of times smaller than what we 
found before, and we can neglect it. It is not worth while calculating more 
accurately, in this approximation; for with He, it appears that on account of 
the lack of orthogonality of the wave functions, the actual depression of the 
energy is very much less than this rough method would indicate, although 
the effect on the wave function is about what we should expect. One can 
reasonably believe for this reason that the polar states in the crystal depress 
the energy only very little. But we recall that their effect on the wave func- 
tion is to introduce free electrons. By our rough method described above, 
we infer that the fraction of free electrons is of the order of 1 percent, for 
reasonable choice of the constants. This could easily be in error by a factor 
of 10 either way; but at least we see that a definite meaning can be attached 
to the number of free electrons, and that there is a definite procedure for 
calculating this number. 

Normalization and orthogonality. We have not considered the lack of 
orthogonality of the wave functions, resulting in factors like the 1/1+S of 
Heitler and London. When one tries to do this, one immediately strikes a 
difficulty which appears insurmountable: the factor in the denominator, 
instead of being like 1+, is like 1+”S+ ---, where m is the number of 
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atoms, so that the term ”S is enormous compared with unity. On examina- 
tion of simple cases, it appears that the remaining terms, coming from other 
permutations than the simple interchange, are also important, in many cases 
the terms almost cancelling each other. Further, in every expression for 
energy, like the simple (J + K)/(1+.S), there are more terms in the numera- 
tor, also of great importance. But the simple cases give no suggestion of how 
to treat the general case. The key to this difficulty comes from Bloch’s 
method. For example, the term of maximum multiplicity has one function, 
which can be expressed either by Heisenberg’s or Bloch’s functions. 
But the difference is that Bloch’s functions are really orthogonal, unlike 
Heisenberg’s, so that we meet no such difficulty. Of course, the same terms 
occur, but now in the normalization of the individual functions. And the 
numerators, and denominators like 1+nS+ ---, appear as products of n 
factors, each of approximately a simple Heitler and London form; further, 
all but one or two of these factors of the denominator cancel against equal 
factors in the numerator, giving very simple results. Essentially the same 
method can be used with the other states; for this method is one for treating 
a determinant of Heisenberg’s wave functions, and converting it into a 
determinant of Bloch’s functions; and all of our wave functions are products 
of two such determinants. When we calculate in our case, it appears that the 
terms S will have small effect; we are roughly half way between the cases 1+.S 
and 1—S, and the effects of S nearly average out. This method at the same 
time gives the proper way of considering more distant pairs, as well as adja- 
cent ones; these contribute the further terms in the numerators, as J/+K 
+ ---. Wesee from the next paragraph that these more distant pairs are 
really quite important. 

Method for space lattice. So far, we have spoken about a linear lattice of 
atoms, rather than a space distribution. We now extend this theory to a crys- 
tal; but we shall not carry it through in the same detail. We consider only the 
problem of zero spin, and use our second order perturbation approximation. 
Let us take the body-centered cubic lattice, which the alkalies have. The 
lowest unperturbed state of this lattice can be set up much as with the linear 
one: we let the electrons at the corners of the cubes have the spin a, those at 
the centers the spin 8. Then each electron is surrounded by eight others of 
opposite spin, so that if we consider only adjacent pairs, the diagonal energy 
of this state is 4nJ, where there are » electrons, 4” pairs. This lowest state 
now has non-diagonal terms, each equal to — K, to the 4m states obtained by 
interchanging an adjacent pair. Each of these states has two misplaced spins, 
each surrounded by 7 spins of the same sign, so that the energy has a term 
—14K. Thus for our perturbation problem, we have a non-diagonal energy 
—K, an energy difference 14K, and 4n non-diagonal terms, so that the 
perturbed energy is 4nJ+(4nK?2/14K) =4nJ+(2/7)nK, 

This formula is rather significant. \We compare the energy with that of 
the lowest state of the diatomic molecule, 2 =2, which is (w/2)J/+3nK. We 
observe that for the crystal the coulomb interaction, the term J, has a coeffi- 
cient eight times as great: each atom has eight neighbors instead of one, each 
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penetrating. On the other hand, the valence term K has a coefficient only 
2/7, instead of 1/2. The valence, so to speak, is spread out among all the 
neighbors, and weakened in the process. It is partly on this account that we 
can say that the coulomb interaction is the more important part of the co- 
hesive force, in metals. 

We have considered only those pairs with smallest separation, and they 
give a definite attraction. But in this lattice, there are not only the eight 
nearest atoms at distance R; there are also six, in directions parallel to the 
edges of the cube, at distance of 1.155R, and these have parallel spins, pro- 
ducing therefore a repulsion. In the diagonal energy, each pair will then 
contribute an energy J—A, a positive amount, so that the diagonal energy 
of the lowest state if 4nJ+3n(J(1.155 R)—A(1.155 R)). The next higher 
diagonal energy also will differ from this not merely by —14 A(R), but also 
by an amount 12 K(1.155 R), because by interchange of two spins some of 
these repulsive terms are removed. Thus the lowest energy level, counting 
also these pairs, is 


4nK2(R) 
14K(R)—12K(1.155R) 





4NJ(R)+3n(J(1.155R) —K(1.155R))+ 


This results, on computation, in a much weakened attraction. If we were 
to consider in succession the effects of pairs at greater and greater distance, 
we should come in succession to attracting atoms with antiparallel spin, 
and repulsive ones with parallel, so that the successive approximations to the 
energy would oscillate, falling first above, then below, the true value. 

A pplication to sodium. For the sodium crystal, approximate calculations 
have been carried out, to test these formulas. These were made by taking a 
simple analytical expression for the wave function of the valence electron 
of sodium, and computing the integrals J and K. The details of the calcula- 
tion will not be given here. The first thing that one notices is that, for Na, J 
is several times larger in proportion to K than in hydrogen. It is this fact, 
taken together with the increased coefficient of the J term, that results in the 
importance of the coulomb term. It is also significant in connection with 
the question, why do the alkalies, and metals in general, form metallic lat- 
tices, while hydrogen does not? We can see the essential answer from our 
energy formulas of the previous page. For substances where J is the impor- 
tant term, the coefficient of J will be greater, and the energy lower, for 
the crystal than for the same number of atoms in diatomic molecules, and the 
crystal will be the stable form. For hydrogen, on the other hand, the valence 
term K is the important one. Here the coefficient in the molecular form is 
greater; and even if the metallic form of such a substance were stable in 
the sense of having a minimum of energy for some definite size, as seems 
quite possible, still the energy in the molecular state would be lower. The 
atoms in the crystal would tend to form pairs, resulting in a molecular lat- 
tice; the molecules would repel each other, and would be held together only 
by van der Waals forces, which have been neglected in this paper. This seems 
to be exactly what hydrogen does. 
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The numerical values for Na are approximately as follows. If we take 
only the adjacent pairs, the minimum comes at R=4.9ay approximately, 
rather seriously less than the correct value 7; this can partly be explained 
by the observation that the best atomic wave function for use in the crystal 
would be more extended than that determined from the free atoms, which are 
here used. The energy at this point comes out about —40 kg cal/gm mol, 
the coulomb term supplying about four fifths of this; the observed heat of 
vaporization is 26.4 kg cal, so that this gives, as we should expect, too large a 
value. If now we consider the repulsive pairs at distance of 1.155 R, the 
situation is quite changed. In the first place, the energy is reduced from 
—40 to about —9 kg cal. When we remember that these two values are the 
first two terms of a series, whose value oscillates on both sides of the answer, 
it seems very reasonable that the final result should be not far below the 
experimental value. The problem of properly computing this energy must 
be done by the method, using Bloch’s functions, described in the preceding 
section. In the next place, the minimum of the curve is greatly broadened: 
for quite a range of values, from R=4.9 (the previous minimum) to R=7, the 
energy stays about constant, the change of the attractive term being just 
about balanced by the relatively more rapid change of the smaller repulsive 
effect. (For smaller R’s, a situation can be found when the denominator 
14 K(R)—12 K(1.155 R) =0, so that the function becomes infinite; but this 
is without physical significance.) No doubt a persistence of this effect in the 
final answer helps to correct the improperly low grating space we have al- 
ready found. It also is interesting in connection with the compressibility. 
The alkalies are remarkably compressible, and if we compute the compressi- 
bility for the case where only adjacent pairs are considered, the result is too 
small by a factor of 2 or 3. On the other hand, considering the next set of 
atoms, our very broad maximum would give much too great a compressi- 
bility. Here again it seems that our result may oscillate, perhaps approaching 
eventually something near the right value. 
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ABSTRACT 


The first sections of this paper are a partial repetition and amplification of Born’s 
theory of homopolar crystal lattices, without, however, making any assumption about 
the law of force between atoms other than that the force is central and that diatomic 
molecules are possible. The mathematical work is supplemented with graphical illus- 
trations. If the lattice is in equilibrium as regards (a) displacement of one atom, 
(6) unidirectional extension of the entire lattice, and (c) simple shear of the entire 
lattice, it is shown that each of these postulates leads to a restrictive condition on 
the symmetry and pattern of the possible lattices, thus accounting for the high 
symmetry actually observed. Finally if the equilibria a, 6, c, are stable, there are 
three series of terms in II’, corresponding to elastic constants, which must be posi- 
tive. Making use of these series, it is shown that the lattice is likewise in equilibrium 
against motion of one portion en bloc relative to another portion and hence that 
“mosaic” crystals cannot be configurations of minimum potential energy. 


INTRODUCTION 


HE primary purpose of this paper is to determine whether there is any 

theoretical basis for the existence of a mosaic structure in crystals of 
metals or other crystals containing atoms of only one kind. This periodic 
inhomogeneity for which “mosaic” seems the best short equivalent, has been 
proposed by various workers! as accounting for the non-agreement between 
calculated and observed tensile strengths of crystals. 

Zwicky! and his collaborators, e.g. Evjen,’ have made numerical cal- 
culations of the energy changes due to the formation of these mosaic struc- 
tures in ionic lattices, and find that they are to be expected. Zwicky’ has 
also published a general thermodynamic proof intended to apply to crystals 
of any kind, which the writer* has criticized. The subject will be dealt 
with more fully in the closing paragraphs of this paper. 

The present paper is an attempt to deal with the same problem, for metal- 
lic lattices only, in a rigorous manner. 

The preliminary steps are elementary and a recapitulation of other work; 
but it is hoped that this feature may be of some value to others who have 
found Born’s® classical treatment “nicht ganz einfach zu lesen.” The method 
and notation here introduced will be used by the writer in other papers. 


* Published by permission of the Navy Department. 

1 Notably Zwicky, Phys. Zeits. 24, 131 (1923). 

* Zwicky, Proc. Nat. Acad. Sci. 15, 11 (1929). 

§ Evjen, Phys. Rev. 34, 1385 (1929). 

* Canfield, Phys. Rev. 35, 1 (1930). 

5 Born, “Atomtheorie Des Festen Zustandes,” Leipzig, 1923. 
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FUNDAMENTAL STATEMENTS 


Let us regard the atoms of a homopolar lattice as centers of mutual forces. 
Whatever wave-functions are in the background of this behavior we do not 
enquire. Neither do we concern ourselves with whether or not the atoms are 
charged. We merely assume that all the atoms are alike and that each two of 
them exert mutual forces on each other, independent of the presence of other 
atoms, such that at great distances they attract, and at very small distances 
they repel each other. At very great distances the force of attraction, and all 
its derivatives, vanishes. Finally, the mutual force has a potential. 

We then represent the mutual energy W of two atoms in terms of their 
distance r by a curve like Fig. 1. In the absence of kinetic energy the two 
atoms will form a molecule, their centers being at a distance r,. If they possess 
relative kinetic energy less than W., they will execute vibrations, the ampli- 
tude being represented by p. Due to the energy curve not being in general a 
parabola (i.e. not obeying Hooke’s Law), these vibrations are neither iso- 
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chronous nor simply harmonic. The theory of the vibration of such an an- 
harmonic oscillator is part of the theory of molecular spectra. 

Our present problem is the application of this energy curve to those polya- 
tomic molecules called crystals. 

We next picture an aggregation of atoms in space occupying the points of 
a regular geometric space lattice and forming a crystal. For the time being 
we are not concerned with the exact configuration of this lattice. We fix our 
attention on a single atom of this lattice. We now construct spherical surfaces 
with this atom as center, such that all of the surrounding atoms fall on one or 
another of these surfaces. We select some distance, say the distance to the 
nearest neighbor atom, as a parameter and call it ro. The number of atoms at 
the distance 7» we call 2”;. We designate the shells in the order of their radii 
by subscripts 1, 2, - - - , and their radii by Ryo, kero, - - - Riro; the number of 
atoms in the 7th shell is 27;. ; 

Referring now to Fig. 2, an energy diagram for our selected atom, we put 
on it the ordinates corresponding to the various groups of neighbors at dis- 
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tances kro and label each ordinate with the corresponding n;. (It is plain that 
since the energy of each pair of atoms is W’;, the energy of each atom due to its 
neighbors at kro is only ;W; instead of 2n;IWV;). Then the energy of one atom 
in the lattice is 


U= ails. (1) 
i=l 


We may quite simply refer to 7; as the number of bonds of length hiro asso- 
ciated with each atom in the lattice. 

Still assuming that the type of the lattice is fixed, it is next of importance 
to find how the crystal attains a configuration of a certain density. We are of 
course assuming that the atoms are at rest. Actually they are in vibratory 
motion but it can be shown that the amplitude of vibration is small com- 
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pared with 7» and thus the fixed positions correspond pretty closely with the 
average positions of the atoms over a period of time. It is plain that a change 
in 79 will result in a dilatation of the entire lattice to correspond; the equili- 
brium configuration will be that for which a small change in ro results in no 
change of energy. We therefore differentiate Eq. (1) as follows: 


AU/dro= Don OW /dro= Yonikyi 0 (2) 


where W;’ = |OW/0r],-:,., and OW/dro=ki Wi’. 

In more graphic terms, the array of ordinates in Fig. 2 is horizontally ex- 
expanded or contracted, until, each ordinate being weighted with the corres- 
ponding number n,, their sum is a minimum. 

It will be of interest at this point to study the diagrams in Fig. 3. Diagram 
A is a plot of the number of bonds of length 7 or less than r for an atom ina 
face-centered cubic lattice; diagram B for a body-centered, and C for a simple 
cubic lattice. The dotted lines represent the same functions calculated on the 
assumption of uniformly distributed matter. 

Several points of interest are at once apparent. There are almost always 
more bonds, shorter than any specified length, in a face-centered than in a 
body-centered cubic lattice, and more in the latter than in a simple cubic. It 
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is just this characteristic that is implied by the words “close-packing” applied 
to the two former types. Thus if atoms were rigid spheres attracting each 
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other according to an inverse power of r, and therefore possessing an energy 
curve like Fig. 4, they would inevitably choose one of the configurations of 


W 








Fig. 4 


closest possible packing, i.e., the face-centered cubic or close-packed hexa- 
gonal (practically equivalent). 











Fig. 5 


On the other hand, if the energy curve were like Fig. 5, they would choose 
the body-centered cubic type, since here there is room for the first two groups 
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of the body-centered lattice, containing 7 bonds and thus more negative 
potential energy than the face-centered type, which would have only one 
group containing 6 bonds in the deep loop of the diagram. 

We thus already see a clue to the explanation of the different types of 
lattices. The expression for dl’ in terms of the second power of dro yields a 
formula for the compressibility of the crystal. This is of course positive, since 
the configuration is one of stable equilibrium. The expression is 


dU =3dr,? Donk2w {”. 


EQUILIBRIUM AGAINST VARIOUS DiIsPpLACEMENTS 


1. Displacement of a single atom. \We have already established the manner in 
which the atoms place themselves on a lattice of a certain type and with a 
certain density, i.e., equilibrium as regards dilatation and change of lattice 
type. It is next in order to study equilibrium against a displacement of one 
atom from its correct position, for if the lattice is not in stable equilibrium 
against such displacements, thermal agitation will destroy it. Referring to 
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Fig. 6, the atom O under consideration is shown at the center of the ith 
sphere. Atom O suffers a displacement dx. The atoms in the shell 7 are ar- 
ranged in groups p for which all the atoms in any group p lie in the same plane 
perpendicular to x. The number of atoms lying in the pth plane of the ith 
shell is ”;, and their bonds to O all make the same angle 6;, with the X-axis. 
Evidently ,”;,=2n;. The change in energy dU due to the motion dx is 
then (only taking account of the atoms in shell 7) 


a= (wy > Nip cos op). 
Pp 
For all the shells 7 we then have 


iv=( (Ww? Nip c0s ip) d= 0. (3) 


But we already have an equation (Eq. (2) ) of the form 


> AW =0, 
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and since these are two linear equations in the W,’ their coefficients must be 
proportional or else the coefficients of Eq. (3) must be zero. The latter is the 
only possible alternative* and thus we have for any 7, 


> nip COS O;,=0 (4) 
P 
as a restrictive condition on the configuration. 

Put into words this means that an atom in the lattice must lie at the 
center of gravity of each of the concentric surrounding shells of atoms. Most, 
if not all, of the known crystalline elements obey this condition, so that it 
cannot be called highly restrictive. 

Having disposed of the terms involving dx we now develop dU; in terms of 
dx. Referring again to Fig. 6,it will be seen that dU; = 3(W;'’2,n;, cos? 4;»)dx’, 
in which the coefficients of W;’’ are always positive. But W;’’ may be either 
positive or negative, so we have no a priori knowledge of whether =dU; is 
positive or negative. If negative the lattice is in unstable equilibrium. Since 
we know, however, that actual lattices are actually in stable equilibrium we 
have the strongest reasons for believing that 


Ewe yo nip cos? tip) >0 (5) 
Pp 


t 


and we can hereafter use this series for comparison. 

2. Equilibrium against uniform extension in one direction. \Ve now assume that 
the whole lattice is uniformly stretched in the X-direction at the same time 
not permitting it to contract laterally. Let the amount of stretch per unit 
length be e. Referring again to Fig. 6, we see that bond 7p is stretched an 
amount e€k;rocos’6;,. Hence 


dU; =kiW ero > nip COS? Bip 


Pp 
where 6;, is now confined to values between —7/2 and +7/2. If equilibrium 
exists 


DkWe Soni, cos? 0;>=0. 
) P 


Comparing the above expression with Eq. (2) as before, it is plain that we 
have the same two alternatives; but since the coefficients of W;’ in the above 
expression are definitely positive, each one must be the same multiple of the 
corresponding coefficient in Eq. (2), i.e. 


> Nip cos? 0;,= an; (6) 
P 
where a is the same constant for every different value of 7, but not the same 
for different directions of the axis of elongation. 
This result may be put into more familiar terms by saying that the planar 
moment of inertia of every shell of atoms with respect to a given YZ plane is 


* For Lpnip cos 0p Si, while nj; becomes a smaller and smaller fraction of «ns. 
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the same fraction of k;?ro’”; for that shell. In particular it will be noticed as a 
result of familiar properties of moments of inertia, that the principal axes of 
inertia have the same directions for all shells, and the principal moments of 
inertia are in the same proportion. This constitutes a very severe restriction 
on the pattern of the lattice. 

The terms in & are now to be evaluated. By a similar reasoning to that 
employed in the previous paragraphs, we find 


dU;=3 wire > Nip cos4 Open? 
Pp 


or 
dU= ero? Dkew Yo nip cos! 0; > 
i P 


and since we have experimental evidence that crystals are in stable equili- 
brium against just such an extension as this, we are justified in concluding 


that 

DkeW i” doniy cost 0:>>0. (7) 

i p 
The series in Eq. (7) is proportional to the elastic constant ¢,; for the Y-axis 
we have selected. 
3. Equilibrium against shear. Imagine the lattice distorted parallel to the 
XY plane, the displacements being given by dy=ex, dx =0, dz=0. This isa 
simple shearing distortion. The first order terms in the energy-change are of 
the type 
dU ;=erokW,’ , a cos 8; Sin 8;p > COS dg 
P q 


where ¢, is the angle between the X Y plane and the plane through the X- 
axis which contains the atom in question, and where the summation with 


respect to g means summation over the atoms lying in the circle of inter- 
section of the plane f and the shell 7. This expression may be rewritten: 


Wi 
dUj=— Dx, 


: kro 
where the summation is over the atoms of the 7th shell. For all the shells: 





dU=ery kiwi > — 


r k; ro- 





and if this latter expression equals zero we have: 
DVews > 
i r 


Comparing this equation with Eq. (2) as before, we find the two only com- 
patible if 


XrVe 
kro? 





> «4, = fniki?ro* (8) 
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where f is the same constant for each value of i, but varies with the orienta- 
tion of the coordinate frame. It turns out that Eq. (8) is a consequence of 
Eq. (6), as a simple calculation will demonstrate. Thus, equilibrium against 
shear introduces no restriction on the lattice geometry beyond what was 
already imposed by the equilibrium against unidirectional elongation. (This 
might have been directly inferred from the theory of elasticity). 

The expression for dU in terms of €& must be greater than zero if the equi- 
librium against shear is stable. The corresponding inequality is: 


ke i” z. sin? 6;, cos? 0;, ys cos* @,>0. (9) 
s Pp q 
The series in Eq. (9) is proportional to the elastic constant cg for the axes 
selected. 
4, Equilibrium against displacement of a large portion of the lattice. The princi- 
pal object in deriving the equations of the preceding paragraphs was not to 
prove that crystals can exist, for we know this from observation, but to ob- 
tain certain series in W;’ and W,’’ of which the sums are known to vanish, or 
be positive respectively, and which can be used for comparison with other 
series in W,’ and W,’’ of which the behavior is not known from observation. 
We now draw another picture of our lattice, Fig. 7, in which we once 
more show all the atoms lying on a series of equidistant planes, each plane 
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containing an equal number of atoms. (For some lattices the planes may not 
be equidistant, or equally densely packed, but there will always be a repeti- 
tion of the same distances in the same order, and the argument is but little 
affected). The X-axis is an arbitrary line in space as before. Though the 
transverse planes may be very close together and contain each but a few 
atoms, nevertheless the lattice can always be resolved into such a series of 
planes with finite separation and perpendicular to an arbitrary axis. The 
distance between planes is ary. 

Now we elect to remove the portion of the lattice lying to the right of 0 
by a distance dx from the portion to the left. In so doing we elongate every 
bond crossing the gap dx by an amount dx cos 6;». 
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We must also notice that if k; cos8;,=a, only the atoms lying in the plane 
through 0 have their bonds affected. But if k; cos 0;,=2a, then the plane 
through 0 and also its neighbor on the left contain atoms whose bonds are 
lengthened; and in fact if k; cos6;, = pa there are p planes to the left of 0, and 
including 0, containing atoms whose 7th bonds are lengthened by the opera- 
tion. So if there are .V atoms in each plane 


dU;=NW/i'dx Donipp cos O:p 








P ’ 
But p=k; cos 6;,/a@ hence 
NRW dx 
dl - TE eee > Nip cos* 0: 
a p 
and 
dU =(Ndx/a) Wik Donip cos? Oip (10) 
i Pp 


where 6;, takes on values between —7/2 and +7/2. 
But Eq. (6) told us that Di [Wi'k:z nip cos? 6;> | =(0. Hence in Eq. (10) dU 
=0. Hence the lattice is in equilibrium against the assumed displacement. 
We now develop dU in terms of dx* as follows: 


1 
dU; =SNde i" Dd Nipp COS? Bip 
a Pp 
1 Ndx* Wyk; 
a See os Ynip cos® Op 
2 a » 


dU =(Ndx?/2a) OkiW ie’ Doniy cos* Bip (11) 
t P 





whence 





The question whether the equilibrium we have been discussing is stable 
depends on whether the series in Eq. (11) is positive. | 
To settle this matter we shall require the series of Eq. (7) for purposes of 
comparison. Let us use the letter S for the latter series and the letter S’ for 
the series of Eq. (11). If we now refer to Fig. 2 we see that for the first few 
terms W’;”’ is positive, and for all the remainder out to infinity W,’’ is nega- ‘ 
tive. Let us say that W,”’ is the last positive term. Then we may write S= 
Si-—Spez where S, is a terminating series ending with W,’’, and Sz is an 
infinite series of positive terms starting with —W,4,'’.. The fact that S>0 
implies that S,4>S,. Similarly we divide S’ into two parts S,4’ and S,’ the 
division being at the same term as for S. Let the 7 term of S be S; and the z 
term of S’ be S,’. Then , 


k?2w dnp cos4 bin 

kW’ Doni, cos? i> 
kW!" Do pniy cos 0;> 

kw” > Nip cos Bip ) 





S; 
5! = 
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Let S;/S;’=t;. Evidently a<t;<aqg; where gq; is the largest integer in 
k./a, i.e., the largest value of p. Also t; <ti41.* 
We may now write the two series as follows: 


u i) 
S = Sy —Spz = Dts _— Dus,’ 
1 


ut+l 


m a 
S’ =S4'—S,'= > Si/- Ea. 
1 


w+ 
Then 
Qa <S4/S4' <t <aq, <= Adyu+1 <Sp/Sp' ‘ 
or 
S4/Sa ‘ <Sp/Sp' ‘ 
But 
Sa>Sp 
hence 
Sp/Sa’ <Sp/Sp’ 
or 


S4'—Sp'>0. 


The expression in Eq. (11) is therefore positive, and the equilibrium is 
stable. In other words, any displacement of a portion of the perfect lattice 
directly away from another portion sets up a restoring force. If the crystal 
does not extend to infinity to right and left, the effect will only be to make 
Sp and S,’ terminating series, and the inequalities will be even more pro- 
nounced. It only remains to carry through a similar argument with respect to 
a lateral displacement of one portion of the lattice relative to the other. It 
turns out that, using a similar notation and procedure to that which pre- 


ceded, 
S= Lk2ws" Y cos? Gy sin? 6ip DY, cos? bg (Eq. 9) 
i Pp q 


=a kW." Yop cos Gip sin? Bip Do cos* $q 
i p ; 


and 
S'= Yk We" Dd cos jy sin? 6;, >, cos? dg. 
i P q 


From which, as before, we find that dU >0 and the equilibrium is stable with 
respect to this displacement too. 


POssIBLE CLASSES OF SYMMETRY IN HOMOPOLAR CRYSTALS 


In this section we shall discuss rather briefly the consequences of Eqs. (4) 
and (6), so far as they affect the symmetry properties of homopolar lattices as 
defined in the opening paragraph. 

The meaning of Eq. (4) is that any atom in the lattice must lie at the cen- 
ter of gravity of any group of atoms which are situated on the surface of a 


* For large values of k, ¢ approaches (4/5) k; there may be isolated values of i for which 
ti > tis, but these are very rare. It is impossible to imagine a W function for which this slight 
departure from rigor would invalidate the general conclusion. 
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sphere with the atom in question at its center. It is plain that the atoms of 
such a group must display at least all the elements of symmetry of the crystal 
as a whole. (They may display additional elements of symmetry, but here 
we are concerned with minimum requirements.) 

The elements of symmetry, so far as they affect our present argument, 
are: (1) a center of symmetry; (2) an axis of symmetry, either digonal, tri- 
gonal, tetragonal or hexagonal; (3) a plane of symmetry. The properties of 
these as they affect the center of mass are as follows: 

1. The center of gravity coincides with a center of symmetry if one exists. 

2. Any axis of symmetry contains the center of gravity. 

3. Any plane of symmetry contains the center of gravity. 

Since all elements of symmetry of points on a sphere contain the center of 
the sphere, it is sufficient that the elements of symmetry shall determine a 
point, to make that point the center of gravity. This then may be accom- 
plished in the following ways: 

(1) A center of symmetry. 

(2) Two axes of symmetry of any order. 

(3) An axis of symmetry and a plane of symmetry which does not contain 
the axis. 

(4) Three non-collinear planes of symmetry. 

Any class of symmetry, then, which displays at least one of the four com- 
binations of symmetry elements listed above, fulfills the requirements of 














TABLE I.f 
Class Symbol Class Symbol 
(Dana) (Schoentlies) (Dana) (Schoentflies) 
I. Triclinic System V. Cubic 
Asymmetric Ci *Tetartohedral i 
*Normal Ci *Pyritohedral T* 
II. Monoclinic System *Tetrahedral T¢ 
Clinohedral C; *Plagihedral U 
Hemimorphic C2 *Normal O* 
*Normal C2 VI. Hexagonal. 
UI. Orthorhombic (Rhombohedral Division) 
Hemimorphic C.” *24 C. 
*Sphenoidal V *Trirhombohedral C3 
*Normal ys Ditrigonal Pyramidal C3” 
IV. Tetragonal Trapezohedral D; 
Tetartohedral S4 *Rhombohedral D;4 
*Sphenoidal ya (Hexagonal Division) 
Pyramidal Hemimorphic C, *23 C3" 
*Pyramidal CH *Trigonotype D;* 
Hemimorphic Cy Pyramidal Hemimorphic C, 
*Trapezohedral Ds *Pyramidal Cs" 
*Normal Dy? Hemimorphic Ce 
*Trapezohedral Ds 
*Normal Deo" 








*Indicates classes which conform to barycentric symmetry condition of Eq. (4) in text. 
t Table copied from Wykoff: “The Structure of Crystals”, N. Y. 1924. 


Eq. (4). In Table I are listed the 32 recognized classes of crystal symmetry. 
Those classes which obey Eq. (4) are designated with asterisks. Of the 32 
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classes, 21 obey this condition. All the known space lattices of crystal ele- 
ments fall into the 21 classes with the possible exception of graphite, of which 
the lattice is still the subject of discussion. 

We still have to discuss the consequences of Eq. (6), which imposes a 
restriction far more severe than Eq. (4). 

Eq. (6) demands that the principal axes of inertia of each sphere of atoms 
concentric with a selected atom shall have the same directions (this also fol- 
lows from symmetry considerations) and that the principal planar moments 
of inertia shall be in the same proportion. So far as the writer can determine, 
this condition is satisfied by the five classes of the cubic system and none 
other. The hexagonal close-packed arrangement may be said to satisfy the 
condition ‘approximately’, i.e. the first departure is in the shell 7=3, and 
involves the transfer of two atoms from a shell containing twelve, to a shell 
of about 6 percent smaller radius. Just what meaning can be attached to this 
‘approximate’ fulfilment of such a condition is not clear. It may mean that 
this small departure from equilibrium can be balanced by a minute polar 
force-field superposed on the central force-field which we have postulated 
in the first paragraph of the Fundamental Statements above. 

It may finally be remarked that if Eq. (2) has a solution ro for a particular 
lattice type, and if this lattice type satisfies Eqs. (2) and (4), and if =",;W; 
is smaller for this configuration than for any of the other configurations 
satisfying Eqs. (1), (2) and (4), then the inequalities expressed by Eqs. (5) 
and (7) show that no lattice obtained by a homogeneous deformation of the 
given lattice can satisfy the given conditions. Moreover, since the validity of 
these inequalities is only increased by leaving out terms beyond the first few 
terms in the summations, all departures from complete geometrical perfection 
in the interior of a homopolar crystal are unstable. 

At points near the surface of a crystal, particularly for points nearer to 
the surface layer than the distance at which W,’’=0, the above statement 
is not true; the writer hopes to make those surface conditions the subject of 
a future paper. 


Mosaic STRUCTURES 


The displacement of one portion of a crystal en bloc relative to the other 
portion may be said to be the purest type of inhomogeneous deformation. 
From the two cases of this kind which we have considered above, in conjunc- 
tion with the postulated stability against homogeneous deformation, it is 
possible to show that any conceivable picture of inhomogeneous deformation 
leads to an increase of energy in the second order terms of the displacements. 
Among others, it is possible to show that the equilibrium is stable against a 
uniform or non-uniform contraction of a single lattice plane. 

Two types of variation from perfect symmetry are conceivable mosaic 
structures. The first is obtained by some sort of inhomogeneous deformation 
of an originally perfect lattice. Our conclusions given above show that any 
such configuration involves a greater energy than the original. 
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The second is obtained by some such inhomogeneous deformation forming 
gaps in the lattice into which additional atoms are inserted. To show that 
this results in a preferred configuration it is necessary not merely to show (a) 
that the new atom has a negative potential energy greater than the positive 
energy created by the deformation, but also (b) that the net decrease in 
energy thereby effected is greater than what would result from adding one 
more atom in a normal way to the structure of the perfect lattice. Zwicky* 
in replying to the author’s‘ criticism his general thermodynamic argument 
for all crystals, recapitulates his calculations for an ionic lattice (rock salt), 
and establishes the condition (a) above; but he does not show that con- 
dition (b) is satisfied as well. 

Evjen® has calculated the case of inhomogeneous deformation of a rock- 
salt lattice and finds as we do here for a homopolar lattice, that the energy 
terms in the second powers of the displacements are positive. He does, 
however, find that when terms of the third and fourth powers are included, 
it is possible to find a secondary minimum in the energy which may be lower 
than the original minimum. 

In concluding the author wishes to express to Professors kK. F. Herzfeld 
and F. D. Murnaghan his appreciation of their comments on certain por- 
tions of this paper. 


* Zwicky, Phys. Rev. 35, 3 (1930) 
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ABSTRACT 
The temperature coefficients of polarization capacity and resistance have been 
determined for gold and plain platinum electrodes in 0.63, 1.50 and 2.42 N solutions 
of sulphuric acid and also at different frequencies for 1.46 N concentration. To 
reduce the electrodes to a condition representing minimum capacity, a method 
of thermal treatment was employed and found quite effective. Values of the angle 


of phase difference have been computed for the various conditions under which 
measurements were made. 

During the process of reducing the electrodes to a condition representin g mini- 
mum capacity, the value of polarization capacity decreased to about one-half the 
original value and the resistance increased to a few times the original value. The 
angle of phase difference increased 15 or 20 degrees. After the electrodes were con 
ditioned, consistent and reproducible values of capacity and resistance could be 
obtained. 

The temperature coefficient of capacity is about twice as great with platinum 
as with gold electrodes and is apparently independent of concentration. Phase 
difference changes only slightly with temperature, but increases with concentration. 
The absolute values of the slopes of the curves, capacity versus temperature and re- 
sistance versus temperature, decrease with increasing frequency. The temperature 
coeflicient of capacity, however, increases with frequency for platinum electrodes 
and decreases in the case with gold electrodes. 


HERE has been a limited amount of investigation of the effects on 

polarization capacity and polarization resistance caused by changes in 
temperature. Sheldon, Leitch and Shaw! determined the temperature coefti- 
cient of capacity for cells consisting of Pt |H2SO, | Pt and Hg |H2SO, | Pt. 
Schénherr? studied the influence of temperature on the reduction of the elec- 
trodes to a condition corresponding to minimum capacity. The work herein 
described was an investigation to determine the nature and magnitude of 
changes in polarization capacity and resistance as influenced by changes in 
temperature, under various conditions with respect to the concentration of 
the electrolyte, the frequency of the applied voltage and the material and 
condition of the electrodes. 


APPARATUS AND EXPERIMENTAL PROCEDURE 


Measurements were made by means of the Wheatstone bridge and 
auxiliary apparatus, diagrammatically represented in Fig. 1. The low- 
power vacuum tube oscillator which had a frequency range from 600 to 4000 


' Sheldon, Leitch and Shaw, Phys. Rev. (1) 2, 401 (1895). 
* Schénherr, Ann. d. Physik u. Chemie 6, 116 (1901). 
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cycles per second was loosely coupled to the bridge by means of a step-down 
transformer in order to maintain the voltage across the polarization cell at 
or below the value corresponding to minimum or initial-capacity. Minimum 
or initial-capacity is that value of polarization capacity characteristic of the 
cell when the applied voltage is sufficiently low that further decrease in 
voltage does not result in any appreciable decrease in the value of polar- 
ization capacity. With respect to variations in capacity as influenced by 
variations in voltage, Wien® found that the polarization capacity was constant 
at the minimum value if the voltage across the cell was 0.05 or less. Wolff* 
found the capacity nearly constant for voltages less than 0.1 volt. In the pres- 
ent work, in order to ascertain whether or not variations in applied voltage 
were influencing the measurements, the voltage across the polarization cell 
was measured by means of a vacuum tube voltmeter and was found to be 
less than the critical value. Also, resistance was inserted in the leads to the 
bridge without decreasing the measured value of capacity. All values of 
polarization capacity recorded in this paper are, therefore, thought to be 
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Fig. 1. Bridge and auxiliary apparatus. 
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minimum values insofar as the influence of the voltage applied to the cell is 
concerned. 

The resistances, R; and Ro, in the ratio arms of the bridge were main- 
tained equal in order to minimize the effects of distributed capacity and 
inductance although with the type of bridge used and at the comparatively 
low values of frequency, effects of distributed capacity and inductance were 
negligible. The condenser, Ci, was a fixed condenser and served to prevent 
any direct current through the cell. 

Fig. 2 is a diagram of the polarization cell. The cell consisted of a con- 
tainer in the form of a test tube, with a pair of similar electrodes mounted 
by inserting through the perforated stopper the glass tubes into which the 
wire electrodes were sealed. The thermometer was similarly mounted with 
the center of the bulb in the same horizontal plane as the faces of the elec- 
trodes, the thermometer and electrodes being securely fastened together so 
that the faces of the electrodes were coplanar and about 3 mm apart. The 
end of any glass tube containing an electrode was ground transversely to the 
length so that the exposed circular cross-section of the wire electrode was 

3 Wien, Ann. d. Physik u. Chemie 58, 37 (1896). 

4 Wolff, Phys. Rev. 27, 755 (1926). 














POLARIZATION CAPACITY AND RESISTANCE 545 


flush with the end of the tube. Gold and plain platinum electrodes were 
used, the gold having a diameter of 2 mm and the platinum a diameter of 
0.623 mm. Under the conditions of the experiment, the capacity of polariza- 
tion cells with electrodes as described was within the range from 0.03 to 0.7 
microfarads and the corresponding range of polarization resistance was from 
4000 to 500 ohms. For any particular test the electrolyte was an aqueous 
solution of sulphuric acid of desired concentration, which had been heated 
to free it of absorbed gases. The temperature of the cell was controlled by 
immersing it in an oil or water bath, the temperature of which could be 
controlled at will. 














LY 


Fig. 2. Polarization cell. E, electrode. T, bulb of thermometer. 


Data were obtained for any one set of conditions by making a series of 
measurements as the temperature was increased through a range from 0°C 
to a maximum of 90 or 95°C and then decreased through the same range. 


THE EFFECT OF THERMAL TREATMENT ON STABILIZATION OF 
THE ELECTRODES 


As have other experimenters,® the writer found that cells having freshly 
surfaced electrodes do not give consistent and reproducible values of capacity 
and resistance. Incidentally, in attempting to duplicate results for rising and 
falling temperatures of a cell with freshly surfaced® electrodes, it was found 
that the effect of any so-called temperature run such as that described above 
was to decrease considerably the values of capacity and increase the values 
of resistance. The value of the capacity at 0°C after the cell had been taken 
through the range of temperatures was reduced to a value which was from 
40 to 60 percent of the original value at that temperature on starting the run. 
More marked changes took place in the resistance, the final value of which 
might be as much as four times greater after the conditioning runs. After 


§ Gordon, Géttingen Diss. 72 (1896-97); Orlich, Diss., Berlin (1896); Wolff, reference 4. 
6 The electrodes were surfaced by filing with a fine-grained file. 
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two or three such runs with a cell, values of capacity and of resistance ob- 
tained with falling temperature were in agreement with those obtained with 
rising temperature. Schénherr,’ also, found by increasing the temperature of 
the cell through a considerable range and then decreasing it through the 
same range that the final value of capacity was much less than the original 
value. 

The manner in which polarization capacity and resistance change for 
freshly surfaced platinum electrodes when the cell is subjected to a series of 
runs is shown by the curves of Fig. 3. Curve A shows the relation between 
capacity and temperature for the initial run and curve B shows the relation 
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Fig. 3. Capacity-temperature and resistance-temperature curves for platinum electrodes 
while in the process of being conditioned. For explanation of curves, see text. Values of capacity 
are in microfarads per sq cm for one electrode and the plotted values of resistance represent 
total resistance of the cell. To obtain resistance per sq cm for one electrode, the values as 
plotted should be multiplied by 0.0015. Open circles, without and with crosses, represent 
values of capacity for increasing and decreasing temperature, respectively, while filled-in 
circles represent values of resistance, circles with projections being values obtained with de- 
creasing temperature. For all curves, the concentration was normal and the frequency was 
1020 cycles per second. 


for the third run, the second having given results nearly the same as the 
third. The corresponding curves showing the development of polarization 
resistance’ are A’ and B’, respectively. It may be noted in curve B (capacity) 


7 Schénherr, reference 2. 

’ Values of resistance are in all cases the measured values of the total resistance of the 
polarization cell without deducting for the resistance of the electrolyte which was assumed 
negligibly small in comparison with the resistance of polarization. A test of the validity of this 
assumption was made by varying the spacing between the electrodes and noting the effect 
on resistance. When the spacing was increased to 16 times the original value, the resistance 
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and in curve B’ (resistance) that the agreement between the values of 
capacity and of resistance for rising temperature and the corresponding 
values for falling temperature is reasonably good. This same fact is also well 
emphasized in the curves of Fig. 4, in nearly all of which the values obtained 
with rising temperature are in very good agreement with those obtained with 
falling temperature. The one prominent exception in the curves of Fig. 4 is 
in the case of curves C (capacity-temperature for 0.63 N solution, with 
plain platinum electrodes) and C’ (the corresponding curve for resistance). 
The lack of agreement in this case is due to the fact that these curves repre- 
sent the first run of the series with the platinum electrodes and there had 
been only one run made to condition the electrodes previous to taking the 
data for these curves. The one run was not sufficient fully to condition the 
electrodes. In fact, in curves B and B’ which are plotted from measurements 
made during a run following that for curves C and C’, there is still lack of 
complete agreement between the rising and falling-temperature branches 
of the curves. 

When cells with gold electrodes were treated similarly, the effect was 
similar, excepting that the gold electrodes were more readily conditioned 
and did not require as many conditioning runs as did the platinum electrodes. 

Success in duplication of results in making measurements of polarization 
capacity and resistance was possible only when the electrodes were fully con- 
ditioned. The writer found the method of thermal treatment as described 
above to be quite effective in accomplishing this and it is, moreover, a 
method which requires much less time than that of leaving the electrodes in 
the electrolyte to age for several days. Just what the condition of the elec- 
trode is when it has minimum capacity for a given set of conditions, is not 
definitely known. It doubtless represents a certain condition with respect 
to absorbed gas and also with respect to the condition of the surface. As 
have other experimenters, the writer observed that the electrodes undergo 
a change in color as they become conditioned, changing from a bright to 
a dull color. It was also observed that the electrodes could be conditioned 
in an electrolyte of given concentration and then be used in an electrolyte 
of different concentration and yet have capacity according to the second 
concentration. 

In curves A and A’, Fig. 3, on the rising-temperature branches of the 
curves, attention is called to the fact that at a temperature of about 31°, 
the capacity underwent a rather sudden decrease and the resistance a sudden 
increase. The cause for these changes was a drop of about 1° in the tempera- 
ture of the cell at this point after which, although the temperature was 





increased only 3.5 percent. Also, the platinum electrodes which were used during the experi- 
mental work were finally platinized and the total resistance was determined. The value thus 
determined, when corrected for concentration, was not in any case more than 5 percent of the 
total resistance of the cell with plain platinum electrodes. Since part of the resistance of the 
cell as measured with platinized electrodes was resistance of polarization for platinized platinum, 
the value of the resistance of the electrolyte thus determined was somewhat too large. Con- 
sequently, 5 percent is a liberal estimate of the outside limit of the approximation. 
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again brought to its former value, the values of capacity and of resistance 
had undergone the changes noted. Later, in order to investigate this effect 
further, the electrodes were again freshly surfaced and the data represented 
in curves C and C’ were obtained. In taking the data for these curves, the 
temperature was purposely permitted to fall 1 or 2 degrees at three different 
points on the ascending parts of the curves. The curves show that apparently 
the reversal of the temperature gradient with reference to the electrode and 
the electrolyte is effective in decreasing the capacity and increasing the re- 
sistance and also that a change in capacity is accompanied by a corresponding 
and closely related change in resistance. 

During the conditioning process in case of the cell with platinum elec- 
trodes the phase difference, i.e., tan~'wRC, where w is 27Xfrequency, R is 
polarization resistance and C is polarization capacity, increased from 27° to 
40° and in case of the cell with gold electrodes the phase difference increased 
from 12° to 35°. The initial low value of phase difference in the latter case 
is probably due to the fact that the gold electrodes had previously been 
heated to the melting point, while in the case with the platinum electrodes, 
the electrodes had been freshly surfaced by filing, only. The electrode, then, 
having the property of polarization capacity, may be regarded as becoming 
a less perfect condenser as the electrode becomes conditioned or as it ac- 
quires the condition corresponding to minimum capacity. 

The writer has chosen in this paper to regard polarization capacity and 
the resistance associated with it as being equivalent to a pure capacity 
and a pure resistance in series. Successful attempts were made, however, to 
balance the bridge both with a resistance in series with a mica condenser 
in the arm of the bridge opposed to the polarization cell, and also with a 
resistance in parallel with the standard condenser. The results were as fol- 
lows: 

Measured values of capacity and resistance when the cell was balanced 
by a capacity and resistance in series were, respectively, 


0.0515 wf and 2840 ohms. 


When C,, and R,,, signifying the corresponding values for the multiple 
relation, are computed from the above values, 


‘s 


a. 
140°C eR? uf 


and 


27. 2R 2 
WU s Ns 


1 
R,.= R(1 +=) =6070 ohms. 


The values of C,, and R,, as measured directly by means of the bridge 
were, respectively, 0.0292uf and 5920 ohms, in reasonably good agreement 
with the;computed values. The balance or null point was quite distinct in 
both cases. 
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TEMPERATURE COEFFICIENT AND CONCENTRATION OF THE ELECTROLYTE 


With gold and plain platinum electrodes, conditioned according to the 
method described above, data were taken for different concentrations of the 
electrolyte as the temperature was varied over the complete range. The data 
so obtained are plotted in the curves of Fig. 4, the variation of capacity and 











Fig. 4. Curves, capacity-temperature (letters not primed) and resistance-temperature 
(corresponding, primed letters) for different concentrations of the electrolyte. Curves A, A’, 
D and D' are for 2.42 N concentration; curves B, B’, Eand E’ are for 1.50 N; and curves C, 
C’, Fand F’ for 0.63 N. The six upper curves are for platinum electrodes and the six lower 
for gold. For the significance of different kinds of points, see explanation, Fig. 3. The fre- 
quency was constant at 1000 cycles per second. All values of capacity are in microfarads per 
cm? for one electrode, but values of resistance represent the total resistance of the cell. To ob- 
tain resistance per cm? for one electrode, the plotted values should be multiplied by 0.0015 for 
the cases with platinum electrodes and by 0.0174 for the cases with gold electrodes. 


the corresponding variation of resistance with temperature for the three 
different values of concentration being shown for the two kinds of electrodes. 

Table I shows resuits of computations based on the curves of Fig. 4 and 
the data relating thereto. As is indicated in the table, the temperature 
coefficient of capacity, defined as 


(1/Cis)(dC/dT))s, 


is about twice as much for platinum electrodes as for gold, other conditions 
being the same. With possibly an exception in the case with platinum 
electrodes for the 2.42 normal solution, the temperature coefficient for both 
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TABLE I. Temperature coefficients of polarisation capacity and polarization resistance for different 
concentrations of the electrolyte. 


Kind of Electrodes | Plain Platinum — Gold 
Concentration of : | i | [ pron 
Electrolyte 0.03 N 1.50 N 2.42N | 0.63 N | 1.50 N 2.42 N 
om | 28.5 44.3 55.3 7.0 10.0 12.6 
(dC dT), | 0.29 0.448 0.466 0.0372 0.053 | 0.067 
(1/Cs) dC dT), | 0.0101 0.0101 0.0084 0.0053 0.0053 | 0.0053 
Ry 4.3° | 3.5 2.9 | 16.7 13.8 | 12.2 
(dR/dT)s —0.045 | —0.032 | —0.026 | —0.007 | —0.043 | —0.040 
(1/Ri) (dR'dT),5 | —0.0092*1 —0.0089 | —0.0089 | —0.004 | —0.0033 | —0.0033 


Notes: Values of capacity are in microfarads per cm? for one electrode. Values of resistance 
are in ohms per cm? of area for one electrode. N signifies normal solution (H2SO, in water), 
i.e., 49 grams of acid per liter of solution. The frequency was constant at 1000 cycles per 
second. 

* This value is somewhat uncertain, due to the fact that the electrodes were not fully 
conditioned. 

Relatively, the values of Cys and Rjs in Table I are accurate to within about 1 percent. 
Values of (dC/dT),s and (dR/dT),s, and, consequently, the values of temperature coetticients, 
are probably not more accurate than to within 5 percent. 


gold and platinum electrodes is independent of concentration for the three 
values of concentration involved. In other words, Cis and (dC d7))\s change 
in the same manner with respect to changes in concentration. 

Also, it may be seen that the values of resistance per unit area of elec- 
trode, Ris, with platinum are about one-fourth the values for the corre- 
sponding concentrations with gold electrodes. On the other hand, values of 
capacity per unit area are about four times as large with platinum as with 
gold electrodes for corresponding concentrations. It may be noted from the 
curves in Fig. 4 that capacity and resistance do not vary linearly with 
temperature in the cases with platinum electrodes as they do with the gold 
electrodes. 

The value of temperature coefficient of capacity as obtained in this work 
for platinum electrodes is somewhat less than that computed from the data 
given by Sheldon, Leitch and Shaw.’ The value computed from their data is 
about 0.02 per degree centigrade when compared to the capacity at a tem- 
perature of 18°C. 

The degree of variation of the product, capacity of cell times total 
resistance of cell (regarded as the polarization resistance), is shown in Table 
II. The corresponding values of the angle of phase difference are also shown 
in many instances. Although in Table II, the capacity and resistance of 
the cell were used rather than the values for a single electrode, nevertheless, 
regarding the capacity and resistance as being in series, the product would 
not be changed in value in considering the results for a single electrode, for 
in that case the value of capacity would be double the total capacity of the 
cell and the value of resistance would be one-half the total resistance of the 
cell. Since the frequency was constant, Table II also indicates the variations 
in the tangent of the angle of phase difference for the different conditions 
represented by the several curves. It may be noted that there is no large 


* Sheldon, Leitch and Shaw, reference 1. 
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variation in the phase difference for any one temperaturerun. With gold elec- 
trodes, the phase difference seems to increase slightly with rising tempera- 
ture, while with platinum electrodes, for the higher concentrations especially, 
the phase difference seems to be slightly lower at the higher temperatures. 
For both kinds of electrodes, the phase difference increases with increase in 
concentration of the electrolyte. 


TaBLe II. Variation in product, C XR, with temperature, frequency being constant. 








Kind of electrodes Plain Platinum Gold 


Concentration 0.63 N (H2SO,) | 1.50 N (H2SO,) | 2.42 N (H2SO.) | 0.63 N (H:SO,) | 1.50 N (H:SO,) | 2.42 N (H:SO,) 








134 (40°10’) 158 (44°48’) 162 (45°34’) 113 (35°30’) 131 (39°27) 142 (41°52’) 
137 159 163 115 133 





























147 
139 160 161 116 136 153 
Temperature 140 (41°20’) 160 (45°15) 162 116 139 156 
increasing 140 158 159 (45°3") 117 (36°20’) 140 (41°30) 159 
139 157 158 118 141 161 (45°22’) 
139 156 158 118 142 160 
138 154 157 118 (36°30’) 143 (41°57’) 159 (45°2") 
Maximum 138 (41°5’) 152 (43°40) 155 (44°15’) 117 142 
temperature — 
139 152 154 111? 141 167? 
141 153 157 115 141 159 
Temperature 142 155 159 114 (35°43’) 140 (41°20’) 157 
decreasing 143 (42°5’) 156 161 (45°27') 114 138 155 
143 158 161 114 137 153 
141 160 163 114 135 (40°22’) 150 
139 161 (45°18") 163 113 133 147 
137 (40°50’) 160 164 (45°51") 114 133 145 (42°20’) 
Average 139 (41°00’) 157 (44°27’) 160 (45°10’) 115 (36°00’) 138 (40°54’) 154 (44°8’) 














- 


re 





Note: The average values of the angle of phase difference are computed from the average values of (C XR). 
In Fig. 5 are curves showing the results of attempted runs after the 


electrodes had been polarized with oxygen or hydrogen, as indicated in the 
figure. The electrodes were polarized by connecting the two together and 
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Fig. 5. Curves, capacity versus temperature, for polarized electrodes. The kind of 
electrodes and type of polarization are indicated in the figure. In all cases, the concentration 
was 1.46 N and the frequency was 1000 cycles per sec. Values of capacity are in microfarads 
per cm? for one electrode. 


to the appropriate terminal of a direct current source of potential while a 
comparatively large third electrode was connected to the other terminal. 
Excepting in the case of gold polarized with oxygen, the capacity increased 
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regularly on starting the runs, but at temperatures of about 35° or higher 
there were sudden and irregular decreases in the value of capacity. Such 
decreases usually occurred simultaneously with the appearance of gas on 
the faces of the electrodes. Continuing the runs throughout the range of 
temperatures, the final values of capacity were much lower, yet still some- 
what higher than the minimum values. 


TEMPERATURE COEFFICIENT AT DIFFERENT FREQUENCIES 

Curves, capacity versus temperature, representing measurements at 
different frequencies, but with the concentration constant, are very similar 
to those of Fig. 4. The slopes of the curves decrease considerably with in- 
creasing frequency. For platinum the temperature coefficient increases with 
increasing frequency, while the opposite seems to be true in the case with 
gold electrodes. The results of such tests are given in Table III. As may be 
seen from the table, both capacity and resistance decrease markedly with in- 
creasing frequency (see under Cis and Rjs). 

The variation of polarization capacity and resistance with frequency has 
been the subject of special investigations by Merritt,!° Jolliffe," and Wolff.” 
The results found in this work for variation of capacity with frequency, 
although representing a comparatively narrow and somewhat lower range 
of frequencies, are in general agreement with the results of their investiga- 
tions. 


TABLE III. Te aperaiane coefficie nt of capacity at diffe rent frequencies. 


















































Kind of Electrode Plain Platinum | Gold 
Frequency 645.0 /1290.0 2581.0 3872.0 645.0 2581.0 
Cis 60.6 42.4 28.2 21.3 14.7 7.6 
(dC/aT):s 0.487 0.356 0.250 | 0.205 0.111 0.048 
(1/Cys) (dC/dT)1s 0.0080 0.0084 0.0088 0.0096 0.0075 0.0063 
4.56 3.18 2.0 1.9 15.5 6.5 
(dR/dT):s —0.034 | —0.028 | —0.015 | —0.015 | —0.083 | —0.029 
(1/Rus) (dR/dT)1s | —0.0075 | —0.0088 | —0.0075 | —0.0088 | —0.0053 | —0.0045 














Notes: Values of frequency are in cy is per ae Values of capacity are in microfarads 
per cm? for one electrode. Resistance is in ohms per unit area of electrode. The concentration, 
1.46 normal, was the same in all cases. 


DIsCUSSION 

It is very well known that temperature affects in a very important way 
most electrolytic processes. The temperature coefficient of conductivity of 
sulphuric acid of normal concentration as given in tables by Kohlrausch 
and Holborn" is about 0.012 per degree centigrade at a temperature of 22°C. 
Also, the temperature coefficient of fluidity as based on data from tables is 
positive with a value of about 0.028 per degree centigrade at 15°C. In con- 
sideration of the facts that the temperature coefficients of capacity for plati- 
num and gold in sulphuric acid are not the same and that the variation of 


1° Merritt, Phys. Rev. 17, 524 (1921). 
1 Jolliffe, Phys. Rev. 22, 293 (1923). 
12 Wolff, reference 4. 

48 Kohlrausch and Holborn, Leitvermégen der Electrolyte (1916). 
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capacity with temperature is not linear in the case with platinum, while it 
is in the case with gold, the temperature coefficients of capacity and resis- 
tance are not fully accounted for in variations of conductivity and fluidity 
with temperature, for if so the coefficients for the two metals would be af- 
fected alike. 

Knobel and Joy found, with smooth platinum, that increasing the 
temperature had the effect of decreasing the overvoltage. Overpotential 
and catalytic activity of metals, according to the work of Rideal, vary 
inversely, i.e., metals arranged in order of decreasing catalytic activity are 
also arranged in order of increasing overpotential. Bancroft,’ also, states 
that the metals which have a lower overvoltage are those that catalyze the 
reaction, 2H—Haz, while those which havea high overvoltage do not. Among 
other things, Bancroft also discusses the action of the so-called catalytic 
poisons in preventing catalysis. Now Wolff'? has shown that catalytic 
poisons have the effect of decreasing polarization capacity. It seems quite 
likely, therefore, that overvoltage, catalytic activity, and polarization 
capacity are related in an important and probably a very complex way. 
The decrease in capacity of the electrode as it acquires its stable condition 
after polishing would mean a corresponding decrease in its catalytic activity. 
The increase in capacity with temperature would mean an increase in 
catalytic activity and a reduction in the conditions favorable to overvoltage, 
as found by Knobel and Joy. It is recognized that in considering polarization 
capacity in relation to overvoltage one is dealing with much lower values 
of voltage in the case of polarization capacity, and for so-called minimum 
capacity, values of voltage which are thought to be strictly reversible. 

If polarization capacity and catalytic activity are intimately related, 
then the results of this paper do not seem to be in full agreement with the 
work of Bowden and Rideal'® who found that the specific catalytic activity 
of the various metals, including platinum, for the electro-deposition of hy- 
drogen suffer only small variations by thermal treatment. Here again, how- 
ever, the full conditions are not the same, for in this work the capacity was 
determined by means of a bridge supplied with alternating current and, con- 
sequently, not only is the deposition of hydrogen involved, but also the 
deposition of oxygen and the effect of one on the other after consecutive 
reversals of the current. 

In conclusion, the writer desires to express appreciation of the excellent 
facilities available at Cornell University, Department of Physics, for doing 
this work. Thanks are extended to Professor E. Merritt, who approved this 
problem and under whose direction the work was done, for his interest and 
many helpful suggestions. The writer is also indebted to Professor C. C. 
Murdock for his interest in the results of the work and a discussion of their 
relation to other problems of polarization capacity. 

14 Knobel and Joy, Trans. Am. Electrochem. Soc. 44, 443 (1923). 

18 Rideal, Jour. Am. Chem. Soc. 42, 94 (1920). 

16 Bancroft, Trans. Am. Electrochem. Soc. 37, 21 (1920). 

17 Wolff, unpublished paper. 

18 Bowden and Rideal, Proc. Royal Soc. of London A120, 80 (1928). 
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(Received January 13, 1930) 
ABSTRACT 


For diffusion in substances for which the conductivity and diffusivity are func- 
tions of temperature, the partial differential equation which determines the tem- 
perature is non-linear in character. For a certain class of boundary conditions the 
introduction of an auxiliary function of time and distance reduces this equation 

. to an ordinary differential equation of the second order, in which the temperature 
and the auxiliary function appear as dependent and independent variables respec- 
tively. As this resulting equation can be solved by methods of approximation, the 
solution to the original equation can be thus obtained. 

The conditions under which such an auxiliary function exists are examined and 
a simple criterion developed by means of which it may be determined whether such 
a function exists for any particular case. The method is shown to apply to the case 
of the semi-infinite solid, and the solution of a particular problem of this type (the 
cooling of a steel ingot) is given. Another case to which the method applies (diffusion 
in an axially heated cylinder) is considered, and shown to afford a possible method for 
determining diffusivity as a function of temperature. It is noted that the method de- 
veloped is available for the solution of problems in the diffusion of matter as well 
as of those concerned with the diffusion of heat. 


HE theory of the diffusion of heat rests upon the assumption that the 
quantity 6Q of heat passing through a plane of area dy 62 in time 6¢ is 
given by the relation: 
6) a9 
— = — k—Syisz, (1) 
ét Ox 
where 06/0x is the temperature gradient in the direction of the X-axis, and k 
is the thermal conductivity. When & is taken as a constant, this assumption 
leads to the well-known differential equation, 


00 (970 06 90°6 
—=h°* ~wLrerna an? (2) 
at Ox? dy? dz* 





in which the diffusivity #2 is the ratio k/cd, c being the specific heat and d 
the density. If h? is taken as a constant, Eq. (2) may be solved for many 
particular cases by well-known methods. 

If k is not a constant, but a function of 0, Eq. (2) no longer applies. 
The differential equations that then apply for linear and radial flow have 
been given by many authors; and the general form is easily obtained in a 
manner analogous to the usual derivation of Eq. (2). The flow into and out 
of a differential block through the two faces normal to the direction of the 
X-axis is evidently, from Eq. (1), given by 
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5Q 1 Ok 0 1 06 
= -(a-> 8x) —(0—— bx) by 
5T (x) 2 dx /dx 2 Ox 


1 Ok 0 1 06 
+(a+ sx) (6 + —éx) by 6 
2 dx /dx 2 Ox 


070) «Ok 00 
=| k—+— — Jp by 62. 
Ox? Ox Ox 


a 


Similar equations apply to the flow into and out of the other two pairs 
of faces, and the total increment of heat in the differential block is given by 
the addition of the increments through the three pairs of faces. As the 
temperature @ is defined by Q/cd. 6x dy 6z, and if k is a function of 6 only, 
there is finally obtained: 


00 (0°60 3076 0°70 1 dk ( /00\? 00\? 06\? 
j—=(—+— 4" “)4— “1 (S) a(S) (“Vb 

Ot \dx*? dy Az? k dé \dx dy 02 
where f(@) is written for 1/h® or cd/k, as c and d may also be functions of @. 
Eq. (3) is the general form of the equation for heat conduction, and reduces 
to Eq. (2) for the special case where & and h* are constants. As is well known, 
the same equations apply to the diffusion of gases and liquids when @ 
represents the concentration, k the diffusion constant, and cd the concen- 
tration-pressure ratio. The material presented below applies therefore to 
problems in the diffusion of matter as well as to those dealing with the diffu- 
sion of heat. 

Due to its non-linear character, Eq. (3) cannot be solved by the methods 
applicable to Eq. (2). In papers! dealing with the experimental determination 
of k and h? will be found discussions of Eq. (3), usually for the case of linear 
flow [00/dy=00/dz=0]. The usual treatment is to employ the solutions of 
Eq. (2) with approximate corrections applicable to the particular problem. 
This procedure is facilitated by the use of formal solutions to Eq. (3) such 
as those obtained by Glage' and by Kirchhoff and Hansemann,? among 
others. 

Problems in the interdiffusion of gases have been similarly treated. 
Boltzmann?’ in particular developed formal solutions to Eq. (3) applicable 
to a number of cases. In the case of the semi-infinite solid, he showed that 
the solution to Eq. (3) could be expressed formally in terms of the quantity 
x*/t only. Advantage has been taken of this relation in the treatment of 
special cases by Wiedeburgt and Germond.* The procedure given below for 
the solution of the semi-infinite solid rests upon this same fact: the method 
will be first discussed in more general terms. 








1 Glage, Ann. d. Physik (4) 18, 904 (1905). This paper gives a critical resume of the 
classical work in this field. 

2 Kirchhoff and Hansemann, Wied. Ann. 9, 1 (1880). 

3 Boltzmann, Wien. Ber. 86, 63 (1882); Ann. d. Physik 53, 959 (1894). 

4 Wiedeburg, Ann. d. Physik 41, 675 (1890). 

5 Germond, Bull. Am. Math. Soc. 35, 460 (1929). 
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II. The discussion will be limited to the solution of Eq. (3) subject to 
boundary conditions of the type 
6=c, when F,(x,y,2,/)=0. (4) 
6=ce when F,2(x,y,2,t)=0. 


If @ is so defined, it may be proved that if there exists a function 
® (x, y, 2, 4) meeting the following conditions: 
@=C, when F,(x,y¥,2,4)=0, 
@=C, when F,(x,y,2,t)=0, 


Ob\? /db\? O?\? eh eh Ah 
Be memes 
Od Ox dy dz _ Ox? dy? dz" 6) 


where f;(®) and f.(®) are functions of ® only, then @ is a function of ® only, 
and, 


(5) 























oe (7) 
° / 
® 


dé d(log k)/ d0\? 
(SO) — fal ®)) 7 — ful @)— (—) = f<®)- : 


dé d® 

For the function 6 obtained by integration of Eq. (7) subject to the bound- 
ary conditions, 9=c,; when ®=(C,, 0=co when ®=Cz, satisfies Eqs. (4) and 
also Eq. (3), as the latter, when @ is taken as a function of ®, reduces to Eq. 
(7) by virtue of Eqs. (6). Hence @ is the required solution, and the theorem is 
proved. 

In what follows, a function satisfying Eqs. (6) for particular boundary 
conditions of the type of Eqs. (5) will be termed a resolving function for those 
boundary conditions. 

III. Suppose now that a function @’ has been determined in the above 
manner for a case of Eq. (3) in which f(6’) is any function whatever of 0’ 
but in which & is a constant, and the second term to the right of Eq. (3) is 
therefore zero. Then, 


a7’ 3°9’ 30’ 
a6’ Ax? dy? dat 


at Ce) 


This is formally identical with the second of Eqs. (6). Substitution for 
® in terms of @ in the first of Eqs. (6) yields: 


a6’\? _/90’\2 /.80’\2\ d® 
Cay) 
36" ax dy az / / do’ 


ot fil ®) 


But as ® is a function of 0’ only, fi(®)/(d®/d6’) is a function of 6’ only. 
Hence this equation is formally identical with the first of Eqs. (6). There- 
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fore 0’ is itself a resolving function for the particular boundary conditions 
involved. 

It follows that if any resolving function exists for particular boundary 
conditions, then any solution of Eq. (3) (i.e., 0’) when & is constant is, what- 
ever the form of f(@’), itself a resolving function for these particular boundary 
conditions. 

Applying this theorem to the special case f(@’) =1, it follows that if any 
resolving function exists, the solution to this special case is itself a resolving 
function. But the solutions to this special case are precisely the well-known 
solutions to Eq. (2) (4? =1). Such solutions satisfy the second of Eqs. (6); it is 
a simple matter to determine by numerical trial if they satisfy the first of 
these equations. If they do, values of f:(6’) vs. 0’ may be computed, f2(6’) = 1, 
and a solution to Eq. (3) may be obtained by integration of Eq. (7). If 
they do not, it is proved that no resolving function exists, and the method 
fails. 

The criteria for the use of the method are then that the boundary con- 
ditions be of the form of Eqs. (4) and that the known solution 6’ for the case 
k=constant, f(@’) =1 satisfies both of Eqs. (6). 

IV. By actual trial it has been found (in the manner described above) 
that the method is not applicable to many important cases, such as the 
diffusion in a slab or cylinder of infinite length when the temperature at the 
surface is kept constant. It is, however, applicable to the important case of 
the semi-infinite solid. If in such a solid the initial temperature is 69, while 
the temperature at the surface x=0 is always 6;, then Eq. (3) becomes: 


1 a0 00, d(log k) Ey ss 
at ax? do \ax/’ 
The problem is then to solve Eq. (8) subject to the conditions, 
6=6 when t=0, 
(9) 
6=0, when x=0. 
For this problem a resolving function is given by, 
ba esl (10) 
Then from Eqs. (6) and (7), @ is given by integration of 
6) log ® dé d(log k)/ d0\? d*0 
(Co g +1) + (log (=) _ (11) 
2 d® dé d® d®? 
subject to the conditions, 
6=6 when #=0, 
(12) 


6=6; when #=1. 


This integration may be accomplished for any particular case by any of 
the well-known methods of successive approximations. 
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This result is of some importance in the solution of physical problems. 
The most common type of problem in the diffusion of heat involves an initial 
constant temperature and thereafter the maintenance of another constant 
temperature on some surface. For all such cases the result for the semi- 
infinite solid serves as a first approximation, applicable within a high degree 
of accuracy for an initial time interval sufficiently small. Thus for diffusion 
in a slab, cylinder, or sphere, the temperature will vary with time and 
distance from the surface in the same way as for the semi-infinite solid from 
zero time up to that which the temperature at the center varies by more 
than the required order of accuracy from its initial value, with the further 
proviso for the sphere and cylinder that the portion in which the temperature 
maintains substantially its initial value extend out to a radius beyond which 
the curvature may be neglected with sufficient accuracy. 

_V. As an actual example of the method discussed in the preceding section 
there will be considered the case of the cooling of a large steel ingot, a problem 
for which the simpler case is discussed by Ingersoll and Zobel, p. 85.° Strictly 
comparable values of the relations between diffusivity and temperature and 
between conductivity and temperature were not found in the literature, 
but data were taken from the Landolt-Bérnstein Tables which were thought 
to be sufficiently representative of a low carbon steel. The diffusivity was 
taken as that given for Bessemer steel from the data of Kirchhoff and Hanse- 
mann,’ giving the diffusivity as: 1/f(@) =h?(1+66), where h?=0.1148 and 
8= —0.00019. The conductivity was taken from the mean values given in 
the Tables’ for a 0.64 percent carbon steel, which gave, for the usual approx- 
imation of a linear relation: k=ko(1+a@), where ky = 0.115 and a= —0.00046. 

Thus, from Eqs. (11) and (12), the temperature in the block is given by 


ntegration of: 





a0 log ® 1 d0 «a /d0\? 
a. a’: 
de? 2n2(1+60) /& db 1+a0\de 


subject to the boundary conditions which were that 6=800 when @=0(¢=0), 
and @=20 when ®=1(x=0). 

The method of solution employed was that of successive approximations 
in which the integrations are mechanically obtained by means of an inte- 
graph, as described in a paper by Fry.* The work was simplified by taking 
as a first approximation the well known solution to the problem when 
h? and k are constant, which may be written in terms of ® as 


2x780 ¢* 2 
6=800—-— f edn. (14) 


gil? (log? /2h) 





In this expression h was taken as having its mean value for the temper- 
ature range involved or (0.1060)'”*._ Due to the discontinuity at ®=0 and 


® Ingersoll and Zobel, Mathematical Theory of Heat Conduction. 

7 Kirchhoff and Hansemann, Wied. Ann. 13, p. 406 (1881). 

8 Landolt-Bérnstein Tabellen, 5th Ed. p. 1291 (1923). 

9 T.C. Fry, Proc. Int. Math. Congress, Vol. II, p. 405, Toronto (1924). 

















—— 











GENERAL EQUATION OF DIFFUSION 559 


various associated peculiarities of the equation, some slight modifications 
in procedure were necessary. Each approximation was made as follows: 
Given approximations for @ and its first derivative, say 0, and 0,’ respectively, 
the corresponding second derivative 6@,’’ was computed from Eq. (13). 
This was then integrated to give 0’,.:, the constant of integration being taken 
such as to give 0’,,; the same value at ®=1 as that found at this point for 
6,,’. In general, this did not make @’,,; approach zero as ® approaches zero, 
a condition which must be satisfied, as is evident from the fact that a 
function of the type of Eq. (14) is a close approximation to a solution for 
low values of ®. Hence @’,,; was corrected to meet this requirement by 
taking values of 6 proportional to the second term of Eq. (14) up to such a 
value of ® as would make this approximation merge into the original curve 
for 6’,,; without any abrupt change in curvature. The @’,,, thus corrected 
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was integrated to give @,,;, the constant of integration being such as to make 
6,41:=800 at @=0. The difference between 20 and the value of 6,4: at 
#=1, or C, was taken as a correction factor, and a term ®C added to 
6,4: to make it meet both boundary conditions, while the term C was 
also added to 6’,,;, in accordance with the procedure developed by Fry.°® 
Then @’,,; was again corrected to make it approach zero as ® approached 
zero, by the same method as described above. The values of 6,4: and 
6,41, thus finally obtained were used to compute 6’’,,; and the operation 
repeated until successive approximations produced no change in @ or 6’. 

The solution obtained for the case cited is shown graphically in Fig. 1, 
in which both @ and @’ are given. The corresponding solution and its first 
derivative for the case in which the conducitivity and diffusivity are con- 
stant (the latter taken as h?=0.1060), as given by Eq. (14), are shown in 
broken lines. It will be seen that the solutions to the two cases differ con- 
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siderably. To emphasize this difference, Fig. 2 has been prepared, cor- 
responding to Fig. 12, p. 86 of Ingersoll and Zobel,® and giving the rate of 
cooling (00/0/) at distances of 0.3 cm and 1.0 cm from the face of the block. 
The broken lines again show the corresponding curves for constant diffusivity 
and conductivity. It is apparent that consideration of the variation in the 
conductivity and diffusivity results in a higher maximum rate of cooling 
than if this variation is neglected. 
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Fig. 2. 


VI. Another case in which the method applies is that of a cylinder of 
infinite length and radius, initially at one constant temperature, of which 
the axis is thereafter kept at another constant temperature. If x is the 
radial distance of any point from the axis, the problem is to solve: 





10) a0 1 a0 8°  d(log k) (“) an 
at ox Ox Ox? dé dx} ’ 
subject to the conditions, 
6=0) when (t=0, 
‘ (16) 


6=6, when x«=0. 


For this case also the ® of Eq. (10) is a resolving function. From Eqs. (6) 
and (7), @ is given by integration of 


f(9) log ® 1 ) dé d(log k)/ 00 \? a*6 
1 )}— p————_[|_ ——_- Pp =0. 17 
( 2 "in .* d® * d@ (3) * dP? (17) 








subject to the conditions of Eqs. (12). 

Brown and Furnas! have discussed the determination of a variable diffus- 
ivity by the graphical integration of Eq. (15) for numerical values of the 
differential terms of this equation in the case of a centrally heated cylinder. 


10 Brown and Furnas, Trans. Am. Inst. Chem. Engineers 18, 295 (1926). 
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The above development suggests a procedure in which the measurements 
and computation would be simpler, though the experimental conditions would 
be more difficult to obtain, as to maintain a constant temperature along the 
axis would require some method of adjusting the current in the central heat- 
ing coil in such a way as to fulfill this condition. It would also be necessary 
to limit the measurements to a temperature range such that the surface never 
varied appreciably from its initial value. If these conditions are fulfilled 
Eq. (17) applies, and (assuming k known as a function of 6*), it is only necess- 
ary to take readings from one thermocouple (one value of x), as @ is a function 


. of ® only, and any value of ® will be obtained at any x for some value of ¢. 


Expressing the values of ¢ in terms of ®, values of @ vs. ® are experimentally 
obtained. Then by graphical differentiation, or numerically by successive 
differences, values of d0/d® and d*6/d®? may be obtained, and Eq. (17) 
solved for f(@), the reciprocal diffusivity, as a function of temperature. 

VII. The development discussed above permits the reduction of the 
general differential equation of diffusion to a single ordinary differential 
equation of the second order, which can be readily solved by methods of 
approximation. The only particular case of major importance to which the 
method applies is that of the semi-infinite solid, which affords, however, a 
good approximation to a solution for many other cases within limited ranges 
of the co-ordinates of time and distance. The form of the differential equation 
resulting when the method applies is such as to make the solution by the 
method of successive approximations preferable, and for this graphical 
methods are more convenient than numerical. 

As is apparent in the example given above, solutions obtained by this 
method when the variations in conductivity and diffusivity are allowed for 
may differ considerably from those obtained on the assumption that these 
quantities are constant. Such differences may be particularly apparent in the 
values obtained for rates of cooling or of temperature gradients, as the 
example again indicates, and it is in the evaluation of these quantities 
that most accurate determinations are desirable (e.g., in the casting of metals 
and glass). It must, however, be admitted that conditions in technical 
operations seldom approach those of the mathematically ideal problem 
sufficiently to make such accuracy in calculation significant, and the 
developments discussed above are more likely to be of practical importance 
in problems concerned with material diffusion than in those involving heat. 

Acknowledgment is due Dr. T. C. Fry of these Laboratories for much 
helpful criticism of the work reported in this paper. 


* This relation would be previously determined by measurements for the steady state, 
as discussed by_Brown and Furnas (reference 10). 
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for the second issue, the thirteenth of the month. The Board of Editors does not 
hold itself responsible for the opinions expressed by the correspondents. 


On the Theory of Electrons and Protons 


In a recent paper,' Dirac has suggested that 
the reason why the transitions of an electron 
to states of negative energy, which are pre- 
dicted by his theory of the electron, do not 
in fact occur is that nearly all of the states of 
negative energy are already occupied. Dirac 
has further shown that the unoccupied states 
of negative energy have many of the proper- 
ties of protons; that, for instance, they may 
be represented by wave functions which 
would be taken to correspond to a particle 
of positive charge and positive mass. He has 
further shown that the mass associated with 
these gaps is not necessarily the same as that 
of the electron, and he has suggested the 
assumption that the gaps are protons. In 
order to account for the fact that the diver- 
gence of the electric field is not, in spite of the 
infinite electron density, everywhere infinite, 
Dirac further assumes that only the depar- 
tures from the normal state in which all 
negative states are filled are to be counted in 
computing the charge density for Maxwell's 
fourth equation 

div E= —4zp. (1) 
Finally, Dirac is able to account for the 
validity of the Thomson formula for the 
scattering of soft light by a free electron, in 
spite of the fact that the derivation of this 
formula on his theory of the electron—a 
derivation which makes explicit use of the 
transitions to states of negative energy which 
are now forbidden, is invalid. According to 
Dirac, the scattering takes place by a double 
electron jump, in which a negative? electron 
jumps up to some state of positive energy, 


1P, Dirac, Roy. Soc. Proc. A126, 360 
(1930). 

2 By a negative electron we mean an elec- 
tron of negative energy. 


and the original positive electron falls down 
into the gap left. 

There are several grave difficulties which 
arise when one tries to maintain the sugges- 
tion that the protons are the gaps of negative 
energy, and that there are no distinctive 
particles of positive charge. In the first place, 
we can easily see that Dirac’s theory requires 
an infinite density of positive electricity; 
and since we should expect the de Broglie 
waves of this charge to be quantized, we 
should expect some corpuscular properties 
for the positive charges. The reason why the 
theory requires an infinite positive charge is 
this: If the explanation of the scattering of an 
electron is to be tenable, a negative electron 
must interact with the electromagnetic field 
in the way predicted by Dirac’s theory of the 
electron; for otherwise the scheme proposed 
would not give the Thomson formula. But 
this means that there must be a term involv- 
ing the current and charge vector of the nega- 
tive electrons in the total energy momentum 
tensor for matter and radiation. Thus by (1), 
the divergence of the electric field will be 
everywhere infinite unless there is an infinite 
density of positive electricity to compensate 
the negative electrons. 

A further difficulty appears when we try 
to compute the scattering of soft light by a 
proton. This difficulty is not unconnected 
with the difference in mass between the elec- 
tron and proton, and makes it seem im- 
probable that this difference can be explained 
on the basis suggested by Dirac. For the 
scattering process must in this case be re- 
garded as a double jump of a single electron, 
in which a negative electron jumps to some 
state of positive energy, and then falls back 
into the hole that is the original proton. Now 
it is easy to see that the probability of this 
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scattering is determined by precisely the 
same matrix components as those which give 
the electron scattering, and that the present 
theory gives equal scattering coefficients for 
electron and proton. Of course, the interac- 
tion between electrons is omitted in this 
computation; but the difficulty is this, that 
such interaction would affect electron scatter- 
ing and proton scattering in precisely the 
same way; whereas the Thomson formula 
requires the latter to be smaller by a factor 
proportional to the square of the ratio of the 
masses, 

Finally, there is a numerical discrepancy 
to be noted. According to Dirac’s suggestions, 
the filling of the proton gaps in the distribu- 
tion of negative electrons should correspond 
to the annihilation of an electron and a pro- 
ton, and should thus, under all normal con- 
ditions, be a very rare occurrence. Now if we 
consider for definiteness a free electron in an 
enclosure in which there are m, free protons 
per unit volume, we may readily compute the 
rate at which the electron should, by the 
Dirac radiation theory, fall into one of the 
corresponding gaps. The conservation laws 
require that at least two quanta be emitted 
in this process; and it is sufficient to consider 
jumps in which no more than two quanta are 
emitted. The details of the calculation will 
be published elsewhere; if we neglect the 
interaction of the electron with the negative 


electrons we obtain for the mean life time of 
the electron: 

T =Gm'ce—4*/n, (2) 
where G is a numerical constant of the order 
of unity, e the charge, and m the mass of the 
electron. Now again it is difficult to see what 
large errors could be involved in the compu- 
tation, since the matrix components which 
give (2) are of precisely the same type as 
those which give correctly the Thomson 
formula and the optical transition probabili- 
ties of the electron; and (2) gives a mean life 
time for ordinary matter of the order of 10- 
seconds. 

Thus we should hardly expect any states of 
negative energy to remain empty. If we re- 
turn to the assumption of two independent 
elementary particles, of opposite charge and 
dissimilar mass, we can resolve all the diffi- 
culties raised in this note, and retain the 
hypothesis that the reason why no transi- 
tions to states of negative energy occur, either 
for electrons or protons, is that all such states 
are filled. In this way, we may accept Dirac’s 
reconciliation of the absence of these tran- 
sitions with the validity of the scattering for- 
mulae. 

J. R. OpPENHEIMER 

The Norman Bridge Laboratory of Physics, 

California Institute of Technology, 
Pasadena, California, 
February 14, 1930. 


The Collision Diameter of the Hydrogen Atom 


Harteck measured the viscosity of mona- 
tomic hydrogen and calculated the collision 
radius to be ry =1.3 X1078 cm (ro of the first 
Bohr orbit is 0.5310-§ cm). We have 
measured the decrease of intensity of a fine 
beam of monatomic hydrogen in passing 
through 3.0 cm of mercury vapor at 0.000185 
mm Hg and find the decrease of intensity to 
be 23 percent. If the collisions are supposed 
to be between elastic spheres the sum of their 
radii is 6.2 X 10-8 or assuming ry, = 1.8 X 1078, 


Excited Radicals in 


The letter of Mr. D. H. Andrews in your 
number of December 15th makes the im- 
portant point it is possible to assign definite 
frequencies to vibrating groups in the mole- 
cule, and that these frequencies are constant, 
to a first approximation, from molecule to 
molecule e.g. the value of the vibration fre- 


the radius of hydrogen atom is 4.4X10~*. 
This is larger than would be assumed even 
on wave mechanics unless there is very con- 
siderable resonance interaction between mer- 
cury and hydrogen. 

E. G. LUNN 

F. R. BicHowsky 


U.S. Naval Research Laboratory, 
Anacostia, D. C., 
February 17, 1930. 


Chemical Compounds 


quency of C-H in any molecule will be of the 
order of 3000 cm™, altering only slightly with 
change of environment. Although some of 
Mr. Andrews’ values are open to question, 
there is a weight of evidence, both from the 
Raman effect and from infra-red spectra, for 
the essential correctness of his statement; 
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the failure to realise it has been the cause of 
false analyses of the spectra of several 
polyatomic molecules, in which the frequency 
of, say, the CO group has been given values 
varing from 1500 to 200 cm™. 

It is possible, however, to go much further. 
First, the frequency of a group, which remain- 
ing constant in molecules where the chemical 
environment is the same, can undergo a 
sharp change to a second value corresponding 
to a second chemical environment. I was led 
to this conception from a study of the 
building-up of the N2O molecule from NO. 
The analysis of the vibration-rotation band 
system of N2O (to be published shortly) 
shows that the molecule is linear and sym- 
metrical. The vibration of N-O along the 
axis of the molecule is 2250 cm™. Now the 
frequency of the NO molecule in its normal 
(Il) state is 1892 cm; the building of the 
N:0 molecule has brought about a very 
definite change. A name can be given to this 
change when it is realized that the frequency 
of the first excited (22) state of NO is 2345 
cm. The numerical resemblance is striking; 
and it becomes more than a coincidence when 
we find that the only other excited state (*I1) 
so far known has a frequency of 1030 cm™, 
and that the Raman effect and infra-red 
spectra of compounds containing the NO 
group give a value c. 1000 without exception. 
With these frequencies it is possible to fit all 
the results upon the NO group at present 
known i.e. a diatomic group within a poly- 
atomic molecule can be described in terms 
of the excited states of a diatomic molecule. 

The results are even more definite when the 
CO group is considered. In a normal ke:one 
or acid the vibration frequency is c. 1600 
cm, The Raman evidence is contradictory 
but the infra-red quite convincing. On the 
other hand, an a-diketone of the form 
R.CO 

| seems, from the careful analysis of the 
R.CO 
absorption spectra of several compounds of 
this type by Louis Light of Ziirich, to be of 
the order of 1100 cm™. 

The first three states of the CO molecule 
are: 


& iz 2155 
a 511? 1724 
d $y? 1105 


the agreement is satisfactory. The effect of 
the conjugated system is particularly inter- 
esting. 

Certain principles may be set down, even 
at this early stage: 


(i) A diatomic group persists throughout a 
large number of molecules with the same 
vibration frequency. As Mr. Andrews sug- 
gests, this establishes the identity of covalent 
bindings. The mass of the atoms to which 
the diatomic groups are attached is of small 
importance. 

(ii) A diatomic group can have several 
sharply distinct frequencies. Each of these 
frequencies behaves as defined in (i). 

(iii) The frequencies of a diatomic group 
correspond to the frequencies of excited 
states of the appropriate diatomic molecule. 
A frequency corresponding to the frequency 
of the normal state of the molecule is not 
usually observed. 


At the moment the attempt to put a 
spectroscopic name to chemical binding must 
stop with groups bound by covalencies. 
Clearly the extension of the idea to such prob- 
lems as the boron hydrides is practicable; 
if in the B.H, molecule there are B—H 
groups bound by singlet linkages, then the 
frequencies of these groups should not 
coincide with the frequency of any of the 
states of the ordinary BH molecule bound by 
a covalency. 

Similarly, the kind of linkage shown to 
exist by Lennard Jones for the normal ('2) 
state of the oxygen molecule should be respon- 
sible for abnormalities in the way outlined 
above. There are indications that the C—C 
binding in C;H, is such a case; if that is 
true, it will be possible to indulge in a 
spectroscopic re-naming of the “‘triple’’ bond. 

C. P. SNow 

Laboratory of Physical Chemistry, 

Cambridge, England, 
January 27, 1930. 


The Glancing Angle of Reflection from Calcite for Silver (Ka,) X-rays 


We have measured the first order glancing 
angle at which the Ka, line of silver is re- 
flected from the cleavage planes of calcite 


by a method independent of any used by 
those who have obtained the most reliable 
results so far. The weighted mean value of 
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fourteen independent observations reduced to 
18°C is 5°17'13”.81 +0. 06. 

This value was obtained from spectrograms 
taken on a new spectrometer specially de- 
signed to utilize the displacement method 
introduced by H. S. Uhler (Phys. Rev. [2] 
11, 1-20 (1918)). A full description of the 
method, instrument and results will be given 
later. 

The value of the wave-length correspond- 
ing to the above glancing angle is 0.558238A + 
0.000002A, obtained by using 3.02904A for 
the “effective” grating space of calcite 
(Siegbahn; Spectroscopy of X-Rays, 1924 
Edition). 

The thermal coefficient of expansion at 
18°C for the grating space of calcite obtained 
from the most reliable sources, without 
approximations, is 0.0000102;/°C. This 
value differs by about 2% from the generally 
quoted value originally given by W. Sten- 
strém (Dissertation, Lund, 1919) as a suf- 
ficient approximation. 

It is interesting to note that our value of 
the glancing angle is 1”.4 smaller than that 
given by G. Kellstrém (Zeits. f. Physik 41, 
516 (1927)) for the first order of calcite,viz. 


5°17'15".2. The probable error of the un- 
weighted mean of Kellstrém’s seven deter- 
minations for this order is +0."04,. The 
unweighted mean of his four second order 
determinations is also larger than ours and 
has relatively four times this probable error. 
With such a difference existing one wonders 
if the grating spaces of the crystal specimens 
used were the same. We employed an unu- 
sually selective piece of calcite which was 
cleaved from an excellent specimen of Iceland 
spar contained in the Marsh collection of 
Yale University. 

Our value is larger than those of K. Lang 
(Ann. d. Physik 75, 489 (1924)), A Leide 
(Dissertation, Lund (1925) and A. P. Weber 
(Zeits. f. Wiss. Photo. 23, 149 (1925)) 
although coming within 0”.11 of Lang’s 
and 0”.36 of Leide’s results when their 
probable errors are added to their respective 
mean values, 


CHARLTON Dows CooKsEY 
DonaLp CooKsEY 
Sloan Physics Laboratory, 
Yale University, 
February 11, 1930. 
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BOOK REVIEWS 


Quantum Mechanics. E. U. Convon anv P. M. Morse. Pp. 250+-xiii, figs. 28. International 
Series in Physics, McGraw-Hill Book Co., New York, 1929. Price $3.00. 

This is the first book of American authorship on the new science of quantum mechanics 
and should find a host of readers both here and abroad. It meets the widely felt need for an 
introductory text making no pretense at completeness, but aiming through numerous illustra- 
tive applications and exercises to give the reader a working knowledge of the methods of the 
quantum theory of today. The authors have kept the work, to use their own words, “on as 
elementary a plane as possible,” but they have not descended to sloppy procedure in their 
anxiety to make easy reading. 

The subject is approached from the standpoint of the statistical interpretation of the 
Schrédinger theory and the first three quarters of the volume is devoted to an exposition of 
results obtainable from the Schrédinger wave equation with a minimum of mathematical 
machinery. A subsequent chapter, however, gives a well written introductory account of the 
general formulation of the theory by Jordan, Dirac, and von Neuman, which should be helpful 
to many who have had difficulty in digesting the original papers on this subject. 

The rather extensive treatment of diatomic molecules will be most useful to workers in the 
field of band spectra as there is at present no comprehensive work on the quantum mechanics 
of molecules. Much of the material in this chapter is derived from the valuable researches of 
the authors. The reviewer notes with especial pleasure the inclusion of a simplified account 
of the basic molecular perturbation theory of Born and Oppenheimer worked out directly for 
the important diatomic case. 

The book also contains an interesting chapter on aperiodic phenomena. 

The authors lay little emphasis on the philosophical aspect of the quantum mechanics and 
ignore almost completely many important phases of the subject such as the “electron spin” 
and the problem of dispersion. These omissions are to be reckoned a virtue rather than a 
drawback, however, since they enable the authors to deal with the fundamentals of the subject 
in a leisurely way which will be welcome to the beginner. 

Epwin C, KeMBLE 


Outlines of Biochemistry. R. A. GortNer. Pp. xv and 793. New York, Wiley and Sons; 
London, Chapman and Hall; 1929. Price $6.00. 

Treatises in biochemistry are the successors of medical textbooks of physiological chem- 
istry; these dealt almost exclusively with mammalian metabolism, foods, enzymes, excreta 
and the blood; and although when vitamins were recognized and the part played by the secre- 
tions of the endocrine glands, these matters were included, the chemistry of the life processes 
of plants was always beyond their scope. Even with such restrictions on the subject matter, 
the “conservation of energy” was referred to when dealing with the measurement of animal 
heat, a formula or two of chemical kinetics got in to ornament the chapter on digestion, and 
mention was made of the freezing-points of serums and of urine; but—with an eye, no doubt, 
to the limitations of the medical student—such topics received in general but incidental men- 
tion. 

It is a far cry from such a book to Professor Gortner’s “Outlines.” He has a section on the 
proteins of course, including chapters on the amino-acids and the polypeptides, an up-to-date 
account of the vitamins written by Professor L. S. Palmer, chapters on the carbohydrates and 
the fats and oils, and others on the tannins and on the colouring matters of plants; but the 
first 290 pages of his book are devoted to “The colloidal state of matter.” Here the preparation 
and characteristic properties of colloids are described, their electro-kinetic behaviour, and their 
behaviour with electrolytes and among themselves; there is a chapter on gels, and one each on 
hydrogen-ion-concentration, surface-tension, osmotic-pressure, and “the Donnan equilibrium.” 
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It is obvious that in this book the generalizations of physical chemistry are not slipped in as 
“explanations” of an isolated fact or two, but have the place of honour; in their light the facts 
of physiology are examined—not mammalian alone, for as might be expected from a professor 
of agricultural biochemistry the plant receives its due. 

To readers of The Physical Review the illustrations here afforded of the scope and power 
of generalizations long familiar should prove intensely interesting. Who that measures vapour- 
tensions will not enjoy the picture of the cactus on the wall in desert Arizona ten months with- 
out rain and sheltered from the dew, or the story of the three-months sojourn of another cactus 
over sulphuric acid, the account of drought-resisting cereals in Utah, winter-hardy wheat, 
and the drying-up of earthworms, or the human-interest paragraphs on the granary-weevils’ 
life without a drink, the new oedema diagnosis, and the rigor of our bodies after death? All 
these and more are brought into the chapter on Gels to illustrate “imbibition” and the 
distinction between “bound” and “liquid” water; and not as incidental illustrations merely, 
for they are examined in the light of tables, formulas and charts, and correlated with freezing- 
points, osmotic-pressures and absorption in the infra-red. Therefore, if any Reader thinks 
he knows more physics though less biology than Gortner, let him go to it; the door is opened 
wide, 

The sample given above could easily be matched; literature citations show where synoptic 
papers can be found; it would be a shame to let mere biochemists monopolize the reading of this 
book. W. LasH MILLER 


Introduction to Physical Optics. JoHN KeLttock Ropertson. Pp. 422, figs 222. D. Van 
Nostrand Co., New York 1929. Price $4.00. 

This is a text-book which is intermediate in grade between the elementary and the ad- 
vanced optics texts. A small amount of differential calculus is used, but otherwise the mathe- 
matics required consists of the usual trigonometry, algebra, and geometry. A large number of 
problems are given, the answers to which are found at the end of the book. Several problems 
are solved in great detail in the text. 

A very commendable feature is the use of photographic illustrations of the various optical 
phenomena such as interference and diffraction bands, and spectra. These show the various 
types of band systems obtained with different types of apparatus, giving the student a clear 
impression of the phenomena under discussion. The photographs are all good except the one of 
interference rings in polarized light. 

In trying to be neither too conservative, nor radical, the author seems to overemphasize 
the conflict between the corpuscular and wave theories of light. Naturally, most of the phenom- 
ena treated require the wave theory. The nature of waves and the composition of vibrations is 
fully and clearly discussed in an early chapter. The principles thus established are later applied 
to the treatment of interference, diffraction, and elliptical polarization. The Cornu spiral is 
explained in a graphical, qualitative manner, and is used to explain some of the Fresnel diffrac- 
tion phenomena. 

The quantum theory and the foundations of Bohr’s theory of spectra are briefly reviewed. 
There are also a few paragraphs on wave mechanics and the diffraction of electron waves, 
which give a rather belated suggestion of the duality of radiation as well as matter. Optical 
phenomena due to motion are discussed in the closing chapter. 

JosEPH VALASEK 


Elementary Differential Equations. THornron C. Fry. Pp. x+255, figs. 42. New York, 
D. Van Nostrand Co., 1929. Price, $2.50. 

This textbook gives a pleasing treatment of its subject of the grade of difficulty that is 
adapted to junior courses in the better American engineering colleges. It is distinguished by 
careful presentation of the fundamentals, recognizing that while applications of the subject in 
physics and engineering are of the highest value, nevertheless, “They are not without their 
danger, however, if they are used unwisely, for should the student come to rely upon his 
physical intuition as a substitute for abstract logic, rather than as an aid to it, he would have 
missed half the benefit of the study.” (Preface). 
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There is a good treatment of operational methods of handling linear equations with con- 
stant coefficients, and systems of them, with emphasis on the applications to electrical circuit 
theory. E. U. Connon 


Die Grundlage der allgemeinen Relativitatstheorie. A. Einstein. Pp. 64. Fifth unaltered 
printing. J. A. Barth. Leipzig, 1929. 

This is a fifth unaltered printing of Einstein’s famous 1916 treatment of the foundations 
of the general theory of relativity. It was first published, except for the table of contents, as a 
single article in volume 49 of the Annalen der Physik, and was the first systematic treatment of 
general relativity, providing the complete mathematical and theoretical foundation leading to 
the three tests of the theory furnished by the rotation of the perihelion of Mercury, the bending 
of light in passing through the gravitational field of the sun, and the red shift for light from the 
surface of the sun. The paper is not only a document of the greatest historical interest, but 
also still furnishes one of the best introductions to general relativity and contains the three con- 
clusions of the theory for which we have the most direct experimental evidence. 

RicHArD C, ToLtMAN 





